GENERAL  DYNAMICS 

Convair  Division 


DISCLAIMER  NOTICE 


THIS  DOCUMENT  IS  THE  BEST 
QUALITY  AVAILABLE. 

COPY  FURNISHED  CONTAINED 
A  SIGNIFICANT  NUMBER  OF 
PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


gd/c -DBE65 -023 


HIGH  ALTITUDE  SOCKET  FLUME  STRUCTURE 


Final  Repo 

Off  let  of  Naval  research 
Contract  No.  Nonr  4-747(00) 
ARPA  Order  No.  562 

September  1965 


This  research  is  part  of  Project  DEFENDER 
sponsored,  by  the  Advanced  Research  Projects 
Agency,  Department  of  Defense,  and  the  Office 
of  Naval  Research. 


GENERAL  DYNAMICS 
CONVAIR  DIVISION 


Distribution  of  this  report,  is  unlimited 


FOREWORD 


This  report  describes  research  carried  out  under  Contract  No. 

Nonr  471*7(00)  during  the  period  December  1964  to  June  1965.  The 
investigation  includes  a  numerical  analysis  of  the  neutral  gas  structure 
of  ionospheric  rocket  plumes  including  diffusive  affects  and  a  principally 
analytic  study  of  the  ion  and  electron  distribution.  The  program  for  the 
numerical  evaluation  of  the  neutral  plume  structure  had  been  under 
deveicpp’.ent  previously  under  Contract  AF19(628) -43bO,  ARPA  Order  No.  363 
and  this  was  continued  concurrently  with  the  present  study.  Although  it 
is  not  really  possible  to  isolate  the  two  efforts,  the  present  study 
has  been  principally  concerned  with  the  incorporation  of  diffusive  effects 
and  with  the  evaluation  of  the  effects  of  the  geomagnetic  field  on  the 
Ion  distribution. 

Because  no  description  of  the  numerical  procedure  is  presently 
available,  w°  give,  in  this  report,  a  fairly  detailed  description  of  both 
the  inviscid  and  viscous  calculations.  Much  of  the  work  on  the  neutral 
flow  structure  will  be  reported  subsequently  under  Contracts  No. 

AF19( 628) -3269  and  AF19( 628) -4360. 

Personnel  who  have  spent  an  appreciable  fraction  of  their  time  on 
this  study  include  J.  R.  Barthel,  I.  J.  Brainerd,  R.  S.  Janda,  M.  S. 
Schoonover,  and  J.  A.  L.  Thomson.  The  Project  Scientist  is  J.  A.  L. 
Thomson.  The  technical  monitor  is  Morton  Cooper,  Fluid  Dynamics,  Office 
of  Naval  Research. 

Figures  11-1  to  II -9  are  classified  Confidential  and  are  submitted 
as  a  separate  addendum  to  this  report;  GD/C-DBE65-O2.5  (Addendum). 


ABSTRACT 


A  calculation  of  the  flow  field  in  the  vicinity  of  a 
liquid-fueled  rocket  travelling  through  the  ionosphere  is 
described.  A  finite  difference  program  is  used  to  obtain 
both  inviscld  and  viscous  approximations  to  the  exhaust  and 


air  flow  field.  The  subsonic  nose  region  is  treated  by  an 
approximate  procedure.  Viscous  effects  are  included  in  the 
so-called  merged  layer  regime.  The  viscous  effects  are 
treated  by  a  modification  of  a  procedure  used  by  Cheng  for 


the  nose  region  of  a  blunt  body  in  hypersonic,  low  Reynolds 
numoer  flow.  Numerical  calculations  of  the  neutral  flow 
field  in  the  viscous  layer  approximation  are  presented  for 
the  flow  at  180  km. 

The  influence  of  the  geomagnetic  field  on  the  motion  of 
the  ions  is  significant  at  altitudes  above  about  150  km. 

Simolified  geometries  are  used  tc  study  the  ion  motion.  At 
low  altitudes  (below  about  150  km)  the  ions  move  with  the 
neutrals.  At  high  altitudes  (above  250  km)  the  ions  are  more 
or  less  constrained  by  the  geomagnetic  field  to  a  one-dimensional 
motion  along  the  field  lines  except  in  a  region  relatively  close 
to  the  missile  where  the  density  is  sufficiently  high  tc  sweep  the 
ions  across  the  field  lines.  The  model  of  the  nigh  altitude  ion 


motion  is  essentially  that  proposed  previously  by  Lighthill 
The  mechanisms  describing  the  motion  cf  the  ions  across  ohe 


(I960). 

field 


lines  are  examined  in  various  limi4 ing  cases. 
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I.  Introduction 


At.  high  altitudes  the  static  pressure  of  the  exhaust  gases  at  the 
exit  plane  of  a  rocket  engine  nozzle  is  considerably  greater  than  that  of 
the  ambient  atmosphere ,  The  exhaust  gases  therefore  expand  greatly  upon 
leaving  the  nozzle.  This  expansion  and  the  resulting  interaction  with 
the  atmosphere  results  in  the  characteristic  high-altitude  missile  "plume". 

As  the  vehicle  and  plume  move  through  the  atmosphere >  the  ambient 
medium  is  displaced  in  much  the  same  way  that  it  would  be  by  a  solid  body 
having  the  shape  and  dimensions  of  the  plume.  The  disturbance  of  the 
external  medium  and  particularly  of  the  ambient  electrons  by  the  passage 
of  the  rocket,  is  the  principal  object  of  the  investigation  described  here. 

If  viscous  effects  were  negligible,  it  would  be  possible  to  describe 
the  external  flow  exactly  as  that  around  a  solid  body  having  the  same  size 
and  shape  as  the  plume.  However,  the  plume  size  and  shape  are  not  known 
a  priori,  and  must  be  described  by  a  method  which  treats  both  the  internal 
and  external  flow  in  considerable  detail.  The  inclusion  of  transport 
effects  makes  the  flow  near  the  exhaust -air  boundary  considerably  different 
from  the  flow  about  a  solid  body,  since  finite  tangential  velocities  at 
the  boundary  and  diffusion  across  the  boundary  must  be  considered  in  the 
case  of  the  plume  at  high  altitudes. 

The  altitude  range  of  interest  in  the  present  study  is  100  to  300  KM. 
The  following  generalizations  about  the  exhaust -ambient  flow  field  apply 
in  this  altitude  range: 

1)  The  free  stream  is  supersonic  (Mi  t),  partially  dissociated,  and 
weakly  ionized.  , 
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2)  The  flow  field  is  time -dependent  due  to  atmospheric  strc tification 
and  vehicle  acceleration. 

3)  The  flow  field  may  not  be  axially  symmetric  for  non-vortical 
trajectories  due  to  atmospheric  stratification. 

4)  The  Reynolds  and  Knudsen  numbers  lie  in  ranges  where  viscous  and 
mixing  effects  are  Important, ,  especially  in  the  nose  region  and 
along  the  plume  boundary . 

5)  The  ion-neutrax  and  ion- ion  mean  free  paths  are  large  enough  that 
the  geomagnetic  field  may  strongly  perturb  the  ion  and  electron 
density  profiles. 

6)  The  mean  free  path  is  large  compared  wivh  missile  dimensions, 
f "  that  the  relative  disturbance  of  the  free  stream  by  the 
missile  itself  is  small. 

The  approach  we  use  to  determine  the  neutral  flow  field  is  discussed 
in  detail  in  Section  IV.  In  general  the  approach  is  to  use  finite  difference 
procedures  to  obtain  numerical  solutions  of  the  appropriate  ir.viscid  and 
viscous  equations  of  motion.  Inviscid  calculations  of  the  air  and  Jet  flow 
are  carried  out  first.  Some  features  of  the  plume,  such  as  the  overall  size 
and  shape,  appear  to  be  quite  well  described  by  the  inviscid  equations. 
Significant  viscous  effects  are  confined  more  to  the  detailed  form  of  the 
density  and  temperature  distributions.  A  brief  preliminary  discussion  of  the 
magnitude  and  importance  of  the  viscous  effects  is  given  in  Section  II-4. 

Since  the  cyclotron  frequencies  of  the  ions  and  electrons  in  earth's 
magnetic  field  may  be  large  compared  to  the  collision  frequencies,  the 
geomagnetic  field  is  expected  to  have  a  significant  influence  on  the  ion 
motion.  The  ionization  level,  however,  is  sufficiently  low  that  the  motion 
of  the  neutral  particles  is  essentially  unaffected  by  the  presence  of  the 
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ions.  The  neutral  flow  field,  therefore,  may  be  evaluated  without  regard 
to  the  magnetic  field.  The  ion  motion  may  be  evaluated  subsequently! 
treating  the  rrotion  of  the  neutral  gas  as  known.  Most  oi*  the  analysis  in 
the  present  study  is  directed  towards  establishing  a  model  for  the  ion 
motion  that,  may  he  used  to  evaluate  the  electron  distribution  in  the 
neighborhood  of  the  vehicle.  The  major  uncertainty  lies  in  the  deter¬ 
mination  of  the  ion  motion  across  the  mag:  atic  field  lines.  Inhomogene ities 
in  the  ambient  ionosphere,  asymmetrical  ion  and  neutral  flow  distributions, 
the  finite  time  dependence  of  the  neutral  flow  may  all  enhance  the  rate  of 
diffusion  of  ions  across  the  field.  The  motion  of  the  .Ions  is  studied  in 
a  number  of  limiting  cases  in  Section  V.  A  summary  description  of 
the  ion  motion  is  given  in  Section  HI  and  ox*  the  neutral  motion  in 


Section  II. 
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II.  SUMMARY  CF  Ifl&JSNT  WETHOD  USED  FOR  CALCUIATING  THE  NEUTRAL  FLCW 


II-l.  General  Description 

The  high  altitude  plume  may  be  divided  into  three  regions:  the 
undisturbed  air  and  jet  streaa  and  the  interaction  region  (see  FigureX"l*«  The 
inner  boundary  of  the  interaction  region  represents  the  greatest  penetra¬ 
tion  into  the  plume  of  the  influence  of  the  external  environment.  This 
inner  limit  of  the  interaction  region  is  the  Jet  shock.  Inside  of  this 
shock  the  Jet  flow  is  identical  to  the  flow  into  a  vacuum.  In  the  inter¬ 
action  region  the  exhaust  gas  flow  is  separated  from  the  perturbed  air 
flow  by  a  slip  surface  (when  diffusion  is  negligible).  The  outer  boundary 
of  the  interaction  region  is  the  air  sho.'k.  At  high  altitudes,  where  much 
of  the  Interaction  region  is  thin  spared  to  the  overall  flow  field 
dimensions.  It  is  convenient  to  think  jf  the  undisturbed  Jet  flow  and  the 
undisturbed  air  flow  as  two  equivalent  regions  separated  by  the  inter¬ 
action  layer.  These  undisturbed  flows  may  be  determined  independently 
(the  Jet  flow  by  evaluating  the  expansion  into  a  vacuum)  and  then  used  as 
boundary  conditions  for  a  calculation  of  the  interaction  region.  In  this 
approach  the  important  characteristic  dimensions  during  the  calculation 
are  those  of  the  interaction  region  Itself  rather  than  the  less  appropriate 
(and  often  very  different)  radial  dimension  of  the  Jet. 

In  the  present  calculation  the  interaction  region  is  analyzed  in 
various  approximations:  first,  Li-rictly  inviscid  flow  end  second,  flow 
with  viscous  effects  included  in  the  so-called  thin  layer  approximation. 

The  supersonic  portion  of  the  inviscid  flow  is  computed  "exactly"  using 
finite  difference  techniques.  At  low  altitudes  the  air  shock  io  attached 
to  the  nozzle  edge  (neglecting  perturbing  influences  of  the  vehicle)  and 
the  sir  flow  is  entirely  supersonic  (described  by  hyperbolic  equations). 

At  sufficiently  high  altitudes  the  air  shock  detaches,  creating  a  subsonic 
region  (described  by  elliptic  equations)  at  the  bow  of  the  plume.  In 
addition  to  exhibiting  the  usual  problems  associated  with  the  numerical 
calculation  of  elliptic  regions  this  portion  of  the  flow  is  difficult  to 
define  for  other  reasons.  First,  FigureX* lb  indicates  that,  at  high  altitudes, 
the  nose  region  will  be  merged  (i.e.,  at  no  point  in  the  flew  is  the  inviscid 
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approximation  useful).  Some  treatments  of  such  a  region  are  available  * 
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but  it  is  not  clear  that  sophisticated,  treatment  is  appropriate  at  present. 
Since  moat  of  the  exhaust  'low  is  directed  backwards  only  &  small  mass 
fractior  of  the  total  exhaust  flow  determines  the  structure  of  the  noBe 
region.  This  flow  may  therefore  be  strongly  influenced  by  the  boundary 
layer  within  the  nozzle  and  by  the  additional  minor  mass  flows  associated 
with  the  normal  operation  of  a  missile:  vernier  engine  exhausts,  gas 
generator  exhausts,  etc.  Thus,  a  detailed  analysis  of  this  region  that 
ignores  these  effects  may  not  bo  warranted. 

Because  of  the  small  mass  flows  involved,  it  is  expected  that  this 
region  will  not  have  a  major  influence  on  the  downstream  flow.  Thus, 
relatively  crude  approximate  techniques  provide  r*\  adequate  description 
of  the  subsonic  region.  The  procedures  used  for  this  purpose  are  outlined 
in  8ectionIV-2.  Calculations  of  the  inviscid  flow  at  altitudes  up  to  lBO  km 
indicate  that  the  inviscid  approximations  of  the  exhaust  plumes  of  typical 
upper  stage  vehicles  Lave  relatively  sharp  noseB  at  all  altitudes.  Although 
the  turning  angle  of  the  flow  at  the  nozzle  lip  increases  monotoni colly 
with  increasing  altitude  and  may  appreciably  exceed  90°  for  some  nozzleB 
at  very  high  altitudes,  the  Jet-air  boundary  inclination  decreases  rapidly 
with  increasing  distance  from  the  nozzle  for  the  first  several  tens  of 
nozzle  diameters.  For  the  50:1  engine  (shown  in  Figure  II -10)  at  l6o  km 
altitude,  less  than  one  per  cent  of  the  air  flow  intercepted  by  the  maxi¬ 
mum  cross  sectional  area  of  the  exhaust  plume  enters  the  air  shock  layer 
at  points  where  the  boundary  inclination  exceeds  45*.  Under  these  can- 
dr  tioiifl  we  expect  rather  small  detachment  distances  and  small  BubBonic 
regions  (see  Figure  II -7) • 

Hozzle  boundary  layers,  gas  generator  discharges,  and  vernier  engine 
mass  flows,  in  addition  to  the  general  merged  layer  effects  in  the  subsonic 
region,  have  a  tendency  to  make  the  nose  of  the  plume  more  blunt.  A 
rough  estimate  of  the  Importance  of  these  mass  flows  and  of  the  assumptions 
made  in  the  calculation  of  the  subsonic  region  may  be  obtained  as  follows: 
The  peripheral  mass  flows  typically  constitute  one  to  three  per  cent  of 
the  total  exhaust  gas  muss  flow.  In  Figure  11-7,  we  no.*,  indicated  the 
positions  along  the  Jet  shock  at  which  various  fractions  (e)  of  the  total 
exhaust  gas  flow  have  been  entrained  by  the  Jet  shock  layer.  We  expect 
that  the  location  and  structure  of  the  shock  layers  in  the  region  close  to 
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the  nose,  where  less  than  five  per  cent  of  the  total  flow  has  been 
entrained,  will  deviate  appreciably  from  those  of  the  actual  vehicle 
because  of  the  neglect  of  these  minor  sue  flows .  far  downstream,  the 
shock  layer  structure  will  be  less  affected.  In  Figures  11-12  and  11-13 
we  have  Indicated,  at  the  dividing  streamline,  the  annular  width  (Ae)  which 
contains  five  per  cent  of  the  exhaust  mass  flow.  In  the  downstream  region 
it  is  apparent  that  this  width  is  small  compared  to  the  total  width  of  the 
mixing  region  (as  would  be  expected  for  a  merged  Jet  shock  layer).  Here, 
we  expect  that  the  major  structure  of  the  flow  will  not  be  strongly 
affected  by  these  minor  mass  flow* . 

In  Flgurei  ll-ll,  12,  13,  we  also  indicate  the  annular  widths  (a  )  required 
to  contain  a  mass  flux  of  air  equal  to  that  which  passes  through  the  sub¬ 
sonic  and  transonic  region  (M^15).  In  the  present  numerical  calculations, 
diffusion  in  this  region  is  not  accounted  for.  This  will  cause  appreciable 
errors  down  to  a  point  where  the  air  mass  flux  in  the  mixing  region  con¬ 
siderably  exceeds  that  passing  through  this  nose  region.  Reference  to 
these  figures  indicates  that  the  effects  cf  neglecting  the  diffusion  in 
the  subsonic  and  transonic  regions  will  not  introduce  much  error  in  the 
^valuation  of  most  of  the  downstream  flow  (for  this  particular  nozzle  and 
at  thet,“  high  altitudes).  In  other  words,  even  though  the  collision 
frequency  and  the  diffusion  rate  in  the  subsonic  region  are  sufficiently 
high  to  result  in  a  completely  merged  nose  shock  layer,  the  amount  of  mass 
involved  is  small  compared  to  the  total  amount  of  air  mixed  with  exhaust 
gases  throughout  much  of  the  flow  field.  When  the  additional  peripheral 
mass  flows  are  included  or  when  smaller  area  ratio  nozzles  are  considered, 
the  plume  tends  to  become  more  blunt-nosed.  This  will  increase  the  size 
and  importance  of  the  diffusive  phenomena  occurring  in  the  subsonic  region. 

At  present,  we  do  not  have  any  quantitative  estimates  of  how  important  these 
effects  will  be  for  configurations  markedly  different  from  the  rather  specific 
one  considered  here. 
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II -2  Results  of  the  Present  Calculations  of  the  Neutral  Flow 

In  Figures  II-2  to  II -8  we  show  the  results  of  inviscid  calculations 
of  the  shock  locations  and  the  contours  of  constant  air  density  at  three 
altitudes  (iOQ,  135y  and  l80  km^.  The  two  lower  altitudes  had  been 
calculated  previously  and  correspond  to  the  flow  resulting  from  the 
25:1  nozzle  labeled  #L  in  Figure  II -10.  The  flow  at  180  km  altitude 
was  evaluated  for  a  different  nozzle  (labeled  "5G:l-Nonuniform"  in 
Figure  II -10).  However,  the  jet  mass  flew  angular  distributions  are 
similar  for  the  two  nozzles  and  little  difference  is  expected  between 
them.  In  Figure  II -9,  we  show  the  shock  locations  and  air  density 
contours  at  l8o  km  evaluated  in  the  viscous  layer  approximation. 

Reference  to  this  figure  shows  that  there  is  surprisingly  little  cnange 
in  the  shock  locations  resulting  from  the  effects  of  diffusion.  In 
Figures  II -11  to  11-15  we  show  radial  distributions  of  density,  pressure 
and  tempo rature  on  three  different  orthogonal  surfaces  at  l80  kin 
(labeled  1,  2,  and  3  in  Figure  II -2).  In  Figures  II -l6  and  11-17,  we  show 
temperature  and  density  histories  along  two  streamlines  (a  and  b  in 
Figure  II -1).  A  number  of  qualitative  conclusions  may  be  drawn  from  an 
examination  of  these  figures; 

(1)  Strong  transverse  gradients  in  both  pressure  and  density  exist 
across  the  Jet  and  air  shock  layers. 

(2)  In  the  inviscid  approximation,  the  density  in  the  jet  layer  is 
very  much  higher  than  in  the  air  layer  although  the  pressure 
is  lower.  There  is  a  slight  entropy  layer  in  the  air  flow 
along  the  boundary  but  it  is  not  marked. 

(3)  In  the  viscous  layer  approximation,  the  jet  shock  layer  appears 
merged  throughout  the  flow.  The  width  of  the  mixing  layer  is 
cam*  rable  to  the  width  of  the  air  shock  layer  only  near  the 
bow  of  the  plume. 

■  5: 

At  low  altitudes,  the  electron  density  distribution  may  be  assumed 

proportional  to  the  air  density. 

^Figures  H- 2  to  II -9  are  presented  in  an  addendum  to  this  report 

(GD/C-DBFi>5-02 j(Addendum) ) . 
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(4)  The  effects  of  diffusion  have  surprisingly  little  influence 
on  the  air  density  profile  although  the  exhaust  gas  profile 
is  markedly  changed.  This  observation  suggests  a  model  of 
the  process  which  imagines  the  jet  as  a  cold  solid  body 
whose  surface  is  evaporating. 

(5)  The  cooling  effect  of  the  exhaust  gases  in  the  mixing  layer 
largely  eliminates  the  high  temperature  entropy  layer. 

(6)  The  pressure  gradients  at  the  jet  and  air  shocks  are 
sufficiently  large  that  the  shock  profiles  are  not  properly 
matched  to  the  interior  flow.  A  more  accurate  treatment 
would  have  to  include  the  effects  of  shoes,  cv-vature  on  the 
shock  profiles. 

(7)  The  nose  of  the  plume  is  relatively  sharp  and  the  subsonic 
region  small. 

(8)  The  Mach  disc  formed  by  the  reflection  of  the  jet  shock  from 
the  axis  is  small  compared  to  the  overall  jet  dimensions. 

A  preliminary  attempt  to  evaluate  the  far  field  motion -at  high 
alt ltudes  has  been  carried  out  in  the  small  disturbance  approximation. 

In  this  approximation  only  radial,  motions  are  considered  important  and 
the  problem  is  reduced  to  obtaining  an  equivalent  one  -dimensional 
unsteady  solution  of  the  Navier-Stokes  equations.  1  .e  numerical  calcu¬ 
lation  procedure  used  is  similar  in  the  type  of  approximat icn s  made  to 

* 

that  used  for  the  two-dimensional  flow.  Results  have  been  obtained  for 
one  very  high  altitude  case  (290  km  for  a  typical  sustainer  mass  flow)  and 
these  are  shown  in  Figures  II -18  to  11-21. 


"¥r 

This  case  is  a  limiting  case  -*  ®)  of  the  two-dimensional  steady 
flow  and  could  be  calculated  with  that  numerical  program.  This  limit 
results  in  several  simplifications  in  the  analysis  and  a  separate 
program  was  written-. 
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The  present  results  'nave  been  calculated  for  specific  vehicles  at 
specific  altitudes.  For  these  highly  under  expanded  jets  the  results  may 
be  scaled  to  other  altitudes,  velocities  and  maos  flows  under  certain 
conditions.  The  high  altitude  sealing  laws  are  discussed,  in  part,  in 
References  3  and  ^  .  Three  different  cases  may  be  distinguished.  In  order 
of  decreasing  generality  they  are:  the  near  field  inviscid  flow  when  the 
free  stream  is  hypersonic,  the  overall  inviscid  flow,  and  the  flow  field 
when  viscous  effects  are  important.  All  of  the  scaling  laws  depend  on 
the  fact  that,  at  these  high  altitudes,  the  undisturbed  flow  1c  source-like 
(p  ~  f(cp)  r-^)  and  the  jet  shock  may  be  considered  very  strong  at  almost  all 
points  in  the  flow. 

a)  The  inviscid  near  field  when  1. 

In  this  case  the  bow  portion  of  the  flow  field  has  a  characteristic 
dimension  R  given  by 


R  = 

where  m  is  the  rocket  exhaust  mass  flow,  its  limiting  velocity  (V^  =  J2h.^ 
where  h^  is  the  stagnation  enthalpy  per  gram},  and  and  the  free  stream 
properties,  'rfhen  the  flow  structure  is  evaluated  in  terms  of  the  dimension¬ 
less  coordinate  (r/R,  x/r),  the  shape  of  the  flow  field  is  independent  of 


A  V.r/oX 


altitude,  velocity  and  mass  flow  and  depends  only  on  the  angular  distri¬ 
bution  of  the  vacuum  flow  (*f  (9  )^)  and  the  specific  heat  ratios  of  the  jet 
and  air  streams.  Here  the  jet  and  air  densities  and  enthalpies  may  be 


represented  in  the  universal  forms 


pa  =  p«  <*/*  V  V  f(9)) 


h»“Vf  VYj'  f (»>) 

V  2  v 

PJ  *Cvj)  Pa  V  V  f(<p)) 

hJ  =  VJ  hQfc  Ya>  Yj.  *«) 

b)  The  inviscid  flow  for  finite  values  of  1<L 

The  scaling  law  Just  discussed  applies  to  that  portion  in  the  T'lume 
where  the  air  shock  may  be  considered  strong.  Far  downstream  the  air 
shock  deviates  appreciably  froc  this  limit.  For  the  flow  field  at  180  km 
nowhere  in  the  flow  is  this  condition  adequately  satisfied.  In  this  case 
the  flow  structure  also  depends  explicitly  on  the  value  of  the  free  stream 
Mach  number  M^.  It  is  expected  that  this  second  scaling  law  will  be  useful 
at  most  ionosphere  c  altitudes  for  correlating  phenomena  such  as  size  and 
shape  that  do  not  depend  on  the  detailed  structure  of  the  shock  layers. 

c)  The  flow  with  viscous  effects 

When  phenomena  which  depend  on  the  detailed  shock  layer  structure  are 
to  be  considered,  the  effects  of  diffusion  may  be  important.  Here,  in 
addition  to  the  free  stream  Mach  number  M^,  the  flow  structure  may  be  scaled 
from  one  altitude  to  another  only  if  the  velocity  ratio  Vj/V^  and  the 
Knudsen  number  ^proportional  to  Jh.  Vj  paj/V^  ^  are  kept  constant. 
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Method  Used  for  Evaluating  the  Inviscid  Flow  Field  In  ihe  Present  Study 


If 


A  finite  difference  approximation  of  the  inviscid  equations  of  motion 
is  used  to  evaluate  the  supersonic  portion  of  the  flow  field.  In  this 
numerical  calculation  tee  equations  of  motion  are  expressed  in  a  curvilinear 
coordinate  syste;m  composed  of  the  streamlines  and  the  associated  orthogonal 
surfaces  (potential  surfaces  in  irrotational  flow).  Of  the  four  differential 
equations  of  motion  for  two-dimensional  steady  inviscid  flow,  three  may  be 
integrated  directly  to  give  the  mass,  total  enthalpy  and  entropy  conserva¬ 
tion  equations.  These  integrated  equations  plus  the  differential  form  of 
the  lateral  momentum  equation  (which  determines  the  streamline  curvature) 
are  used  to  evaluate  the  flow.  Data  are  evaluated  on  successive  orthogonal 
surfaces. 

The  streamline  coordinate  system  was  chosen  for  the  finite  difference 
formulation  for  several  reasons.  First,  in  the  limiting  case  of  very  large 
Mach  number,  where  characteristics  programs  have  difficulty  because  the 
characteristics  often  do  not  intersect  in  a  diverging  flow  unless  the  mesh 
spacing  is  very  small,  the  streamtube  procedure  reduces  to  the  condition 
*'  the  streamline  curvature  is  zero  and  the  solution  is  obtained  trivially. 

Second,  by  using  the  integrated  conservation  equations,  it  can  be  assured 
that  these  are  rigorously  satisfied  no  matter  what  mesh  spacing  is  used. 

Third,  because  we  are  following  fluid  elements,  chemical  reactions  or  internal 
relaxations  may  be  easily  incorporated.  Also  slip  lines  are  treated  without 
difficulty.  Fourth,  because  of  ohe  marching  procedure  used,  lateral  diffu¬ 
sion,  heat  and  momentum  transfer  can  be  included  without  excessive  difficulty 
in  the  thin  layer  approximation. 

* 
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In  the  formulation  cf  the  finite  di^e-ence  procedure,  an  effort  has 
been  made  to  retain  a  physically  meaningful  Interpretation  even  when  the 
mesh  sizes  are  sufficiently  large  that  the  relation  between  derivatives 
and  finite  differences  is  somewhat  obscure.  For  this  purpose  a  streamline 
or  Iagrangian  coordinate  systems  offers  several  advantages  over  an  Eulerian 
coordinate  system  for  the  evaluation  of  hypersonic  flow  fields  even  though 
it  is  often  more  complex  and  less  efficient  in  terms  of  computing  time  pe1" 
mesh  point.  For  example,  in  the  computation  of  the  inviocid  flow,  the 
integrated  forms  of  the  contin’iity,  parallel  momentum  and  energy  equations 
are  used  so  that  these  are  rigorously  satisfied  no  matter  how  large  a  step 
size  is  used.  Truncation  errors  are  limited  to  the  lateral  momentum 
equation: 


3cp  _  1  dp 

“5s  ~  2  §y 

pu  J 


(1) 


This  equation  for  the  streamline  curvature  is  treated  as  follows :  Each 
streamtube  is  characterized  by  pressure  p  and  a  gas  velocity  |u|.  These 


are  related  to  the  streamtube  area  by  the  integral  energy,  momentum  and 

mass  relations.  The  curvature  C  of  the  streamline  separating  the 

v  ds  y12 

adjacent  streamtubes  1  and  2  is  evaluated  from  the  relation 


bn 

1  m  l'  V  ^*0  Pr>  ”  ^*1  Pi  ) 

ml  U1  +  m2  ^  ^  ^  c  1  1  y 


(2) 


where  rp  and  r1  are  the  radii  of  the  midpoints  of  the  two  streamtubes 
measured  along  a  surface  orthogonal  to  the  streamlines,  and  and  m^  the 
(constant)  mass  flows  in  these  tubes.  Thus  reasonably  accurate  evaluations 


of  the  flow  can  be  obtained  even  when  the  values  of  the  various  flow  quantities 


vary  considerably  from  tube  to  tube.  This  is  particularly  useful  in  the 
high  altitude  Jet  shock  layer,  where  t'r <*  total  pressure  may  vary  by  orders  of 
magnitude  and  the  density  by  a  factor  of  5  to  10  across  the  shock  layer 
while  the  curvature  and  flow  velocity  are  virtually  constant.  For  this 
case  Eq.  (2) would  give  accurate  results  even  if  the  entire  shock  layer 
were  represented  as  one  tube.  Since  the  computation  time  is  inversely 
proportional  to  the  square  of  the  tube  width  and  appreciable  computing  times 
are  involved,  even  on  the  7090,  it  is  important  to  make  the  mesh  size  as 
large  as  possible. 

In  the  present  calculations  relatively  simple  marching  procedures 
have  been  used  (at  present  tie  method  is  essentially  second  order  Runge- 
Kutta).  Thus  truncation  errors  can  arise  when  flow  variations  parallel 
to  the  streamline  are  large.  Although  the  procedures  can  be  readily 
improved  to  be  considerably  more  accurate  (i.e.,  fourth  order  Runge-Kutta) 
the  longitudinal  flow  variations  generally  encountered  are  usually  quite 
small  compared  to  the  marching  rate  permitted  by  the  stability  of  tl«e 
difference  scheme  so  that  the  present,  calculations  appear  to  be  stability 
limited  rather  than  accuracy  limited. 

An  attempt  has  been  made  to  incorporate  the  lateral  diffusion  in'*  the 
calculation  in  a  manner  similar  to  the  centrifugal  momentum  transport.  Here 
more  difficulties  are  encountered.  The  diffusion  effects  arc  incorpo¬ 
rated  by  evaluating  diffusion  fluxes  of  heat,  momentum  and  specie  at  the 
surface  of  each  stream '.ube.  Truncation  errors  enter  into  relation  between 
the  fluxes  and  the  differences  of  enthalpy,  velocity  and  concentration  in 
adjacent  streamtubes.  However,  the  calculation  is  arranged  so  that  gain  of 
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heat,  momentum  or  specie  by  one  tube  are  exactly  balanced  by  losses  in 
adjacent  tubes.  Thus,  so  long  as  the  longitudinal  variations  of  the  flov 
quantities  are  small  compared  to  the  stepping  rate,  most  of  the  truncation 
errors  can  be  thought  as  being  equivalent  to  using  slightly  erroneous 
effective  difiusion  coefficients. 

Since  the  missile  is  acce’  rating  and  the  atmosphere  is  not  uniform 
the  actual  plume  expansion  is  not  a  steady  phenomenon.  However,  the 
exhaust-a.4r  velocity  difference  far  downstream,  where  this  time  dependence 
is  expected  to  have  a  significant  effect,  is  smaller  than  the  missile 
velocity.  Thus  an  approximate  account  of  this  effect  for  a  vertical  tra¬ 
jectory  may  be  obtained  by  assuming  that  the  flow  disturbance  at  a  given 
altitude  is  the  same  as  that  for  a  rocket  moving  through  a  uniform  atmos¬ 
phere  at  constant  velocity.  Here  the  values  of  the  density  and  missile 
.elocity  must  he  taken  equal  to  those  corresponding  to  the  altitude  of 
interest.  (This  is  really  equivalent  to  a  small  disturbance  ipproximation 
far  from  the  vehicle.) 

Because  of  atmospheric  stratification,  the  disturbance  created  by  a 
vehicle  traveling  on  a  non-vertical  trajectory  will  not  be  axially  symmetric. 
An  approximate  account  of  this  effect  may  be  obtained  if  it  can  be  assumed 
that  the  azimuthal  flow  velocities  are  negligible.  In  this  approximation 
the  flow  distribution  at  a  given  azimuthal  angle  will  be  the  same  as  that 
in  an  "equivalent  axisymmetric  flow"  in  which  the  atmospheric  density  varia¬ 
tion  with  radial  displacement  is  assumed  to  be  the  same  as  tt  1 n  the  real 
atmosphere  in  the  given  azimuthal  direction.  In  other  words,  for  a  hori¬ 
zontal  trajectory,  the  stream  density  and  Mach  number  are  assumed  to  vary 
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with  radi-i?,  displacement.  For  an  exponential  atmosphere  the  "eqni valent 
axisymmetric  free  stream"  should  have  a  density  which  varies  with  radial 
displacement  from  the  nxi::  according  to 


/  v  -rcos  o'u 

p  v  r  )  ~  p  e 
“  ¥o 

where  6  is  the  azimuthal  argle  of  interest,  71  the  local  scale  height  and 


Pq  the  local  ambient  atmospheric  density  on  the  plume  axis.  Since  the 
finite  difference  formalism  is  the  same  for  uniform  as  for  nonuniform 


external  flows  this  method  is  easily  implemented. 
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II -h ,  Method  Used  for  Evaluating  the  Viscous  Flow  Field 

s 

Probetein  and  Kemp  classify  different  types  of  viscous  flows  into 
seven  regimes:  boundary  layer,  vorticity  interaction,  viscous  layer, 
incipient  merged  layer,  fully  merged  layer,  first  order  collision  theory 
and  free  molecular  flow  regimes.  In  tin*  first  three  regimes  the  shock  is 
treated  as  a  discontinuity  obeying  the  classical  Rankine-Kugoniot  relations. 

In  the  fourth  (incipient  merged  layer)  the  shock  structure  has  to  be  included. 
Here  suggested  treatments  are  either  tc  use  the  complete  N&vier-Stokes 
equations  to  describe  both  the  shock  profile  and  the  viscous  interior  as 
one  continuous  pattern,  or  tc  treat  the  shock  as  a  discontinuity  and  modify 
the  Hugoniot  relations  so  that  the  tangential  shear  stress  and  the  enthalpy 
flux  conducted  into  the  shock  are  properly  accounted  for.  Cheng  has  used 

Z 

this  latter  approach  successfully  to  treat  the  hypersonic  blunt  body  problem. 

In  treating  rarefied  gas  flows  a  very  useful  approximation  that  is 
often  introduced  is  the  so-called  thin  layer  assumption.  Basically  the 
assumption  says  that,  in  the  coordinate  system  chosen  to  evaluate  the  solu¬ 
tion,  the  viscous  stresses  and  heat  conduction  are  determined  by  tne  velocity 
and  enthalpy  gradients  in  a  direction  known  a  priori,  such  as  parallel  to  one 
of  the  coordinate  axes.  In  many  applications,  such  as  the  blunt  body  shock 
layer  or  the  rypersonic  slender  body  flow  in  the  small  disturbance  approxi¬ 
mation,  this  may  be  satisfied  if  the  shock  radius  of  curvature  is  much  greater 
than  the  lateral  dimensions  of  the  flow.  It  is  expected  that  this  condition 
applies  to  the  interaction  region  of  the  plume  as  long  as  the  free  stream 
Mach  number  is  sufficiently  high. 

We  have  incorporated  viscous  effects  into  the  numerical  calculation 
using  the  thin  layer  approximation .  Using  this  approach  in  the  incipient 
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merged  layer  regime,  Cheng  has  developed  an  elegant  treatment  of  the  hyper¬ 
sonic  blunt -body  problem  at  j.cw  densities. The  same  sophistication  of 

solution  is  presently  possible  for  the  high  altitude  Jet.  A  complete  evalua- 

* 

tion  of  the  jet  flow  field  would  proceed  as  follows;  first,  the  flow  is 
evaluated  in  the  tnviscid  approximation,  where  the  supersonic  region  is 
treated  exactly  and  the  subsonic  region  approximately.  Second,  viscous  and 
heat  transport  terms  are  included  in  the  viscous  layer  approximation  for  the 
supersonic  flow  (i.e.,  the  shocks  are  still  treated  as  discontinuous  a..d 
obey  Lhe  classical  Kugoniot  relations).  The  third  step  is  to  treat  the 
supersonic  flow  assuming  it  is  merged.  Here  Cheng's  two-layer  approach  is 
appropriate.  This  amounts  to  replacing  the  Hugoniot  relations  by  modified 
forms  which  allow  the  stagnation  enthalpy  and  the  tangential  velocity  to 
vary  through  the  shock.  At  the  same  time  the  shock  profile  itself  (shock 
transition  zone)  may  be  treated  independently  in  the  thin  layer  approxima¬ 
tion.  Here  the  shock  is  assumed  to  be  locally  approximated  by  a  plane 
shock  whose  profile  may  be  evaluated  by  a  direct  integration  of  the  one- 
dimensici?!  Navier-Stokss  equations,  using  the  local  conditions  evaluated 
at  the  edge  of  tne  merged  layer  as  init ial  conditions.  This  profile  is  then 
"tacked"  on  to  the  merged  layer  calculations.  This  is  essentially  the 
approximation  in  which  Cheng  treats  uhe  blunt-body  problem.  The  equations 
are  detailed  in  Section  IV-3a>b,c. 

Up  to  this  point  an  approximate  inviscid  description  of  tue  subsonic 
region  such  as  that  detailed  in  Section  IV-2c  is  used. 

The  logical  fourth  step  ir  this  procedure  would  be  to  use  Cheng's  method 
to  evaluate  the  subsonic  and  transonic  region.  Tne  only  difference  from  Cheng's 

During  the  present  study,  calculations  have  been  carried  out  only  in  the  first 
two  approximations. 
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procedure  would  be  that  instead  of  a  solid  boundary  for  the  innermost  stream¬ 
line,  a  boundary  condition  would  have  to  be  used  which  would  Join  the  elliptic 
air  region  to  the  hyperbolic  Jet  region.  A  merged  layer  calculation  would  be 
used  for  the  flow  inside  this  dividing  sti'eamline .  Diffusion  of  heat, 
momentum  and  species  across  the  dividing  streamline  would  be  allowed  for. 

The  equations  and  method  for  this  step  are  detailed  in  Section  IV-jd. 

During  the  present  study  it  v&r,  found  that  the  centrifugal  pressure 
gradients  across  both  the  Jet  and  the  air  shock  layer  were  large  (see 
Fig.  II -1^  throughout  most  of  the  flow.  Except  for  a  small  region  near  the 
nose  of  the  pi  me  the  pressure  gradient  in  the  air  shock  layer  at  180  km 
altitude  is  not  very  much  smaller  in  absolute  value  than  the  peak  gradient 
within  the  shock  transition  zone  itse1^  (Figs.  11-11,12  and  13).  Conse¬ 
quently  it  is  expected  that  the  curvature  of  the  shock  front  will  have  an 
appreciable  influence  on  the  shock  profile  and  it  will  be  necessary  to 
include  this  effect  in  order  to  get  a  consistent  fully  merged  layer  calcu¬ 
lation.  A  self  consistent  procedure  has  not  yet  been  formulated  for  this 
purpose.  Because  the  shock  layer  tangential  shear  stress  and  the  stagnation 
enthalpy  flux  are  relatively  small  near  the  shock  Interface,  except  near  the 
nose  of  the  plume  wlierc  the  deflection  angles  are  large,  it  is  expected  that 
the  influence  of  the  shock  curvature  (really  of  the  curveture  of  the  stream¬ 
lines  behind  the  shock)  on  the  shock  profile  will  be  the  dominant  source  of 
differences  between  a  viscous  layer  calculation  and  a  truly  merged  layer 
calculation.  Expressed  in  other  words,  the  difficulty  is  that,  in  the  air 
shock  transition  zone,  the  streamline  curvature  and  pressure  gradient  are 
of  opposite  sign  to  those  in  the  shock  layer  and  a  plane  shock  approximation 
results  in  a  discontinuity  in  these  quantities.  In  practice  this  means  that 
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the  peak  pressure,  density  and  streamline  deflection  calculated  for  a  plane 
shock  will  not  be  achieved  in  the  actual  flow.  The  detailed  calculations 
to  date  have  been  limited  to  the  viscous  layer  approximation.  The  density 
contours  shown  in  Figures  Il-j  to  8  in  the  neighborhood  of  the  shock  inter¬ 
face  were  determined  by  fairing  the  shock  layer  profile  smoothly  into  a 
shock  transition  zone  profile  calculated  assuming  zero  shock  curva+"re. 

For  the  jet  shock  this  difficulty  is  less  apparent  since  the  curvatures 
in  the  transition  zone  and  in  the  shock  layer  have  the  same  sign.  In  this 
case,  it  would  be  consistent  to  locate  the  shock  interface  at  the  point  in 
the  shock  transition  profile  where  the  streamline  curvature  matched  the 
value  in  the  shock  layer.  However,  the  viscous  layer  calculations  carried 
out  during  this  study  have  not  included  this  sophistication  and  it  has  been 
assumed  that  the  normal  pressure  and  density  jumps  were  achieved  across  the 
shock.  The  errors  induced  are  expected  to  be  relatively  small,  perhaps 
equivalent  to  a  shift  in  the  position  of  the  interaction  layer  by  an  amount 
somewhft  less  than  the  width  of  the  shock  transition  zone. 
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Figures  II-2  to  II -r'  appear  in  an 
addendum  to  this  report  -  GD  C-DBE65-023 
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Fig.  II-20.  Air  pa:  iai  den  ity  distribution  a’  various  tUes  (in 

seconds)  la  the  small  disturbance  approximation  at  290  km. 
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FIG.  E-21,  FLOWFIELD  AT  290km  IN  THE  SMALL  DISTURBANCE  APPROXIMATION 


VELOCITY  (km /sec) 


III.  TUe  Effect  of  the  Geomagnetic  Field 


ITI-1.  General  Description 

The  ultimate  purpose  of  the  present  study  is  to  determine  the  electron 
density  distribution  in  the  vicinity  of  a  powered-rocket  in  the  ionosphere. 
Since,  throughout  much  of  the  altitude  range  of  interest^the  Larmor  radius 
of  both  ions  and  electrons  may  be  considerably  smaller  than  a  mean  free 
path,  the  geomagnetic  field  may  considerably  influence  their  motion. 

Lighthill^  I  ^  has  considered  the  disturbance  that  would  be  induced  by  the 
passage  of  a  satellite  through  the  ionosphere  and  concluded  that  the 
perturbation  produced  would  be  primarily  a  one- dimensional  motion  of  the 
ions  parallel  to  the  magnetic  field.  Although  we  are  concerned  with  a 
somewhat  wider  range  of  altitudes,  the  neutral  wind  created  by  the  passage 
of  a  rocket  provides  a  similar  but  much  stronger  perturbation  to  the  ion 
fluid. 

It  is  known  that  the  motion  of  a  plasma  in  a  magnetic  field  can  be 
quite  sensitive  to  the  boundary  conditions  imposed.  The  effects  are  primarily 
due  to  the  long  range  of  the  electromagnetic  body  forces .  In  our  initial 
studies  we  have  limited  ourselves  to  two  simplified  approaches .  The  first 
is  to  consider  sufficiently  simple  geometries  for  both  the  neutral  flow  field 
and  the  electromagnetic  field  that  it  is  possible  to  obtain  "exact"  numerical 
solutions  to  the  non-linear  equations  of  motion.  The  second  is  to  consider 
approximations  of  aquations  of  motion  that  will  allow  consideration  of  more 
general  symmetries  and  allow  at  least  a  semi -quantitative  evaluation  of  the 
importance  of  properties  of  the  neutral  flow  and  the  plasma  (such  as  non- 
uniform  distributions)  that  it  is  difficult  to  include  in  the  detailed 
solutiors . 
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For  a  detailed  solution,  we  have  chosen  to  c;  .sider  firs  +  the  case 
of  a  cylindrically  symmetric  neutral  motion  in  an  axial  magnetic  field. 

The  basic  a  priori  assumption  that  is  made  is  that  the  revel  of  ionization 
is  sufficiently  lew  that  the  motion  of  the  neutral  species  is  completely 
unaffected  by  the  presence  of  uhe  ionized  species.  In  this  ca”e  the 
neutral  motion  may  be  calculated  without  regard  to  the  presence  of  the  ionic 
species.  The  ion  motion  may  then  be  determined  by  considering  the  ion  fluid 
as  a  fully  ionized  plasma  which  moves  under  the  influence  of  a  strong  body 
force  (the  net  momentum  transfer  from  the  prespecified  neutral  wind  to 
the  ionized  fluid)  in  addition  to  the  electromagnetic  field  and  its  own 
pressure  gradients.  In  addition,  appreciable  heat  transfer  to  the  ion  fluid 
may  take  place  as  a  result  of  the  viscous  stresses  generated  between  the 
ions  and  the  neutrals . 

The  condition  for  the  neutral  fluid  to  be  independent  of  the  ion 

motion  can  be  estimated  as  follows .  Suppose  that  the  magnetic  field  is 

* 

sufficiently  strong  that  «  V.  Then  the  force  exerted  on  the  ions  by 
the  neutrals  per  unit  mass  of  the  neutrals  is  given  by 

f 

Thus  the  neutral  velocity  will  be  appreciably  reduced  by  this  drag  in  a 
time  of  the  order 

t  w  P/P 


Thus 


,  for  frequencies  less  than  -p-  v^,  the  neutral  and  ionized  species  move 


together  as  a  single  fluid.  This  frequency  does  not  vary  very  strongly  with 
altitude  and  at  180  km  altitude  is  of  the  order  of  10  J  sec  \  Tims,  for  the 
much  shorter  time  of  present  interest,  we  may  ignore  the  influence  of  the  ions 


Here  u,  and  V  are  the  radial  components  of  the  ion  and  neutral  velocity, 
respectively.  The  ion  mass  density  is  Pi,  the  neutral  density  is  p  arid  the 
ion-neutral  collision  frequency  is  v^. 


Ill -2 


on  the  neutrals . 


If  there  were  no  magnetic  fieid  present,  the  ions  would  be  pretty  much 
swept  along  with  the  neutrals.  When  the  magnetic  field  is  considered,  we 
must  allow  for  the  possibility  that  the  ion  motion  across  the  field  lines 
may  be  strongly  hindered,  A  rough  picture  of  the  phenomena  involved  may 
be  obtained  by  considering  the  expansion  of  a  cylindrically  symmetric 
neutral  volume  of  gas  into  the  ionosphere  when  the  magnetic  field  is  aligned 
parallel  to  the  axis  of  symmetry.  If  the  outward  velocity  of  the  neutrals 

is  V  and  the  ion  velocity  is  u^,  then  the  force  per  unit  volume  exerted  by 

the  neutrals  on  the  ions  due  to  the  ion  slip  will  be  of  uhe  order 

fc  (V  -  vr) 

where  v  is  a  quantity  of  the  order  of  the  ion-neutral  collision  frequency 
(and  may  also  depend  on  V).  This  force  will  cause  an  azimuthal  drift  of 
the  ions  and  electrons.  For  a  weakly  ionized  gas,  this  force,  even  when 
exerted  over  a  depth  of  several  kilometers,  is  small  compared  to  the  magnetic 
pressure  B  /8n ,  For  example,  if  n^  -  10  cm  ,  vin  =  10  see  ,  V-u^  =  2  x 

10^  cm/sec  over  a  5  km  distance,  then  the  force  per  unit  area  is  of  the  order 

of  3  *  10  dynes/cm  for  an  ion  mass  of  16.  Thus  a  slight  gradient  in  the 
magnetic  field  (.0002  gauss/km)  is  sufficient  to  balance  this  viscous  stress. 

Let  us  first  consider  the  case  when  the  ion  mass  and  electron  mass  are 
equal,  as  are  their  interactions  with  the  neutral  fluid.  Here,  because  of  the 
symmetry,  we  may  conclude  that  no  charge  separation  will  tend  to  occur.  Thus 
both  the  axial  and  the  radial  electric  field  vanish  identically.  The  azimuthal 
drift  velocity  will  be  of  the  order 

-  CWi  (v~u?) 

~  ft  3 
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/Wain,  due  t,o  collisions  with  neutrals,  there  wilJ  be  an  azimuthal  drag 
l'orce  of  the  order 


vin 


Thus,  the  outward  drift  velocity  will  be 


«*,  .  (  vin/wt. )’"(  V-u,) 


or 


V, 


v£ 


V 


0) 


where  u'  -  eB/m ,  c  . 
c  i 

1 

rirus  strong  hindrance  of  the  motion  across  the  field  lines  is  expected 

when  v.  <  uj  (~  249  sec  1  for  B  =  0.4  gauss  ana  an  ion  mass  of  l6). 

in  c. 

When  the  ion  and  electron  masses  are  unequal ,  radial  fields  due  to 

charge  separation  may  occur.  These  may  considerably  influence  the  motion. 

A  very  slignt  difference  between  the  ion  and  the  electron  density  can  create 

large  electric  fields.  Thus  we  expect  that  n  ng  to  a  goci  approximation. 

This  condition  implies  that  the  ions  and  electrons  will  move  radially  with 

* 

the  same  velocities.  This  strong  coupling  of  the  two  motions  results  in 
nighly  distorted  Larmor  orbits.  If  the  ions  were  acted  on  by  a  velocity- 
independent  radial  force  and  the  collision  frequency  were  small,  the  two 
orbits  would  be  azimuthally  drifti’ig  ellipses  having  a  ratio  of  semi-minor 

f% 

to  semi -major  axis  of  J  —  .  The  electron  orbit  ..ould  be  elongated  in  the 

f  m . 
i 

azimuthal  direction  whereas  the  on  orbit  is  elongated  in  U  radial 


This  conclusion  depends  specifically  on  the  assumption  that  neither  the  neutral 
motion  V  nor  the  plasma  properties  vaiy  in  the  direction  of  the  mejoetic  field. 
As  is  shown  in  Section  V-3  this  requirement  c_n  be  approximated  in  practice  only 
if  toe  f Jew  field  dimensions  are  extraordinarily  large, 
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In  this  case  it  can  be  shorn  that  the  outward  drift  velocity  has 


the  form 


•**  ^ 


V 


(2) 


where  nr.  is  the  ,  ybrid  Larmor  frequency  (id..  =  uj  uu  )  and  v  is  a 
K  K  c.  c  en 

i  e 

quantity  cf  the  order  of  the  electron-neutral  collision  frequency.  Thus 
vp.en  the  negative  ion  (electron)  mars  is  small,  the  coupling  to  the  magnetic 


field  is  stronger.  Taking  v 


en 


.  F 

v  m 


v  .  .we  find  that  appreciable  ion  slip 
m  m 

i/4  _  -1 


will  occur  when  <  (m^/r^)  ■  u)c  (=“  3100  sec  for  =  16).  In  Table  III-l 
we  give  the  ion-neutral  collision  frequency  in  the  ambient  io-  sphere  at 


various  altitudes  assuming  v 


m 


-1^ 

6  X  10  n  where  n  is  the  neutral  number 
a  a 


density.  Also  in  Table  III-l  we  indicate  the  collision  frequency  that  might 
be  expected  behind  a  moderate  strength  shock  (taken  to  be  &  10  times  that  in 
the  ambient).  It  will  be  shown  later  that,  because  of  the  finite  extent  of 
+he  neutrl  flow  in  the  direction  of  the  magnetic  field,  the  appropriate 

Reference  to  Table  III-l  then  shows  that 
appreciable  slip  effects  will  occur  at  altitudes  above  120  km  in  the  less 
dense  parts  of  the  flow  field  and  at  altitudes  above  150  km  even  the 
denser  part <5  cf  the  flow  f-'eld. 


coupling  parameter  is  \r  n  /oj^ 


Table  III-l 


Alt  j.tude 

km 

na(cn:"')) 

V  . 

m 

ambient 

(sec  ■*") 

shock 

Ui 

,  \ 
(sec  ) 

m4 

(— )  » 

vm  '  Cj 

e  i 

100 

8  x  IO1" 

4800 

48000 

240 

5100 

120 

3  x  1011 

180 

1800 

150 

4  x  KT0 

24 

240 

200 

8  x  109 

4.8 

48 

250 

3  x  109 

1.8 

18 

300 

Q 

1  X  10- 

0.6 

6 

I 
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The  discussion  given  In  the  preceding  section  indicates  that,  at  high 
altitudes,  l  geomagnetic  field  strongly  hinders  the  transverse  motion  of 
the  ions  in  the  less  dense  portions  of  the  flow  field.  The  neutral,  flow, 
which  is  essentially  independent  cf  the  presence  of  the  magnetic  field, 
acts  as  a  driving  wind  which  stresses  and  heats  the  ion  fluid.  When  the 
suppression  of  the  transverse  motion  of  the  ions  is  strong,  the  ions  move 
one -dimensionally  along  the  magnetic  field  lines  under  the  influence  of 
this  three -dimensional  wind.  Except  in  the  region  cf  the  neutral  flow 
where  there  are  strong  density  gradients  (shock  fronts),  the  characteristic 
times  for  the  motion  are  lonD  compared  to  the  time  between  ion-neutral 
collisions.  In  this  approximation  it  is  expected  that  the  ion  and  electron 
temperature  and  pressure  distributions  will  be  reasonably  isotropic  and 
Maxwellian  (see  Section  III -3  although  the  individual,  values  for  the  ion 
and  electron  temperatures  may  differ  considerably  hath  from  i'ne  neutral 
gas  temperature  and  from  each  other. 

Also  when  the  motion  of  the,  ions  across  the  magnetic  field  is  strongly 
suppressed,  the  lateral  transport  of  heat  and  momentum  is  markedly  reduced. 
In  this  situation  the  problem  of  evaluating  the  three-dimensional  electron 
density  distributions  may  be  reduced  to  a  number  of  independent  one- 
dimensional  flow  problems  once  the  underlying  neutral,  flow  field  has  been 
evaluated.  Here  the  magnetic  field  acts  simply  to  confine  the  ionized 
fluid  to  essentially  rigid,  frictionl ess  channels.  The  ionized  gas  then 
executes  a  one -dimensional  motion  under  the  influence  cf  the  momentum  and 
heat  input  from  the  predetermined  neutral  wind.  In  this  motion  the  electro¬ 
magnetic  fields  play  r.o  role  (except  for  ambipolar  diffusion  effects). 

At  high  altitudes  we  expect  that  this  model  of  the  motion  will  apply  to 


\ 
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much  oi  the  flow  field.  Close  to  the  vehicle,  however,  will  always  be 
a  region  where  the  gas  density  is  suff iciently  high  to  ensure  that  the 
ioi.  j  from  the  ambient  ionosphere  are  swept  along  with  the  neutral  exhaust 
gases  regardless  of  the  magnetic  field.  For  typical  vehicles  it  is 
expected  that  this  region  may  be  of  the  order  01  a  kilometer  in  diameter. 

This  model  of  the  perturbations  induced  in  the  ionosphere  by  the 
passage  of  the  vehicle  was  first  proposed  by  Lighthill^  J"  ^  to  describe 
the  perturbation  to  the  F-region  electron  density  distribution  created 
by  the  passage  of  a  satellite.  The  only  difference  lies  in  the  magnitude 
of  the  effect  and  the  nature  of  the  source.  Lighthill  was  primarily 
concerned  with  the  far  field  disturbance  induced  by  a  relatively  small 
body  travelling  in  the  u-  er  regions  of  the  F  layer.  In  our  case,  the 
source  of  the  disturbance,  the  neutral  wind,  is  distributed  over  a  large 
volume  and  is  expanding  at  a  speed  somewhat  greater  than  the  local  sound 
speed.  Thus  thq  acoustic  type  of  motion  induced  In  the  ionized  fluid  will 
occupy  roughly  the  same  volume  as  the  neutral  flow  field  although  the 
velocities  and  density  distribution  of  the  ions  may  be  .  arkedly  dif * .rent 
f x'Oin  tho^e  of  the  neutials. 

Since  the  ion  motion  is  essentially  one-dimensional,  there  is  little 
geometrical  "attenuation"  of  the  disturbances  induced  in  the  ion  fluid  as 
compared  to  the  l/r  or  l/r~  stter.uation  of  the  neutral  flow.  Thus,  if  there 
were  negligible  viscous  dissipation  of  the  motion,  the  disturbance  generated 
in  the  ion  fluid  could  travel  long  distances  while  retaining  significant 
amplitude.  However,  this  effect  is  not  expected  jn  the  present  case  since 
the  ion  velocity  is  strongly  coupled  to  the  component  of  the  neutral  flow 
velocity  parallel  to  the  magnetic  field,  A  disturbance  in  the  ion  fluid 
which  propagates  into  a  region  where  the  neutral  flow  velocity  is  small 
would  be  damped  within  a  few  ion-neutral  mean  free  paths. 

IIT-7 


Some  non-acoustic  types  of  disturbances  may  propagate  beyond  the 
confines  of  the  neutral  disturbance.  The  high  thermal  conductivity  of  the 
electron  gas  will  cause  the  heat  generated  by  the  friction  between  the 
neutral  wind  and  the  electron  gas  to  be  spread  throughout  a  volume  which 
may  be  much  larger  than  that  of  the  neutral  flow  field.  The  effect  is 
similar  to  the  broad  thermal  layer  in  the  elect. ron  gas  which  occurs  in  front 
of  a  shock  travelling  through  a  fully  or  partially  ionized  gas.  (.Jaffrin  (i!) 
Probstein  and  Jaffrin  (j)  ).  In  the  piesent  case  it  is  expected  that 
these  conduction  lengths  will  be  considerably  larger  than  the  dimension 
of  the  neutral  flow  field  for  most  altitudes  of  interest.  This,  together 
with  the  fact  that  the  mechanism  of  thermal  energy  transfer  from  the 
neutrals  to  the  electrons  is  rather  inefficient  because  of  the  small  electron- 
neutral  mass  ratio,  implies  that  the  electron  smperature  will  not  be  signi¬ 
ficantly  altered  by  the  friction  of  the  neutral  wind.  Because  of  this 
high  conductivity,  compression  or  expansion  of  the  ionized  fluid  is  also 
not  expected  to  alter  the  electron  temperature  appreciably.  Thus  it  appeal’s 
that  the  electron  gas  may  be  considered  isothermal  with  a  temperature  equal 
to  that  of  the  ambient  ionosphere. 

Adthough  the  high  thermal  conductivity  of  the  electron  gas  parallel  to 
the  magnetic  field  dees  not  appear  to  give  rise  to  any  effects  outside  the 
region  of  the  neutral  flow,  significant  effects  ma’  arise  as  a  result,  of 
the  correspondingly  high  parallel  electrical  conductivity.  Since  the 
accelerations  imparted  to  the  positive  ions  due  to  their  interaction  with 
the  neutral  wind  are  generally  different  from  those  given  to  the  electrons, 
charge  separation  will  tend  to  occur  within  the  confines  of  the  neutral 
flow  field.  Because  of  the  high  longitudinal  conductivity,  these  charges 
will  spread  rapidly  along  the  field  lines  until  they  are  eventually 
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neutralized  by  the  small  cross  field  leakage  current.  It  is  expected  that 
tnis  spreading  constitutes  an  Alfven  motion  guided  along  the  field  lines 
in  a  manner  related  to  the  magnetonydrodynamic  whistler  mode.  However, 
the  characteristic  frequencies  of  this  disturbance  are  considerably  smaller 
than  ohe  collision  frequency.  Thus  the  motion  will  be  more  diffusive  than 
wave -like. 

Since  the  velocities  associated  with  the  neutral  flow  are  very  much 
less  than  the  ALfv4n  speed  (the  Mach  number  based  on  the  Alfven  speed  is 
of  the  order  of  0.01),  the  amplitude  of  the  disturbance  will  be  small. 

The  nature  of  these  far  field  disturbances  is  discussed  further  in 
Section  V-6.  A  preliminary  estimate  of  the  disturbance  expected  in  the 
steady  state  for  an  axially  symmetric  neutral  flow  in  an  aligned  field 
indicates  that  the  perturbation  to  the  electron  density  beyond  the  confines 
of  the  neutral  flo-T  may  be  of  the  order 

+  x  v,y ( v,1  * -£ ) 

<¥V*  J 

where  f  is  the  reduction  factor  defined  in  Section  V-3*  This  probably  is 
too  small  to  be  noticeable.  In  a  transverse  magnetic  field,  larger  effects 
are  expected.  Also  the  time -dependence  of  the  Alfv4n  motion  and/or  inhomo¬ 
geneities  in  the  density  of  the  ambient  ionosphere  may  give  more  pronounced 
effects,  especially  when  the  collision  frequency  is  small.  Further  analysis 
will  be  required  to  evaluate  the  level  of  the  field  motion. 

In  this  connection  it  is  interesting  to  note  that  Drell,  Foley 

(  4  ) 

and  Ruderman  have  recently  discussed  a  similar  problem.'  1  They  were 
concerned  with  the  interaction  between  a  large  satellite  and  the  ionosphere 
and  invoke  an  Alfven  type  of  motion  radiated  by  the  satellite  to  account 
_ jr  its  observed  drag.  Except  for  the  fact  that  collisions  are  less 
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important  in  their  cane  because  of  the  much  higher  altitude  and  that 
the  source  (ECHO  l)  is  small  and  of  fixed  size,  the  phenomena  are  quite 
closely  related:  in  both  cases  the  source  acts  as  a  battery  which 
drives  currents  alcxig  the  conducting  field  lines.  These  authors  expect 
a  perturbation  in  the  electron  density  as  large  as  15$  to  be  propagated 
along  the  field  lines  up  to  several  hundred  kilometers  from  the  satellite. 
The  present  analysis  has  not  been  carried  sufficiently  far  to  give  any 
real  indication  whether  or  not  any  perturbation  approaching  this  magnitude 
may  be  induced  by  the  wind  created  by  a  large  rocket. 
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Relative  motion  of  ions  and  neutrals  at  various  inclination 
of  trajectory  to  magnet i  ■  field. 


Ill  -1-2. 
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IV  Detailed  Description  of  the  Neutral  Flow 
IV-1.  Background 

During  recent  years  an  extensive  literature  dealing  with  the 

1-3 

structure  of  high  altitude  rocket  plumes  has  developed.  Experimental 

measurements ,  "exact"  numerical  methods  and  approximate  analytical 

methods  have  been  used.  For  example,  schlieren  and  shadowgraph  techniques 

have  been  U'ed  to  study  underexpanded  jets  expanding  into  still  atmospheres 

by  Iatvala,^,'>  Wilcox  et  al.  /  D’Attorre  and  Harshoa’-ger  J  because  of 

experimental  difficulties  laboratory  data  are  generally  limited  to  exit 

to  ambient  pressure  ratios  less  1000  and  in  most  cases  less  than  IOC. 

These  are  much  less  than  the  values  for  a  missile  at  ionospheric  altituu  s. 

The  laboratory  data  for  jets  in  highly  supersonic  streams  is  meager. 

Rosenberg  has  used  photograph! c  techniques  to  define  tne  outline  of  the 

plume  at  high  altitude  during  a  number  of  missile  launches.  From  his  data 

he  has  determined  an  approximate  size  and  shape  of  the  plume  boundary  (as 

\  9 

defined  by  its  luminosity;.  Hill  and  Babert  have  obtained  approximate 
correlation  of  Rosenberg’s  data  using  blast  wave  theory. 

Ko  attempt  at  a  complete  review  of  the  existent  teclniques  will  he 

IO 

made  here  in  view  of  the  several  studies  presently  availaole.  Adamson 

has  reviewed  much  of  the  unclassified  literature  on  the  inviscid  aerodynamic 

II 

structure  of  the  rocket  exhaust  plume.  Rosner  gives  a  more  recent  detailed 

survey  of  the  state  of  the  art  for  predicting  plume  structure  at  various 

altitudes  with  and  without  viscous  effects.  In  general,  reliable  procedures 

-ire  available  for  certain  portions  of  the  plume  but  not  for  others.  The 

expansion  into  vacuum  when  the  nozzle  flow  is  reasonably  uniform  has  been 

treated  by  a  number  of  people  using  uhe  method  of  cV.vacterist. c0  for  the 
i 

i 
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flow  close  to  the  nozzle  and  an  asymptotic  analysis  for  the  source-like 
flow  at  great  distances.  In  the  presence  of  a  low  but  finite  pressure  free 
stream,  a  complex  sncck  structure  develops.  Upstream  of  the  point  where 
the  jet  shock  returns  to  the  axis  the  flow  lias  a  simple  structure  and  has 
been  treated  by  both  exact  and  approximate  methods.  The  reflection  of  the 
jet  shock  from  the  axis  of  symmetry  typical"^  results  in  the  formation  of 
a  shock  triple  point  and  a  Mach  disc.  The  flow  downstream  of  this  point 
contains  both  supersonic  and  subsonic  regions.  Little  detailed  analysis 
has  been  successfully  applied  to  this  region  for  highly  unuerexpanded  jets. 

The  analytic  techniques  may  be  divided  into  three  categories:  linea¬ 
rized  methods,  approximate  methods  and  "exact"  or  detailed  numerical  methods 
for  obtaining  solution  of  the  equations  without  approximation.  The 
linearized  theories  are  concerned  with  small  degrees  of  under  and  over¬ 
expansion  and  are  nut  of  primary  interest  to  the  present  problem.  The 
approximate  methods  attempt  to  obtain  information  about  the  plume  structure, 
such  as  boundary  and  shock  locations  or  the  far  downstream  average  properties 

without  requiring  a  detailed  calculation  of  the  entire  flow  field.  The 

Q 

boundary  location  procedures  are  reviewed  in  detail  by  Adamson.  In  general 
the  approximations  Involve  a  method  for  specifying  a  jet  pressure  and  a 
stream  pressure  at  the  boundary  and  demanding  equality  of  these  pressures. 
Newtonian  or  Prandtl  Meyer  expansion  approximations  have  been  used  to  specify 
these  pressures.  These  procedures  appear  to  be  capable  of  predicting, 

withir  a  reasonable  approximation,  the  boundary  location  not  too  far  from 
the  nozzle,  but  the  approximations  typically  worsen  downstream.  Hill  and 

9 

Habert  have  been  able  to  obtain  reasonable  correlation  of  the  high  altitude 
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plume  shape  using  blast  wave  theory.  Thomson  has  used  an  asympto+ic  analysis 

to  correlate  the  properties  of  the  far  downstream  pressure  equilibrated  jet 

«3 

to  the  nose  shape.  Alden,  Hill  and  Habert  treat  the  shock  layer  as  thin 

and  obtain  a  universal  plume  shape  at  high  altitudes  for  hypersonic  free 

14 

streams.  AxMni  uses  an  inviscid  two  shock  layer  formulation  for  the  very 
high  altitude  jet  which  includes  centrifugal  relief.  He  also  assumes  thin 
shock  layers  and  a  hypersonic  free  stream.  It  is  expected  that  these  approxi¬ 
mate  techniques  will  yield  useful  information  about  the  gross  plume  structure 
but  how  accurately  they  represent  the  detailed  profiles  is  questionable. 

The  basic  difficulty  with  maiiy  of  the  approximate  techniques  that  have 
been  developed  is  that  the  validity  of  the  various  approximations  is  difficult 
to  ascertain  a  priori.  In  many  cases  a  detailed  assessment  of  the  accuracy 
of  the  approximations  has  also  not  been  made  a  posteriori,  often  because  the 
results  sometimes  appear  to  agree  well  with  certain  features  of  observed 
plume  shapes  or  more  detailed  numerical  calculations  of  the  plume  shape.  The 
difficulty,  of  course,  is  that  the  agreement  is  often  limited  to  one  quantity, 
such  as  shape.  A  classic  example  is  that  of  the  blast  wave  similarity 
approximation  of  a  cylindrical  or  spherical  explosion.  Here  the  location 
of  the  shock  front  is  well  described  by  the  approximate  theory  long  after 
the  similarity  approximations  break  When  the  details  of  the  internal 

structure  are  compared,  however,  disparities  often  appear.  For  the  present 
problem  the  detailed  structure  is  oi  primary  interest.  For  this  purpose  it 
appears  that  an  "exact"  numerical  solution  is  required  in  order  to  obtain 
reliable  information. 

Several  numerical  programs,  using  the  method-of-characteristies,  have 
ISM7 

been  developed.  Gome  of  these  are  isentropic  programs  and  are  most  useful 
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for  nozzle  or  vacuum  expansion  calculations  or  for  nearly  balance. 1  nuzzles 
where  the  shocks  are  weak.  Others  treat  the  air  and  jet  i.  he 'k  as  non- 
isentrcpic  discontinuities  assuming  the  air  shock  is  attached  to  t.ue  nozzle 
lip.  When  the  free  stream  is  hypersonic  and  the  air  shock  detaches,  creating 
a  subsonic  region,  the  jet  flow  alone  may  be  evaluated  by  using  the  Newtonian 
approximation  to  specify  the  boundary  pressure.  When  the  Newtonian  approxi¬ 
mation  is  valid  the  jet  flow  may  be  computed  first, yielding  a  blunt  body 

ia 

contour  around  which  the  air  flow  may  be  evaluated  subsequently.  In 
the  present  problem  the  shock  may  be  detached  so  that  the  Newtonian  approxi¬ 
mation  is  useful  at  the  forward  part  of  the  plume,  but  the  Mach  numbers  are 
not  high  enough  that  this  is  a  good  approximation  for  the  downstream  portion 
of  the  plume.  At  present  there  does  not  appear  to  exist  a  complete  "exact" 
numerical  procedure  for  treating  this  problem. 

There  are  very  few  treatments  of  a  highly'  underexpanded  jet  including 

diffusion  effects.  Vasiliu  lias  formulated  the  problem  in  the  boundary  layer 

19 

approximation  (.very  thin  mixing  layer).  However,  the  procedure  is  an  implicit 
one  and  is  awkward  to  apply.  Some  calculations  have  been  carried  out  for  a 
turbulent  thin  boundary  layer  along  +he  plume  boundary.  In  any  case,  the 
strong  lateral  gradients  that  occur  at  high  altitude  make  the  boundary  layer 
approximation  suspect. 

The  approach  used  in  the  present  study  is  a  combination  of  procedures 
that  have  been  developed  for  axi symmetric  jet  flows  and  for  hypersonic  flows 
over  solid  bodies.  In  the  following  paragraphs  we  discuss  some  of  the 
characteristics  of  the  available  procedures  as  they  pertain  to  the  present 
S  problem. 

Within  the  limits  of  continuum  theory  the  Navier-Stokes  equations  are 
usually  assumed  to  provide  a  reasonable  description  of  moderately  rarefied 
.  flow  fields.  Questions  exist  concerning  the  validity  of  these  equations 
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for  the  structure  of  very  strong  shocks  tut  at  present  no  sui cable  alternative 
lms  been  presented.  Many  different  useful  methods  have  been  developed  for 
obtaining  exact  or  approximate  solutions  to  these  equations  for  various 
boundary  conditions.  One  powerful  tool  has  been  the  use  of  the  similarity 
assumption.  The  similarity  technique  assumes  that  a  separation  of  variables 
is  possible  in  the  differential  equation  and  the  boundary  conditions.  With 
this  assumption  the  partial  differential  equation  is  reduced  to  two  sets  of 
ordinary  differential  equations  which  may  be  integrated  by  standard  techniques. 
In  general  similarity  may  be  obtained  only  for  certain  types  of  boundary 
conditions.  The  similarity  technique  has  proved  a  powerful  tool  for 
treating  the  hypersonic  flow  over  slender  bodies.  Here,  when  the  small 
disturbance  approximations  are  valid,  (6  «  1  and  5  «  1  where  6  is  a 
characteristic  body  slope  and  tlie  free  stream  Mach  number)  solutions 

70 

to  an  equivalent  unsteady  one  dimensional  flow  may  be  sought.  Similar 
21-23 

solutions  are  obtained  for  flows  over  bodies  where  radius  r,  varies 

b 

as  power  of  the  axial  distance: 

n 

r,  -  cx  . 

0 

2 

Here  n  -  1  is  the  wedge  or  cor.e,  and  ~J~+y  <  n  <  1  yields  convex  bodies. 

Here  v  =  0  for  planar  bodies  and  v  -  1  for  axisymmetric  bodies.  The  special 

o 

case  n  -  yields  the  so-called  blast  wave  solutions  corresponding  to  the 

one  dimensional  unsteady  problem  of  an  intense  energy  release  along  the 
axis  of  symmetry  at  t  -  0.  For  axially  symmetric  bodies  the  blast  wave 

solution  co.  >sponds  to  a  parabolic  body  (n  =  l/2).  Blast  wave  theory  has 

q 

been  used  to  correlate  empirical  data  for  the  plume  shape  at  high  altitude. 

The  ease  of  application  of  the  theory  justifies  its  application  for  obtaining 


order  of  magnitude  results.  In  order  to  use  blast  wave  theory  for  the  flow 
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over  a  peraiolic  body  the  effect  of  the  blunt  nose,  where  ohe  slender  bod^ 

24,25 

assumption  fails,  must  be  considered.  xt  han  been  shown  that  nose 

blunting  of  otherwise  slender  bodies  can  have  a  disproportionately  large 
effect  on  the  flow  field  about  the  body.  Here  the  crong  and  highly  curved 
shock  at  the  nose  of  the  body  produces  a  strong  gradient  in  the  entropy  near 
the  body  surface.  Blast  wave  theory  may  possibly  be  used  for  the  flow  far 
from  the  body  if  an  approximate  account  of  the  entropy  layer  is  included 
(see  Sychev  and  "Jakuva  ').  When  the  similarity  assumption  may  not  be 
invoked  th.  partial  differential  equations  must  be  treated  directly. 

Normally  this  means  a  numerical  procedure  is  .,evuired. 

Numerical  Methods  for  Inviscid  Flow 
Supersonic  P-1  ow 

One  of  the  most  useful  methods  for  obtaining  numerical  solutions  of 
the  partial  differential  equations  for  supersonic  f Low  has  proven  to  be 
the  method  of  characteristics.  Prop  ly  used  this  method  yields  an  "exact" 
solution.  In  the  high  Mach  number  regions  of  tie  plume  flow,  however, 
difficulties  arise  in  the  ‘'nrlication  of  this  method  in  a  numerical  procedure 
Since  the  characteristic  lines  make  an  angle 


a 

the  characteristic  grid  is  elongated  in  the  streamwise  direction  when  the 
Mach  number  is  xarge.  As  a  result,  errors  will  be  magnified  in  the  stream 
direction  and  the  Inaccuracies  in  locating  tie  interaction  of  characteristics 
will  be  large.  To  overcome  these  inaccuracies,  small  grids  (i.e.,  small  in 
the  direction  transverse  to  the  streamlines)  would  be  necessary. 

In  the  plume  problem,  this  shortcoming  of  the  method  of  characteristics 
becomes  serious  many  nozzle  exit  diameters  from  the  nozzle  exit.  Fortunately 
the  regions  in  which  the  characteristic  angles  are  small  are  regions  in 
which  other  methods  (e.g.,  finite-difference  techniques)  become  quite  useful. 
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Under  tne  heading  of  finite-difference  techniques  we  include  all  numerical 
methods  for  direct  numerical  integration  of  the  partial  differential  equation 
other  than  the  method  of  cnaracteristics.  Basically  two  types  may  be  dis- 
tirguished:  those  that  use  a  prescribed  set  of  mesh  points  and  those  that 
include  the  determination  of  the  location  of  the  mesh  points  in  the  calculation. 
The  distinguishing  property  of  supersonic  flows  is  that  the  differential 
equations  are  of  the  hyperbolic  type.  Consequently,  finite  difference 
marching  techniques  may  be  directly  applied. 

Subsonic  Flow; 

Nose  Region  Solutions  for  the  External  Shock  Layer 

The  metnods  discussed  so  far  are  inapplicable  in  the  nose  region  of 
the  stream  shock  layer  where  the  flow  is  subsonic.  Here  the  elliptic  nature 
of  the  equations  prevents  the  direct  application  of  finite  difference  marching 
tecimiques.  For  the  present  problem  the  subsonic  region  is  expected  to  he 
small  and  to  exert  a  relatively  minor  influence  on  the  overall  flow  field. 

In  order  to  obtain  the  initial  values  for  the  calculation  in  the  supersonic 
region,  however,  an  estimate  of  the  subsonic  flow  in  the  nose  region  is  required. 

The  evaluation  of  the  hypersonic  flow  over  a  blunt,  body  has  been 
studied  extensively  during  recent  years  and  several  methods  for  inviscid  flow 
are  available,  including  the  Newtonian  solution,  •  the  thin-layer  approxi- 
mation,  *  1  tne  integral  relations  method,  and  the  various  numerical* 

methods  of  Swigart,  \  Van  Dyke/*'*  Garabedian  and  Lieherstein,  etc. 

Since  the  free  stream  is  hypersonic,  the  Newtonian  approximatioi  may  be 
used  to  specify  the  pressure  on  the  jet  boundary  in  the  subsonic  nose  region. 

This  provides  an  appropriate  boundary  condition  for  the  solution  of  the  jet 
now  in  the  inviscid  approximation.  Given  this  boundary  one  of  the  various 
approx iaiat ion  procedures  may  then  be  used  to  estimate  the  subsonic  flow. 
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IV-2.  The  In/iseid.  Flow  Field  Evaluation 

IV-2a.  Finite  Difference  Procedure  Used  For  the  Supersonic  Flow 
The  momentum  equation,  for  inviscid  steady  flow  is 


div  pu  -  0  .  (3) 

The  equation  of  state  is  p  ■  ffC£  . 

The  momentum  equation  may  be  rewritten  in  the  f  o.  m 


u2  is  +  u  |H  s  .  .  ESS_E  (It) 

ds  ds  p 

where  u  ■  un,  n  is  a  unit  vector  parallel  to  the  streamline  and  s  is  the 
distance  along  the  streamline.  For  two-dimensional  flow  we  may  write 


l  d£ 
p  ds 


0 


(5) 


and 

u2  |2  +  i  |E  »  0  .  (6) 

ds  p  dy 

Here  cp  is  the  angular  inclination  of  the  streamline  and  y  is  the  distance 
along  the  surface  normal  tj  the  streamline.  Equation  (6)  may  be  rewritten 


-TS  +  i  <h  +  iu2>-°-  co 

where  the  entropy  S  is  given  by  S  ■  R  [/r>  (p/p^)  ]/(y-l)  . 

Thus  the  flow  is  isentropic  along  streamlines. 


Here  h  r 
to  give 


Equations  2,  3,  7> 


may  be  integrated  along  streamlines 


P/p7 

puA 


■  constant 

■  constant 
constant 


P  "  pRT 

Here  A  is  the  area  of  an  individual  streamtube.  The  remaining  equation 
is  the  centrifugal  equation: 


where  we  have  set 


dcp  1  5 

38  V  37 


Un  p) 


(9) 


The  finite  difference  approximation  is  carried  out  as  follows.  The  flow 
is  divided  into  a  finite  number  of  streamtubes.  Initial  duta  is  fed  in  along 
an  orthogonal  surface.  The  initial  data  consists  of  the  local  streamtube 
properties  and  the  coordinates  of  the  dividing  streamlines,  on  the  input 
surface.  The  calculation  proceeds  by  predicting  the  flow  properties  on  a 
new  orthogonal  surface  using  the  flow  data  on  the  previous  surface  plus  the 
boundary  condition  at  each  edge  of  the  surface.  The  curvature  k(K,T'  of  the 
streamline  K  at  the  point  where  it  intersects  the  surface  L  is  evaluated  from 
the  lateral  pressure  gradient  determined  from  the  pressures  in  the  adjacent 
streamtubes  (K,  K  +  1)  and  (K  -  1,  K)  on  the  surface  L  (see  Figure  IV-3).  This 
curvature  is  used  to  extend  the  streamline  K  to  the  surface  L  +  1.  This 
procedure  is  repeated  for  each  streamline  in  the  flow.  From  the  thus 
determined  otreamtube  areas  on  the  surface  L  +  1  the  fluid  dynamic  properties 
for  each  tube  are  then  evaluated  from  the  three  conservatic  equations. 


IV -9 


It  is  known  that  the  straightforward  application  of  the  marching 
procedure  described  above  is  unconditionally  unstable.  However,  it  may 
be  rendered  '.table  by  the  following  procedure.  Having  determined  the  fluid 
properties  on  the  surface  L  -  1  he  curvature  h #  (fe,  L  *  1)  of  the  stream¬ 
line  K  at  the  surface  L  +  1  is  evaluated  from  the  pressure  gradients  on 
this  surface.  Then  the  surface  L  +  1  is  recalculated  using  for  the  mean 
curvature  of  the  streamline  K  between  the  surface  L  and  L  +  1  the  expression 

k  =  (l -a)  k(k.  l)  +  a  k'  (k,l+i). 

Perturbation  theory  (acoustic  approximation)  shows  that  this  procedure 
is  stable  when  Of  >  l/2>  is  enstabi-  when  ft  <  1/2  and  is  neutrally  stable 
when  a  =  l/2.  In  addition  to  this  restriction  the  step  length  between 
the  surfaces  L  and  L  +  1  must  be  chosen  to  be  less  than  the  value  jf  the 
separation  between  the  streamlines  K  and  K  ±  1  multiplied  by  a  number  of 
the  order  of  J  vf  -  1.  Calculations  with  various  test  flows  show  that  when 
ft  >  l/2,  perturbations  (sound  waves)  are  damped.  In  practical  calculations 
it  is  usually  convenient  and  sometimes  necessary  to  choose  a  slightly 
greater  than  l/2  in  order  to  smooth  cut  perturbations. 

As  presently  programmed,  the  procedure  marches  through  an  axisymmetric 
or  planar  feuw  generating  the  orthogonal  suri  .res .  On  each  surface  it 
calculates  the  flow  between  two  boundaries.  These  boundaries  may  be  of 
various  t.-pes:  a  shock,  a  solid,  or  a  free  (preesure  given  as  a  function 
of  boundary  inclination)  boundary.  The  shock  boundary  condition  uses  tne 
Rankine-Hugoniot  relations  to  connect  the  interior  flow  to  the  external 
flow.  The  external  flow  may  ’>e  non-uniform  and  ma  be  given  in  analyt ic 
or  tabular  form.  The  input  tabular  data  is  required  to  be  in  the  same  form 
a:-  the  data  the  program  calc  -late 
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In  the  evaluation  of  the  high  altitude  flow  field,  vc  first  calculate 
the  undisturbed  now  into  a  vacuum.  This  calculation  is  carried  sufficiently 
far  from  the  nozzle  that  an  analytic  extrapolation  to  infinity  may  be  carried 
out,  it  is  convenient  to  represent  this  flow  in  an  analytic  form  to  be  used 
as  the  flow  into  which  the  Jet-intercepting  shock  propagates.  A  procedure 
for  doing  this  is  outlined  in  Section  IV-i-.b. 
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IV-2b.  Represent at io a  of  he  Undisturbed  Jet  Flow 


The  flow  into  a  vacuum  is  relatively  easily  calculated  aid  there 

« 

are  a  number  of  formulations  available.  Here  viscous  effects  (except  for 
nozzle  and  combustion  chamber  boundary  layers)  are  expected  to  be  small. 

In  Figure  P/-2  we  show  the  streamlines  calculated  with  the  ,te 
difference  streamtube  program  for  the  50:1  sustainer  engine  expanding 
into  vacuum.  In  this  par  ticular  calculation,  the  flow  was  pre>'  -p  . .  Yom 
expanding  more  than  85°  by  an  artificial  wall.  At  the  present  ti.  e  che 
program  is  not  set  up  to  calculate  a  flow  which  turns  through  more  than 
90°.  Although  this  type  of  flow  car.  be  allowed  for  with  relatively  minor 
changes,  it  is  not  anticipated  that  this  will  be  necessary.  For  this  nozzle 
the  mass  flow  which,  in  a  vacuum  expansion,  turns  through  more  than  even 
50°  is  small  fraction  of  the  total  flow  (less  than  1$).  Thus,  in  this 
care,  we  do  not  expect  this  large  angle  flow  to  have  much  effect  on  the 
downstream  flow.  Also,  in  practice,  this  portion  of  the  flow  will  be 
largely  nozzle  boundary  layer  flow  and  is  not  well  defined  at  present 
(Fig,  IV-k),  (The  flow  that  turns  more  thar.  60°  is  confined  to  an  annulus 
less  than  1  mm  thick  at  the  nozzle  lip.) 

For  the  finite  difference  calculation,  it  is  convenient  to  express 
the  undisturbed  flow  in  an  analytic  form.  Since  the  flow  distribution  is 
singular  at  the  nozzle  lip,  the  zero  pressure  streamline  and  the  point  at 
infinity,  it  is  convenient  to  choose  an  analytic  form  that  gives  the  correct 
limiting  form  of  the  flow  at  these  singular  points.  A  useful  form  is  the 
following  Taylor  series  representation  of  the  densPy: 
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where 

po  is  the  exit  plane  density  a,  the  nozzle  edge, 

V  =rJ  r+VK->o) 
b  =  tan  G  r/x 

ft  =  the  turning  angle  for  the  zero  pressure  streamline,  and 
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x  is  measured  from  the  intercept  of  the  7  ~o  pressure  line  with  the  axis. 

_2 

This  form  gives  the  correct  dependence  c  £  far  from  the  nozzle  (p>  ~  r  ) 
and  the  correct  dependence  on  6  close  to  the  zero  pressure  line  (p  ~  (0^-  6 


The  coefficients  a 

mn 

Let 


may  he  evaluated  as  follows; 
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Then 
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cos  ff(n  +  ~)  <10.  (3) 


The  quantities  n^rj)  may  be  evaluated  numerically  at  various  values  of  77 
fiom  tne  vacuum  flow  field  calculation.  The  values  of  the  a^  may  then  be 
found  by  standard  least  squares  fit  procedures.  Standard  computer  routines 

are  available  for  this  latter  fit. 

A  procedure  similar  to  this  has  been  used  previously  and  was  found 
to  be  a  good  fit  to  the  flow  throughout  most  of  the  field  with  a  28  term 
series. 


2/y-1). 
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In  addition  to  tne  density  we  also  need  the  flow  inclination  angle 


cp.  Here  an  appropriate  representation  is: 

M' 

max 

l  1'“  cos !  G  C>!  ]  eo 

oo 

m=0 


The  same  procedure  as  that  outlined  above  may  be  used  here.  In  Figs. 

IV- 3  and  4,  we  show  the  comparison  between  the  numerical  calculation  of 
the  density  profile  and  streamline  locations  and  the  analytic  representa¬ 
tion  obtained  with  a  oxo  polynomial -Fourier  representation. 

The  outer  parts  of  the  vacuum  flow  field  are  appreciably  affected  by 
the  peripneral  mass  flows  associated  with  the  rocket  and  by  the  turbulent 
boundary  layer  within  the  nozzle.  In  Fig.  IV- 5  we  show  the  effect  of  the 
boundary  layer  on  the  vacuum  flow  for  the  50:1  engine.  (The  boundary 
layer  was  calculated  using  a  Von  Karm&n-Pchlhausen  method.)  In  Fig. 

IV-6  we  show  the  effect  of  the  sonic  exit  exhausterator  which  surrounds 
the  25:1  engine.  About  %  of  the  total  mass  flow  is  released  through  this 
exhausterator. 

_2 

In  the  far  field,  the  flow  from  any  nozzle  becomes  source-like  (o  ~  R  ). 
In  Fig.  IV-T,  we  show  the  mass  flow  angular  distribution  for  a  number  of 
different  nozzles.  In  the  fax  field  limit  the  analytic  representation  is 

2  f  Zli] 

a  Fourier  series  representation  of  H  p*  2  > .  In  the  present  calculation 
we  have  typically  used  a  5  or  b  term  sum.  For  approximate  purposes,  it 
often  suffices  to  use  a  one-term  representation: 
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dfi 


a  cos 
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In  Fig.  I V-8  we  compare  the  detailed  mass  flow  distributions  for  various 

nozzles  with  this  simple  form.  In  Table  IV -1  we  list  the  corresponding 

values  of  the  parameter  a.  Reference  to  this  figure  indicates  this 

form  gives  a  reasonable  representation  of  the  density  distribution 

obtained  (within  a  factor  of  2).  The  advantage  of  using  this  form  is 

that  it  may  be  evaluated  for  any  nozzle  without  requiring  a  detailed 

numerical  evaluation  of  the  vacuum  flow.  Mirels  and  Mullens  have  suggested 

[  2'  0 

a  similar  approximation:  vl  -  (8/0^)  instead  of  cos  ^  0”  *  Either 

00  I 

formula  gives  a  reasonable  rough  representation  of  the  numerical 
calculations. 

In  Figure  11-10  we  show  the  far  field  distributions  for  three  nozzles: 
the  50:1  nonuniform  exit  nozzle  and  two  25:1  nozzles,  #1  and  #2.  Nozzle 
§2  represents  the  far  field  flow  for  the  S-4  engine  whose  exit  plane  pres¬ 
sure  profile  is  shown  in  Fig.  IV-9-  This  distribution,  denoted  Nozzle 
#1,  was  intended  to  represent  this  flow  but,  due  to  a  mistake  in  the 
analytic  fit  procedure,  the  far  field  flow  distribution  actually  was 
incorrectly  computed  at  large  angles  to  the  axis.  To  a  certain  extent 
the  errors  are  of  a  n°ture  similar  to  a  boundary  layer  effect.  Unfortunate 1 
this  flow  distribution  was  used  to  evaluate  the  flow  fields  at  100  and 
135  km.  The  effects  on  the  external  flow  fields  are  expected  to  be 
relatively  small  and  confined  to  the  bow  region  of  the  plume.  Thus  we  do 
not  know  exactly  what  shape  nozzle  the  flow  labeled  corresponds  to. 

Howe"er,  differences  In  the  flow  are  probably  less  than  the  absolute 

accuracy  of  the  overall  calculation. 

T" - 

The  two  nozzles  shown  in  Fig.  II -10  have  a  sufficiently  non-uniform  exit 
plane  profile  that  a  precise  calculation  would  predict  a  significant  Mach 
disc  shock  structure  near  the  nozzle  exit.  In  the  present  calculation 
this  difficulty  was  circumvented  by  choosing  the  mesh  size  large  enough 
that  no  shock  structure  was  apparent.  The  resulting  errors  are  expected 
to  be  confined  principally  to  the  flow  near  the  axis. 
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XV-2c .  APPROXIMATE  TREATMENT  OF  THE  NOSE  HWIftCID  SHOCK  LAYER 

The  inviscid  flow  abound,  a  blunt  body  moving  at  hypersonic  speeds  has 
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been  treated  extensively.  Hayes  and  Probsteir*  give  an  excellent  review  of 
most  of  the  techniques  used.  For  the  present  problem  the  subsonic  region 
is  not  expected  to  exert  a  major  influence  on  c.he  downstream  flow.  For  high 
area  ratio  nozzles,  the  air  shock  may  remain  attached  (i.e.,  no  subsonic  region) 
to  rather  high  altitudes.  Thus,  at  least  for  the  initial  evaluations,  it 
is  not  expected  that  a  sophisticated  treatment  of  this  region  is  necessary. 

Of  the  several  procedures  available  those  that  assume  a  simple  relation 
for  the  pressure  distribution  more  or  less  a  priori  and  use  the  conservation 
laws  to  roughly  define  the  subsonic  region  are  particularly  suitable  to  be 
used  with  the  marching  ;echnique  used  here  for  the  supersonic  flow.  Here 
the  evaluation  of  the  subsonic  calculation  may  essentially  be  reduced  to 
a  boundary  condition  at  air-exhaust  dividing  streamline.  One  procedure 
that  will  be  considered  is  outlined  in  the  following  paragraphs. 

Here  the  thin  layer  approximation  is  used  to  estimate  the  pressure 
gradients  normal  to  the  boundary  streamline.  In  this  approximation  and  away 
from  the  stagnation  point  the  pressure  at  a  height  y  above  (i.e.,  normal  to) 


the  boundary  is 


p(y)  =  PB  + 


-r 


2, 

pu  dy 


(1) 


where  K  is  the  bcundary  curvature.  Assuming  that  the  streamlines  are  approxi- 


nately  narailel  to  the  surface,  the  mass  flux  per  unit  area  is  given  by 


Here  G  is  h/RT  and  p  io  the  total  pressure  for  a  given  streamline.  This 
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total  pressure  may  be  expressed,  in  terms  of  the  shock  inclination  at  the 
point  where  the  streamline  enters  the  bow  shock.  Mass  conservation  requires 


dy  =  p  u  R  dJR/rpu 


where  R  is  the  radial  displacement  of  the  streamline  at  the  point  where 


it  enters  the  shock. 


y  =  y(R)  = 


Thus  we  may  write 
/n  .  1  In  . 


R'd  fi' 


R  1/G  1/G v 
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P(y)  =  K  Pa)  um  um  \\l-  1  |L  dR'  +  Pfi.  (5) 

Also  the  local  radial  displacement  of  the  streamline  is 

r  =  r  +  y  cos  9  (6) 

B 

where  tan  9  =  dr/dz. 

B 

At  eech  station  the  shock  location  is  determined  by  solving 

Rs  =  +  y(Rg)  cos  9  (7) 

where  y(R  )  is  given  by  Equation  4. 
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In  order  to  evaluate  the  pressure  p,  the  boundary  pressure  p  must  be 
evaluated  in  the  Newtonian  approximation.  A  somewhat  improved  value  may 
be  obtained  by  requiring  that  the  pressure  behind  the  shock  computed  from 
Equation  (5)  agree  with  the  value  obtained  from  the  Rankine-Rugoniot 
relations;  hcvevor,  this  has  not  been  incorporated  in  the  present  calculation. 

This  procedure  is  not  valid  near  the  symmetry  axis.  However,  tb.e  effect 
on  the  downstream  flow  involved  in  using  it  into  the  axis  is  expected  to  be 
small.  To  start  the  calculation,  we  need  a  value  for  the  shock  stand -off 
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distance  that  is  compatible  wl Lh  the  above  procedure.  If  the  Newtonian 
form  for  the  boundary  pressure  close  to  the  axis  is  assumed,  and  the  jet 
flow  evaluated  using  this  boundary  pressure,  then  an  average  radius  of 


curvature  R  for  the  nose  may  be  defined.  The  boundary  curvature  should 
o 


be  averaged,  over  radii  of  the  order  of  the  expected  stand-off  distance  f . 
Assuming  that  p/p(  =»  1  and  that  6-^  «  1,  Equation  k,  in  the  limit  of 
r/r  «  1,  may  be  reduced  to  the  form 


o  u  u 

CO  CO  P] 


For  hypersonic  flow 


6  =  y(°)  =  ^fv 


-§(^)y^Tv 


(8) 


(9) 


These  relations  give  a  stand-off  distance  that  is  compatible  with  the 
procedure  described  by  Equations  4  to  7. 

In  practice,  however,  a  simpler  procedure  is  more  convenient:  the 
stand-off  distance  is  first  assumed  zero  and  the  entropy  gain  through  the 
shock  calculated  from  the  slope  of  the  line  joining  the  axis  starting  point 
and  the  first  shock  point  given  by  Equation  1  (unless  this  inclination 
exceeds  90°).  If  this  latter  condition  occurs,  the  entropy  gain  is  set 
equal  to  that  for  a  normal  shock. 

In  using  this  procedure,  it  was  found  that  the  shock  front  tended  to 
be  unstable.  The  method  was  altered  to  a  mixed  marching-relaxation 
procedure  in  which  the  flow  is  calculated  a  number  of  times.  The  stream¬ 
line  total  pressure  and  the  streamline  curvature  K(r,z)  (assumed  now 
to  vary  linearly  across  the  shock  layer)  for  the  (N+l)th  iteration  are 
evaluate"1  from  the  Nth  iteration.  For  the  first  Iteration,  tl  e  shock  slope 
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is  assumed  to  be  the  same  as  the  boundary  slope  and  the  streamline  curva¬ 
ture  uo  be  the  same  as  the  local  boundary  curvature.  This  procedure 
appears  to  converge  within  3  or  4  iterations  when  used  with  typical 
boundaries.  It  may  be  convergent  only  in  an  asymptotic  sense  in  that 
it  is  not  clear  that  small  irregularities  in  the  boundary  contour  do  not 
result  in  an  increasingly  fluctuating  shock,  profile.  However,  in  the 
present  calculations,  the  amplification  factor  is  sufficiently  small  that 
it  does  not  cause  any  difficulty.  (See  Figure  IV-10). 


IV-?d.  The  Mach  Disc  Region 


It  is  well  known  that  the  reflection  of  the  jet  snock  from  the 
symmetry  axis  is  of  the  irregular  type  wit n  triple  3hock  point  and  Mach 
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disc  formation.  The  detailed  solution  for  this  transonic  flow  requires 
the  solution  of  a  set  of  mixed  elliptic -hype  rbo.'.ic  equations.  Various 
approximate  criteria  have  oeen  suggested  to  determine  the  Mach  disc 

*V° 

location  in  a  simple  fashion.  Of  these,  that  suggested  ly  Bowyer, 
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D'Attorre  an.’.  Yoshihara  seems  to  give  the  best  agreement  with  experiment. 
Here  the  inclinations  of  the  Mach  disc  and  the  reflected  shock  are  calcu¬ 
lated  at  each  point  along  the  jet  shock  using  the  condition  that  the 
pressure  and  flow  direction  downstream  of  the  triple  point  are  the  same 
for  streamlines  passing  on  either  side  of  the  triple  point.  The  actual 
triple  point  is  supposed  to  occur  at  the  station  where  the  Mach  disc  is 
locally  normal,  to  the  incoming  flow.  In  other  words,  the  flow  passing  Just 
inside  the  tripie  point  is  assumed  to  s\tffer  only  a  slight  deflection 
(see  Figure  IV-U) . 


Fig.  IV-11.  Mach  Disc  Geometry 
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When,  as  in  the  high  altitude  jet,  the  Mach  number  in  front  of  the 
shock  is  high  (~  100),  this  criterion  for  the  triple  point  formation  may 
be  expressed  rather  simply.  In  this  hypersonic  limit,  it  may  be  shown 
that  the  Mach  uisc  would  form  when  the  wave  angle  of  the  jet  stack  reaches 
a  critical  value  (0  )  which  depends  only  on  the  specific  heat  ratio: 

( 1  -  e(0*)  j  -1'  sin  2e*  =  (y+1  - 2  si“2  e*)  z(0*)  - 1 


where 


z(e  ) 


(l - sin2  0*1  / 

\  (yi)2  / 


and 


e(9*)  = 


yi) 


4  2 

sin  9  . 


For  y  =  1.3,  9  takes  the  v«.lue  19 .2  .  For  (y-*-l)/(y-l)  »  1,  9  may  be 
approximated  by  the  much  simpler  expression 


# 

9 


v  y+1 


For  y  =  1.3,  this  approximate  relation  yields  a  value  of  9  close  to  the 
exact  one  (20.7  )* 

The  Mach  numbers  in  the  neighborhood  of  the  Mach  disc  indicate  the 
nature  of  the  flow  in  this  region.  Behind  the  incident  shock,  the  Mach 
number  is  7*3,  behind  the  reflected  shock,  4.64,  and  behind  the  Mach  disc, 
0.32.  The  pressure  behind  the  Mach  disc  is  about  20*  above  the  ambient 
value . 
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We  expect  that  this  criterion  (the  Mach  disc  locally  normal  to  the 

flow)  will  somewhat  overestimate  its  size.  Since  the  subsonic  i'low  must 

contract  in  order  to  lower  its  pressure,  the  streamline  passing  just  below 

the  triple  point  should  be  deflected  slightly  inwards.  If  the  slip  line 

o 

were  renuired  so  have  a  downward  inclination  of  only  0=)  ,  the  Mach  disc 
radiu  .  would  be  reduced  by  about  40$.  Although  D'Attorre  nd  Harshbarger 
have  compared  the  locally  normal  criterion  to  a  variety  of  experimental 
data  and  conclude  that  toe  triple  point  occurs  close  to  this  condition, 
the  sensitivity  of  the  size  of  the  disc  to  the  particular  criterion  chosen 
would  make  a  more  detailed  .=  tudy  of  the  flow  in  this  region  desirable. 

In  any  case.  Ute  disc  diameter  appeal’s  to  be  sufficiently  small  in  the 
present  cases  treated  t..at  its  presence  can  probably  he  ignored  for  most 
practical  purposes.  Reference  to  Figures  II-l  and  II-2  shows  that,  at 
180  km  and  M  4,3,  the  disc  diameter  is  less  than  one-tenth  the  maximum 
Jet  diameter.  For  these  highly  underexpandea  jets,  the  scaling  laws  show 
that  the  diameter  of  the  Mach  disc  (for  a  given  nozzle)  depends  only  on 
the  values  of  the  specific  heat  ratios  and  the  value  of  the  free  stream 
Mach  number  (M^) ,  Since  4,3  tends  to  be  <_  lower  limit  to  the  Mach  number 
during  sustainer  flight,  the  Mach  disc  will  be  even  less  important  at 
other  altitudes  (see  Figures  II- 4  and  t>) . 

These  conclusions  are  subject  to  the  qualification  that  the  vacuum 
flow  near  the  axis  is  properly  represented.  For  highly  non-uniform  nozzles 
the  near  axis  flow  nay  be  improperly  evaluated  in  the  present  calculation 
due  to  the  smearing  cut  of  the  Mach  disc  structure  near  the  nozzle  exit 
plane.  A  more  re- fistic  evaluation  of  this  part  of  the  flow  would  be 
required  to  evaluate  in  detail  the  influence  of  the  nozzle  non-uniformities 
on  the  down,  tream  Mach  disc  structure. 
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IV-3.  Viscous  Effects 


Since  the  altitude  range  of  interest  for  the  preset',  study  is  100 
to  300  kin,  the  viscous  and  non-continuum  effects  may  have  a  significant 
influence  on  certain  region)  of  the  flow  •f’ield.  Hie  free  stream  mean- 
free-path,  for  instance,  is  about  the  length  cf  a  typical  vehicle  (20m) 
at  135  km  altitude.  Hius  the  vehicle  itself  is  in  free -molecular  flow 
over  most  of  the  altitude  range.  Hie  Reynolds  number  and  Knudsen 
number  based  upon  free  stream  conditions  and  the  plume  nose  radius  vary 
over  quite  large  ranges . 

42. 

Hie  photographic  measurements  of  Rosenberg  and  his  coworkers 
indicate  that  the  outline  of  the  visible  plume  is  roughly  parabolic  in 
shape  and,  at  an  altitude  of  about,  180  km,  the  radius  of  curvature  of  the 
nose  region  for  an  Atlas  or  Titan  vehicle  is  about  1  km.  For  a  vehicle 
velocity  of  4  km/sec  at  this  altitude  the  Reynolds  number  based  on  the 
free  stream  velocity  and  density,  the  plume  nose  radius  and  the  viscosity 
at  the  stagnation  temperature  (—  9000°K  for  no  shock  induced  dissociation) 
is  of  the  order  of  30*  (Here  have  assumed  a  square  root  variation  of 
viscosity  with  temperature). 

Cheng  h*iS  treated  the  evaluation  of  the  flow  field  in  the  nose  region 

Z 

of  a  hypersonic  solid  blunt  body  at  low  Reynolds  number.  He  has  developed 
a  chin  two  layer  model  for  this  calculation  in  which  the  evaluation  of  the 
shock  layer  is  carried  cut  independently  of  the  evaluation  of  the  "shock 
transition"  region  and  the  two  profiles  matched  at  their  common  boundary. 
At  very  high  Reynolds  numbers  the  "shock  transition"  region  may  be 
identified  as  the  shock  itBelf  and  the  shock  layer  as  the  inviseid  flow 


i 
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behind  it.  fi,1  low  Reynolds  number  (merged  flow)  the  separation  is  a 
mathematical  artifice. 

In  Fig.  x~3  we  have  reproduced  Cheng's  Figure  5*3  in  which  he  compares 
the  thin  two  layer  evaluation  of  the  flow  with  that  obtained  by  Levinsky 
and  Yoshitiara  along  the  symmetry  axis .  Except  for  the  wall  boundary 
conditions,  the  flow  conditions  ai  similar  to  those  expected  at  about 
180  km  altitude  for  nose  regions  of  the  plume  of  an  Atlas  or  Titan  vehicle. 
Reference  to  Fig.I-c.  and  to  Fig.  1-3  shows  that  viscous  effects  will  be 
quite  important  in  the  nose  region  of  the  plume.  Although,  as  Cheng  points 
out,  the  good  agreement  between  the  two  calculations  is  not  an  indication  of 
the  absolute  accuracy  of  Cheng's  method,  it  does  give  us  confidence  in 
extending  the  thin  two  layer  procedure  to  the  evaluation  of  the  flow  field 
induced  by  the  rocket. 

The  importance  of  viscous  effects  far  from  the  nose  region  is  more 
difficult  uo  evaluate.  In  the  following  paragraphs  we  go  into  seme 
detail  to  get  an  estimate  of  the  magnitude  of  the  effect  since  this  pretty 
well  determines  the  approach  required  for  the  numerical  calculations. 

One  measure  of  the  importance  of  the  viscous  effects  is  the  ratio  6g/6  > 
where  6^  is  a  boundary  layer  thickness  -.ad  6  is  the  thickness  of  the  air 
shock  layer  assuming  inviscid  flow,  xf  p(r)  and  u(r)  are  mean  values  for 
the  density  and  velocity  across  the  layer  at  the  point  where  the  boundary 
is  displaced  radius  r  from  the  axis,  then 


In  order  to  estimate  the  value  of  6/r  we  note  that,  in  the  supersonic 

-  /C-l 

portion  of  the  flow,  u  varies  between  -J-3— s-  u  and  u  where  u  is  the 
y  ’  Y+l  mm  m 


rr-2k 


limiting  velocity  u  for  hypersonic  flow).  For  y  =  1.*+,  0.41  <  u/u^  <  1. 

The  mean  density  may  be  estimated  as  follows.  The  density  at  any 
point  in  the  inviscid  air  layer  may  be  written  in  the  form 

P  =  P  (p/ps)1^  (1) 


where  p  is  the  local  pressure  and  the  subscript  s  refers  to  conditions  on 

the  given  streamline  just  behind  the  air  shock.  When  the  air  shock  is 

strong  (i.e,,  when  M  sir.1' 9  »  1  where  tan  9  (=*  dr/cLz)  is  the  boundary 

00  B  B 

slope) 

2  ,  2a 
p  ~pu  sin  9_ 
rs  ■  ®  ®  B 


and 

Rosenberg's  field  observations 
roughly  parabolic : 


P 


p 

s  y-1  at> 

indicate  that  the  plume  bo'indary  is 


(2) 


rB(z)  ~/2Tp  (3) 

where  R  is  the  nose  radius.  Thus,  at  the  station  z,  the  shock  layer 
n 

pressure  is  of  the  order 


P 


p  u '/  (1  +•  r 
00 '  ' 


-IK)- 


(M 


As  a  mean  value  of  the  initial  pressure  p  we  choose  the  value 

s 

Ps  ~  +  I  (r/RN)2)‘  (5) 

This  is  the  value  of  p  for  the  median  streamline  (on  a  mass  flow  basis) . 
Thus,  in  the  supersonic  region  (roughly  r  >  Rq), 

|  *  p/ps  *  |  (6) 

or  3«65  k  p/p  k  4.9  for  y  =  1.4. 
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Near  the  sonic  line,  therefore,  pu  ~  2.0  p  u  whereas  far  downstream 

7  7  cd  on 

ou  ~  3*6  o  u  .  These  values  apply  only  to  the  region  where  the  shock  is 
strong.  When  M  sin  0  =*■  1,  the  shock  will  begin  to  diverge  from  the 

exhaust  plume  and  the  shock  layer  width  wi1!  begin  to  grow  linearly  with 
distance  downstream.  Also  here  the  parabolic  form  of  the  boundary  shape 
will  become  a  poor  approximation.  Assuming  a  parabolic  shape  so  long  as 
dr^/dz  >  1/M^,  we  estimate  the  value  of  the  boundary  radius  at  the  point 
where  the  shock  begins  to  diverge  from  the  exhaust  plume  as 


r'  ~  M  R  . 

sr>  n 

Thus  we  expect  the  preceding  approximations  to  be  useful  for  boundary  radii 

lying  between  R  and  M  R  . 

n  oo  n 

Choosing  an  average  value  of  2.8  p  u  for  pu  we  obtain  as  a  rough 

00  oo 

estimate  of  the  air  shock  layer  thickness  the  value 

6  ~  C.2r 


for  1  £  r/R  £  M  .  In  the  downstream  supersonic  portion  of  the  shock  layer 
n  a* 

we  may  roughly  expect  the  boundary  layer  thickness  no  grow  according  to 


6B  ' 


i.e.,  for  a  parabolic  plume  6,.  ~  r.  Thus  the  ratio  of  boundary  layer  thick- 
ness  to  shock  la^er  thickness  (6g/6)  appears  to  be  roughly  independent  of 
radius  for  R  £  r  £  M  R  .  For  r  >  M  R  ,  however,  we  may  expect  the  ratio 
of  the  boundary  layer  width  to  the  shock -boundary  separation  distance  to 
decrease  roughly  according  to 


V* 


Another  measure  of  the  importance  of  viscous  effects  is  the  ratio  cf 
the  shock  transition  zone  thickness  to  the  shock  layer  width.  For  a  normal 
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V 


shock,  a  measure  of  the  shock  thickness  is 


A*  = 


V  -Y 
2  1 


'S  "  ( dv/dy) 


max 


(7) 


1  2 

For  Pr  =  Le  =  1,  h  +  ^  v  is  constant  through  the  shock  and  dv/dy  is 
maximum  at  the  sonic,  point.  In  this  case  the  shock  thickness  may  be 
expressed  in  the  form 

* 


(8) 


Here  p  is  the  viscosity  at  the  sonic  point  and  =  v^/a  where  a*  is  the 
sound  speed  when  M  =  1.  For  an  oblique  shock  having  a  wave  angle  6> 
(assuming  p  ~/T)  the  shock  width  may  be  expressed  in  the  form 


is~c 


M  sin  0 
00 


sin2  0  /mJ+1 

CTi 


(9) 


Here 


Ml  = 


M  sin  8 
00 


V1  *  U ^ 


sin  0-1) 


and  C  is  e  constant  dependent  on  the  free  stream  conditions  Cc  = 

V  y+1 

In  a  Mach  8  free  stream  with  y  =  1.40,  the  shock  transition  widths  and 
the  sh'  ck  layer  widths  sire  compared  in  the  following  table  for  a  parabolic 
shock  (tan  0  =  dr/dz  ~  r^/r): 


■** 


"For  air  with  y  =  1.4,  C  =  5.9  X  10~^/p  cm  where  p  is  in  gms/ 


cm 


(at  180  km,  C  ~  104  cm). 
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R/Rn 

^  sin  9 

* 

M1 

At  IbO  km  with 

R  =  1  km 
n 

6 

km 

AS 

km 

,AS, 

(r) 

0 

8,0 

2.71 

.16 

.10 

0.62f 

l.C 

5.69 

2.27 

.2 

.12 

0.62 

2.0 

3.58 

2.07 

.4 

.15 

0.38 

4.0 

1.94 

1,60 

.8 

*29 

0.37 

Thus  it  appears  that  the  ratio  of  shock  transition  width  to  shock  layer 

width  varies  relatively  slowly  in  the  strong  shock  region  (1  s  r/r  <  M 

'  n  cd  ' 

In  the  outer  region  (r/r^  >  M^),  where  the  shock  is  weak,  the  shock 
width  increases  proportional  to  — j—  .  In  other  words,  when  the  shock 
i 

is  weak,  Ag  ~  —  where  Ap  is  the  pressure  jump  across  the  shock.  Far  from 
the  exhaust  gas  structure,  the  wave  is  expected  to  develop  an  N-wave  form. 
The  strength  of  the  trailing  shock  of  the  N-wave  depends  on  the  extent  to 
which  the  boundary  pressure  drops  below  ambient.  The  amplitude  of  the 
wave,  neglecting  the  viscous  effects,  decreases  according  to 


whereas  the  characteristic  width  of  the  wave  form  increases  as  R1^ 
(Whitham,  1952;  Bethe,  1948) .  Thus  the  shock  thickness  varies  as  r3/4 
and  the  ratio  of  shock  thickness  to  N-wave  width  increases  as  /r?  At 
sufficiently  large  radii,  therefore,  the  N-wave  structure  will  be  engulfed 
by  the  shock  transition  region.  At  larger  radii  the  wave  appears  as  . . 
smooth  acoustic  pulse  traveling  at  the  local  sound  sp?ed.  If  the  viscous 

*For  the  nose  region,  we  estimate  6  as  0.16  R^  (see  Section  IV-2c). 
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attenuation  were  negligible,  the  wave  shape  would  be  maintained  and  the 
amplitude  would  decay  as  R-1/2.  However,  because  of  the  viscosity,  the 
wave  is  both  distorted  and  attenuated.  Since  the  attenuation  depends 
sti'ongly  on  wavenumber,  the  change  of  the  wave  shape  is  best  examined  by 
considering  its  Fourier  components,, 

For  weak  sound  waves  the  viscous  equation  of  motion  may  be  written 

in  the  form. 

ii  -  <.2  is.  „  iat  i-  1 do) 

where  6  is  the  local  displacement.  For  5  =  6Q  ^  we  have  the 

dispersion  relation; 


Assuming  the  viscous  term  is  small,  we  may  write 

k  ~  k  _  i  k2 
-  o  3pc  o 

where  k  =  uj/c.  Thus  the  absorption  coefficient  is 
o 

l"1  =  ikik2/3pc  =  V(2w)2/3pc1q 


where  X  is  the  wavelength  of  the  disturbance.  An  examination  of  the  geometry  of 
o 

N-wave  formation  indicates  that  an  N-wave  wavelength  (\Q)  of  the  order 

of  M  may  be  expected.  The  numerical  calculations  of  the  flow  at  l6o 

km  for  M  =  4.3  indicate  that  an  estimate  of  about  4  times  MR.  is  more 

realistic.  At  180  km  this  correspords  to  about  20  km  for  a  typical  Atlas 

or  Titan  vehicle,  implying  an  N-wave  damping  distance  of  about  100  km. 

(9) 

Since  considerations  of  the  scaling  of  the  plume  shape  with  al  ude 


(see  also  Section  II -2)  indicate  that  the  term  pJyT  is  expected  to  vary 
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little  with  altitude,  the  N-wave  dumping  distance  is  also  roughly  independent 
of  altitude.  Note  that  these  considerations  really  only  apo’y  tc  j,  vertical 
trajectory  where  atmospheric  str unification  effects  are  not  too  important. 

At  these  altitudes  the  atmospheric  scale  height  is  about  ^0  km.  For  non¬ 
vertical  trajectories  the  mode  of  decay  of  the  N-wave  is  more  complicated. 

The  motion  in  the  horizontal  plane  will  be  roughly  similar  to  that  in  a 
nor-stratified  atmospher  but  propagation  in  the  vertical  direction  will 
be  distorted  by  the  atmospheric  stratification.  In  any  case,  at  very  high 
altitudes  where  the  plume  dimensions  become  appreciable  compared  to  the 
N-wave  damping  length,  N-wave  formation  is  not  expected.  At  the  lower 
altitudes  with  which  we  are  primarily  concerned,  a  disturbance  will  propa¬ 
gate  away  from  the  plume  but  the  high  frequency  components  will  be  largely 
absorbed  within  a  few  tens  of  kilometers  from  the  plume  axis.  Only  the 
low  frequency  components  will  persist  to  large  radii. 


( 


IV- ja.  THE  SUPERSONIC  MERGED  IAYF.H  CAICUIATION 

In  the  two  layer  flow  model  the  flow  field  is  divided  into  two 

X. 

regions:  the  shock  layer  and  the  shock  transition  zone.  Following  Cheng, 
we  denote  the  interface  separating  these  regions  as  the  "shock  interface." 
In  general  the  shock  layer  is  the  region  between  the  air  shock  transition 
zone  and  the  jet  shock  transition  zone.  We  first  consider  the  equations 
of  motion  in  the  shock  layer. 

The  basic  assumption  (equivalent  to  the  thin  layer  assumption)  is 
that  heat,  momentum  or  species  diffusion  parallel  to  the  streamlines  may 
be  neglected  compared  with  the  lateral  transport  (j.e.,  »  |j).  With 

this  assumption  the  Navier-Stokes  equations  may  be  reduced  to  the  form: 


div  pu  =  0 


u 


u  + 1 1£  =  i  .£(;!£) 

os  p  os  p  oy  V  oyy 

H  =  ;r sIC*  s| tE+Cs-  -  0 h] +  p  rdf) 


(1) 

(2) 

(3) 


u2|£=  .  ±|£  +  u|£) 

5s  p  oy  3p  dy  V  oy/ 


(«0 


Equation  2-2  may  be  written  in  the  form 


T  M  =  SH  _  1  _d  f  du\  ( 

b s  ds  p  dy  v  byJ 

When  the  viscous  effects  are  significant  only  the  continuity  equation  may 
be  directly  integrated  along  streamlines: 

puA  =  constant.  (6) 
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Equations  3>  and  6  plus  the  equation  of  state  form  a  convenient  set 

for  the  finite  difference  procedu  e.  Here  the  marching  procedure  is  to  be 
carried  out  in  a  manner  similar  to  the  invibcid  case  except  that  Eqr.  3 
and  5  '  re  used  to  evaluate  the  change  of  the  stagnation  enthaLpy  and  the 
entropy  along  streamlines. 

For  the  present  problem  ve  will  assume,  at  least  initially,  that  the 
diffusion  may  be  described  in  terms  of  a  single  binary  diffusion  coefficient. 
With  this  assumption  the  continuity  equation  for  each  species  is 


where  c..  is  the  mass  fraction  of  the  ith  species  and  on,  the  mass  rate  of 

"l  I 

formation  of  the  ith  species  per  unit  volume  per  second  due  to  chemical 
reaction.  In  the  present  study  chemical  reactions  are  expected  to  have  a 
negligible  influence  on  the  flow  field,  so  that  setting  =  0  is  expected 
to  be  a  reasonable  approximation.  With  these  approximations  Equation  7  is 
decoupled  from  the  other  flow  equations  and  thus  may  be  used  to  determine 
the  species  distribution,  but  is  not  required  to  determine  the  total 
density  and  temperature  distribution.  This  is  strictly  true  only  when  the 
local  mean  molecular  weight  does  not  vary  through  the  mixing  layer.  How¬ 
ever,  since  the  purpose  of  monitoring  the  species  diffusion  is  primarily  to 
determine  the  mean  specific  heat  and  the  air  partial  density,  this  approxi¬ 
mation  is  considered  to  be  a  reasonable  first  approach. 
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IV_3b.  THE  MODIFIED  RAJKINE-HUCrOHIOT  RELATIONS 

The  shock  transition  zone  equations  (given  in  the  following  section) 
when  evaluated  at  the  shock  interface,  yield  a  modified  net  of  RurJcine- 
Hugortot  equations  which  are  to  be  used  as  the  boundary  conditions  for 
the  calculation  of  the  shock  layer.  When  the  external  flow  is  a  uniform 
stream  having  a  density  pro,  a  pressure  velocity  and  stagnation 
enthalpy  H^,  these  modified  equations  have  the  form: 


pm  Uo  sin  *[ug  cos( -  cpp )  -  cos  iji  j  -  ~  [d  cos(Y  -  o)  )j 

C 

P2  -  pffl  +  Pro  U,  sin  sin(t  -  <*>2)  "  u«»  sin  ♦] 

=  ^  '|(u  Sin(t  ■ 


U 


2 

sin  t[H2  -  Hj  -  -  [d  j|  (»  +  +  5"  Bin2^  "  'f’Olg 


(1) 

(2) 

O) 

(M 


and 

p  U  sin  f  p0  u  sinU  -  q>?)* 

Here  tan  *  is  the  slope  of  the  shock  interface  and  the  subscript  2  refers  to 
conditions  in  the  shock  layer  evaluated  at  the  shock  interface.  For  a  non- 
uniform  str<^  such  as  the  internal  Jet  flow,  appropriately  modified 
equations  a,  ,quired.  In  the  transport  terms  of  these  equations  some 
factors  of  the  order  of  cos(\jf  -  q>2)  have  been  set  equal  to  unity  (here 
\jf  -  cp2  is  the  angle  between  the  shock  and  the  local  f  an'-ine  direction). 
Thus  these  equations  are  useful  for  shock  angles  such  that  the  angle  (  $-<£>) 2 
is  small.  This  latter  condition  is  not  too  restrictive  in  the  downstream 
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flow  when  the  fiee  stream  h  hype r sonic  since,  for  example,  fcr  y  =  1.4 

and  =  6.0,  V  -  o0  is  less  than  19°  when  the  Mach  number  Dehind  the  shock 
is  greater  than  l.y. 

When  m,  is  set  equal  to  zero  in  Equations  1  to  5,  they  reduce  u  the 
normal  Hugoniot  relations. 
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IV-3c.  THE  CHOCK  TRANSITION  ZONE  CALCULATION 

The  basic  approximation  made  in  the  two  layer  model  is  that  diffusion 


fluxes  are  significant  only  in  one  direction,  which  is  essentially  speci¬ 
fied  to  be  one  of  the  coordinate  directions.  In  the  supersonic  shock 
layer  this  direction  is  assumed  to  be  the  orthogonal  to  the  streamlines. 

In  the  shock  transition  region  the  appropriate  direction  is  normal  to 
the  shock  interface.  In  order  -chat  a  reasonable  approximation  be  obtained 
these  two  directions  should  differ  by  only  a  small  angle  at  the  matching 
boundary,  the  shock  interface.  In  the  nose  region  the  shock  layer  co¬ 
ordinates  are  chosen  as  the  normals  to  the  body  surface  and  the  associated 
orthogonal  surfaces.  The  small  angle  condition  requires  that  the  shock  and 
the  body  surface  be  approximately  concentric.  This  will  normally  be  satis¬ 
fied  in  hypersonic  flow.  In  the  downstream  supersonic  flow  we  require  that, 
at  the  shock  interface,  the  streamline  normal  makes  only  a  small  angle  with 
the  interface  normal.  This  condition  is  usually  approximated  in  hypersonic 

flow.  Far  downstream,  where  the  shock  is  weak,  the  angle  between  these  two 

_  2  2 
normals  approaches  (M^)  .  In  the  intermediate  region  where  Mw  sin  t  »  1 

but  ('  is  not  too  large  ($  is  the  wave  angle  of  the  shock)  this  difference 

angle  is  of  the  order  of  Thus  the  two-layer  approximation  is  expected 

to  be  useful  for  hypersonic  flow  for  small  values  of  (y-l)/2. 

For  the  shock  transition  region,  Cheng  chooses  a  coordinate  system  in 
which  y  =  0  corresponds  to  the  shock  interface.  Here  gradients  parallel  to 
the  shock  interface  are  to  be  neglected.  With  these  approximations  we  may 
write  (Cheng's  Equation  (5.7)) 

•|:(p  v)  =  0  (1) 
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T - 


I 


P  V 


§p. 

dy 


m  _ 

ay 


uU  d  f  Su  \ 

c  v  5y  ~  ay  ^  dyV 

(2) 

dv  4  9  f  dv 

+  p  v  ^  =  3  ay  V  *) 

(3) 

4(*  |iE  +  (r  ■  0  h  +  r  J ) 

(M 

Here  v  is  the  velocity  normal  to  the  shock  interface  and  y  the  distance 
along  the  normal. 

The  free  stream  boundary  conditions  are:  P1  =  p^,  =  oa,  h.,  =  h^ 

cos  t  v^  =  Um  sin  t  where  tan  $  the  local  slope  of  the  shock 
interface . 

Integrating  these  equations  across  the  shock  transition  zone  we  get 
for  arbitrary  strength  shocks  the  relations 


P  v  =  Pi  \  (5) 

p  f  "i  vi  T  ’  5  “  S  '  pi  +  pi  Ti2  (6) 

C1  V1  u  '  p  §  '  P1  T1  U1  (7) 

,  .  2 

pi  vi H  =  ^  &  CH  +  (r  ■ 1)h  +  ?" )  +  pi  vi  h  (3) 


This  set  of  ordinary  differential  equations  may  be  directly  integrated  to 
give  the  shock  transition  zone  structure.  Here  the  method  would  be  to 
start  at  the  shock  interface  with  the  values  of  the  flow  parameters  that 
were  obtained  from  the  shock  layer  calculation  and  integrate  outwards  along 
the  normal  to  the  shock  interface. 
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When  the  Prandtl  number  has  the  particular  value  3 /**  these  equations 
may  be  reduced  to  a  single  quadrature.  Since  expected  values  are  of  this 
order  (0. 70-0.72),  ve  ha/e  utilized  this  approximation  to  evaluate  the 
shock  profiles. 
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XV-3d.  the  nose  merged  iayeh 


This  is  the  problem  that  Cheng  treats  for  tile  rose  region  of  a  blurt 
body.  The  shock  transition  zone  equations  are  those  given  in  Section  IV-3c 
and  are  the  same  as  for  the  downstream  portion  of  the  flow.  There  are 
two  essential  differences  in  the  shock  layer  calculation  between  the  nose 
region  and  the  downstream  region.  First,  in  the  ncse  region  it  is  assumed 
that  the  diffusive  fluxes  occur  in  the  direction  par  illel  to  the  local 
normal  drawn  out  from  the  body  surface  instead  of  being  normal  to  the  local 
streamline  direction.  Second,  since  the  flow  equations  are  elliptic  in 
the  subsonic  region,  the  momentum  equations  are  altered  to  convert  the 

equations  to  a  set  of  parabolic  equations.  The  resulting  set  of  equations 

.  3  Ik 

(given  by  Hayes  and  Prob stein  )  is 


^(p  u  r)  +  ^(p  v  r)  =  0 


S-okcuv,<u|h  +  v|)  =  K,|; 

|+pKcu^<ug  +  v|)=^|) 


(1> 

(2) 

(3) 

(*) 


Here  u  is  the  x  component  of  velocity,  v  the  y  conq>onent,  y  the  distance 
along  the  outward  normal  and  x  the  distance  along  a  surface  parallel  to 
the  body  surface.  Kc  =  Kc(x)  is  the  curvature  of  the  body  surface  (of  the 
dividing  streamline  in  the  jet  case).  In  the  case  of  the  jet  we  would  take 
one  of  the  streamlines  which  originate  at  the  nozzle  lip  as  the  dividing 
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streamline.  For  the  flow  inside  this  streamline  we  use  the  su^croonic 
formalism  of  Sec.  lV-2h  and  the  above  formalism  for  the  flow  outside  this 
region.  An  appropriate  matching  condition  is  required  at  the  dividing 
streamline;  equal  pressure,  equal  and  parallel  velocity  and  equal 
transverse  fluxes,  ^rom  this  boundary  condition  would  come  the  equation 
that  specified  Kc(x) . 

Since  the  symmetry  axis  contains  a  singular  point,  the  stagnation 
point,  a  method  for  starting  the  matching  process  is  required.  For  the 
blunt  body  problem,  Cheng  expands  the  equations  in  a  Taylor  series  about 
x  =  0  and  solves  for  tie  zeroth  and  first  order  terms.  A  similar  procedure 
may  be  possible  for  the  jet. 

An  alternate  approach  is  to  solve  the  inviscid  Jet  flow  assuming  the 
dividing  streamline  pressure  is  that  given  by  Newtonian  theory.  Then, 
using  the  thus  determined  nose  rs.ius  for  the  exhaust  plume,  obtain  the 
starting  conditions  for  the  marching  process  usinc,  Cheng's  expansion 
procedure . 

The  change-over  from  the  subsonic  nose  calculation  to  the  supersonic 
downstream  calculation  does  not  have  to  be  made  right  at  the  sonic  point. 

It  should  be  made  downstream  at  a  point  where  the  angle  between  the  shock 
and  the  adjacent  internal  streamline  is  small  but  not  so  far  downstream 
that  the  streamline  curvature  varies  signif icantly  across  the  shock  layer. 
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IV-3e 


Small  Disturbance  Approximation 


An  approximation  which  has  proved  very  useful  in  the  treatment  of 
the  flow  around  hypersonic  slender  bodies  is  the  so-called  small  disturbance 
approximation.  Here  it  is  assumed  that  the  steady  two-dimensional  axi- 
aymmetric  flow  about  the  body  is  equivalent  to  a  one -dimensional  unsteady 
flow  that  would  be  seen  by  an  observer  at  rest  in  the  ambient  atmosphei'e . 

The  requirement  for  the  validity  of  tnis  approximation  is  that  the  perturba¬ 
tions  in  the  axial  component  of  the  velocity  be  negligibly  small.  If  a 
characteristic  deflection  angle  for  the  flow  is  t,  then  the  perturbation 

in  the  axial  component  of  the  velocity  is  of  the  order  of  u(l  -  cos  t)  or 

o 

.  In  order  that  this  velocity  perturbation  induce  only  a  slight  pressure 
variation,  it  should  be  small  compared  to  the  local  sound  speed.  In  other 

words,  the  small  disturbance  approximation  will  prove  useful  when 

2 

Mt  «  1 

where  T  is  a  characteristic  flow  inclination.  The  free  stream  Mach  numbers 
experienced  by  a  typical  missile  traveling  through  the  ionosphere  range 
between  4  and  12.  Thus  we  expect  that  the  small  disturbance  approximation 
will  prove  useful  when  the  flow  inclinations  are  less  than  about  5  degrees. 

In  air  the  shock  wave  from  a  vehicle  traveling  at  a  Mach  number  of  5  will 
give  rise  to  a  deflection  of  less  than  ?  degrees  at  distances  downstream 
where  the  pressure  jump  across  the  shock  is  less  than  a  factor  of  2.  Thus 
the  small  disturbance  approximation  will  be  strictly  valid  only  for  the  far 
field  part  of  the  flow.  For  the  near  field  the  more  complicated  two-dimensional 
steady  equations  are  required. 
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The  one-dimensional  flow  equations  may  be  written  in  the  form; 


1  D£  =  _L_ 

o  Dt  "  yJ 


(3) 


where 

H  =  h  +  |  u2  (M 

and 

P  =  p|h  "  \  )  *  (5) 

Here  J  =  0  for  planar,  J  =  1  for  socially  symmetric,  and  J  =  2  for 
spherically  symmetric  flows.  These  equations  may  be  recast  in  a  form 
more  convenient  for  numerical  evaluation! 


c  ,  t  3f  ,  f  i  . 

Dh  p  f  1  q  h  p  Sul 

Dt  v B 1 »/ ay  cp  p/ I  ™ 

(6) 
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Here  y  is  the  mass  point  coordinate  and  Ay  tne  spacing  between  two 
adjacent  points.  These  equations  have  been  programmed  for  the  IBM  7090. 
Representative  flow  distributions  calculated  for  a  simulated  sustainer 
engine  at  290  km  altitude  are  shown  in  Figures 


The  Air-Jet  Mixing  layer  -  Boundary  Layer  Approximation 

At  sufficiently  low  altitudes  the  mixing  layer  at  the  Jet-air 
Interface  will  remain  thin  compared  to  the  shock  layer  thickness  if  it 
remains  laminar.  In  order  to  estimate  the  importance  of  the  viscous 
dissipation  at  the  Jet-air  interface  a  much  simplified  model  has  Deen 
used  to  evaluate  the  mixing  layer  profile  at  low  altitudes.  Although  the 
mixing  layer  is  unstable  to  become  turbulent,  it  is  expected  that  the 
time  required  for  transition  to  take  place  will  be  long  compared  to  the 
flow  time  (through  the  major  part  of  the  jet)  at  least  at  the  higher 
altitudes . 

In  order  to  obtain  an  easily  applied  analysis  two  approximations 
have  been  made:  first;  that  the  mixing  layer  remains  sufficiently  thin 
that  the  variation  of  the  in-iscid  flow  field  properties  across  the 
mixing  layer  is  negligible;  second,  that  the  rate  of  change  of  the 
inviscid  flow  field  with  distance  downstream  is  sufficiently  slow  that 
the  mixing  layer  profile  may  be  assumed  to  be  self -preserving .  In  this 
approximation  it  can  be  shown  that-  the  air- Jet  mixing  layer  profile  is 
equivalent  to  that  between  two  parallel  and  uniform  streams .  It  turns 
out,  however,  that  even  for  the  lower  altitudes  of  our  interest  neither 
of  these  two  approximations  may  be  expected  to  be  very  good  -  particularly 
the  former.  In  Figure  IV-1  we  show  density,  pressure,  and  temperature 
distribution  in  the  neighborhood  of  the  boundary  calculated  for  the  25:1 
engine  at  135  ten  altitude  assuming  no  viscous  transfer.  Also  in 
Figure  IV-1  we  show  the  self -preserving  mixing  layer  profile  calculated 
assuming  that  the  conditions  at  the  jet -air  boundary  were  maintained 
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across  the  mixing  layer.  Reference  to  this  figure  shows  that,  especially 
in  the  Jet  shock  layer,  the  inviscid  pressure  and  density  both  vary 
strongly  across  the  mixing  layer. 

Because  of  the  strong  variation  of  the  temperature  and  density 
across  the  Jet  sr.ock  layer  the  calculated  self -preserving  profile  cannot 
be  expected  to  be  too  meaningful.  Techniques  ore  available  which  will 
allow,  in  an  approximate  manner,  for  some  variations  in  the  fiee  stream 
flow  outside  the  mixing  layer.  However,  since  it  appears  that,  at  least 
above  135  toe  altitude,  the  Jet  shock  layer  may  be  completely  merged, 
it  appears  that,  at  the  very  least,  a  viscous  layer  ci lculation  will 
be  required  to  evaluate  the  entire  shock  layer. 

Thus,  for  the  present,  it  does  not  appear  too  useful  to  try  and 
refine  this  mixing  layer  calculation  in  the  boundary  layer  approximation. 
Tne  description  of  the  analysis  of  the  mixing  layer  keeping  the  self- 
preserving  assumption  will  be  reported  subsequently.  In  brief  the 
procedure  is  to  treat  the  equivalent  mixing  layer  between  two  uniform 
free  streams.  Here  the  equations  admit  a  similarity  solution  in  the 
classical  compressible  boundary  layer  approximation.  The  Von  Karman- 
Fohlhausen  method  is  used  to  obtain  3n  approximation  of  the  solution  to 
these  equations.  The  procedure  has  been  programmed  for  the  IBM  7090 
computer. 

Bie  basic  conclusion  of  this  analysis  is  that  the  Jet  shock  layer 
will  be  strongly  affected  by  the  diffusive  fluxes.  In  our  present 
description  of  the  electron  distribution,  the  neutral  density  distribution 
within  the  Jet  region  may  exert  a  strong  influence  on  the  electron 
distribution.  Thus  it  appears  necessary  to  obtain  a  considerably  more 
sophisticated  description  of  the  Jet  flow  than  that  resulting  from  a 
simple  self -preserving  boundary  layer  approach. 
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I V - 3 & •  Evaluation  of  the  Initial  Flow  Distribution  at  Very  High  Altitudes  in  the 


Collin ionlebs  Approximation 

'Jlie  mixing  layer  calculations  carried  out  at  135  altitude  indicate 
that  the  jet  shock  Layer  may  be  completely  merged  even  at  this  relatively 
low  altitude.  Although  at  this  altitude  the  shock  widths  are  relatively 
thin  compared  to  the  shock  layer  thicknesses,  this  will  not  be  the  case 
at  very  high  altitudes.  At  sufficiently  high  altitudes  that  the  shock 
layer  width  is  comparable  to  a  mean  free  path,  the  flow  will  be  poorly 
described  by  the  continuum  equations.  This  situation  is  approached  in  the 
bow  region  at  altitudes  above  200  km. 

An  alternate  approach  to  a  complete  continuum  calculation  is  to  compute 
the  early  non-continuum  part  of  the  flow  assuming  zero  interaction  of  the 
exhaust  gases  with  the  ambient  atmosphere  (cross  sections  for  collisions 
between  exhaust  molecules  and  air  molecules  vanish),  and  then,  at  some 
more  or  less  aroitrary  time,  start  a  one-fluid  continuum  calculation.  This 
ic  a  particularly  useful  procedure  at  high  altitudes  from  a  numerical  point 
of  view,  since  a  large  amount  of  computing  time  is  normally  required  to 
compute  the  flow  close  to  the  vehicle  when  the  continuum  equations  are  used 
throughout . 

We  first  of  all  consider  this  collisionless  approximation  of  the  initial 
flow  in  the  small  disturbance  approximation.  The  small  disturbance  approxi¬ 
mation  gives  a  useful  descrip  .on  of  the  flow  field  when  the  axial  gradie'’fc 
of  the  various  flow  properties  are  small  compared  to  the  radial  gradients 
and  may  be  neglected.  If,  in  addition,  the  variation  of  the  axial  component 
of  the  velocity  is  small  across  the  entire  flow  field,  a  Galilean  transformation 
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may  be  used  to  transform  the  steady  state  two  dimensional  flow  to  a  one 
dimensional  unsteady  flow  problem.  We  restrict  the  present  discussion  t<o 
this  case. 

Suppose  the  vehicle  velocity  relative  to  the  atmosphere  is  V  .  We 

assume  therefore  that  the  axial  component  of  the  exhaust  gas  velocity 

relative  to  the  vehicle  is  also  V  .  At  high  altitudes  V  is  'therefore 

m  m 

related  to  the  enthalpy  release  in  the  combustion  chamber  by  Lhe  expression 


V 

m 


(1) 


where  hc  is  in  units  of  ergs  per  gram.  At  very  high  altitudes  the  exhaust 
flow  close  to  vehicle  is  source-like.  Thus  we  may  write  (see  Section  IV-2b) 


where 


__  m  f(cp)  sin  9 


m 


(2) 


(3) 


and 

2ir  Acp)  sin  cp  ip  =  1  9  (4) 

o 

r  is  the  r.vdial  displacement  from  the  symmetry  axis  and  z  tnu  axial  distance 
(?  w  V^t).  The  radial  component,  v,  of  the  velocity  is  sin  cp.  If  m  is 
the  total  exhaust  mass  f loc  rate,  then  m/V  is  the  mass  deposited  per  unit 
length  along  +he  -trajectory.  The  fraction  of  the  exhaust  flow  that  has  a 
radial  velocity  between  v  and  v  +  dv  is  given  by,  approximately. 


1  = 
m  dv  V 


an  9 


ra 


(5) 
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where  tin  cp  =  v/V  . 

T  '  m 

For  the  vacuum  expansion  the  radial  velocity  as  a  function,  of  radius 
and  tiii.-e  is 

v  -  r/t  (6) 


The  exhaust  gas  density  as  o.  function  of  radius  and  time  is  (for  R  <  V  t) 


P 


=  m  f(cp) 

3  2 
VJ  t 


where  cp  =  sin  ^  •  A  useful  approximation  for  f  (cp)  is  a  ^cos  tt,  cp/0^  ^"*^f or 

cp  <  9ro  (see  Section  IV' 2b)  where  9^  is  the  limiting  expansion  angle.  In 
this  small  disturbance  approximation  9m  is  assumed  less  than  unity. 

In  order  to  ir'tiate  a  continuum  calculation,  we  require  the  mean 
fluid  properties  to  be  specified  in  terms  of  a  one  fluid  model.  The  mean 
fluid  velocity  will  be  given  by 


v  -  p  v  /lP„+P  J 
e  e'  a  e 


(8) 


where  pa  is  the  ambient  density. 

The  fluid  density  is  given  by 

P  =  Pe+Pa  (9) 

and  the  enthalpy  by 

h  =  [Paha  +  \  Peva  -  i  P  v2]  i  .  (10) 

In  Fig.  IV-lj,  we  shew  the  collisionless  approximation  to  the  flow  field  at 
an  altitude  about  290  km. 

This  model  gives  an  accurate  description  of  the  flow  only  when  the 
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mean  free  path  for  exhaust-air  collisions  at  all  points  in  the  flow  field 
is  larger  than  the  flow  field  dimension.  Although  this  condition  is  not 
even  approached  in  the  rocket  flow  field  at  any  altitude  of  interest,  never¬ 
theless  this  approximation  provides  a  useful  representation  of  certain 
aspects  of  the  flow  while  retaining  a  very  simple  analytic  stricture.  In 
the  remainder  of  this  section  we  examine  the  validity  of  the  collisionless 
approximation. 

It  should  be  emphasized  at  this  point  that,  in  these  approximations,  we 
are  primarily  concerned  with  the  early  part  of  the  flow  (i.e.,  the  bow  region 
of  the  plume).  As  will  be  apparent  in  the  ensuing  discussion,  the  specific 
criterion  is  that  the  exhaust  gas  density  on  the  axis  of  symmetry  be  much 
larger  than  that  of  the  ambient  atmosphere.  In  order  to  get  a  better  under¬ 
standing  of  the  approximations  involved  we  will  first  consider  another 
simplified  but  more  accurate  description  of  the  early  part  of  the  flow. 

Here  we  assumt  that,  in  the  region  where  collisions  of  ambient  molecules 
with  exhaust  molecules  occur,  the  exhaust  gas  density  is  much  larger  than 
the  ambient  gas  density.  As  long  as  the  effect  of  the  ambient  atmosphere 
is  small,  the  exhaust  gas  density  rises  from  zero  at  the  edge  of  the  vacuum 
plume  with  a  high  power  (— j-)  of  the  distance  from  the  edge  (see  Section  IV-2b). 
Close  to  this  edge  the  exhaust  gas  temperature  is  near  zero  and  the  velocity 
attains  its  limiting  value,  V  .  In  a  coordinate  system  moving  with  the 
exhaust  gas  the  ambient  atmosphere  appears  to  impinge  on  the  stationary 

exhaust  molecules  with  a  velocity  V  . 

m 

Let  us  consider  first  the  sequence  of  events  which  occur  during  the 
expansion  when  the  molecular  weight  of  the  atmospheric  species  is  much  larger 
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than  that  of  tne  exhaust  gases.  This  would  apply  to  a  rocket  exhausting 
pure  hydrogen,  for  example.  In  this  case  the  atmospheric  molecules  penetrate 
deeply  into  the  exhaust  gases  losing  small  amounts  of  momentum  and  energy 
at  each  collision.  The  fractional  energy  loss  per  collision  is  of  the  order 
e,  where  e  is  twice  the  mass  ratio^and  the  average  scattering  angle  is  e/2. 
The  kinetic  energy  is  reduced  by  an  order  of  magnitude  after  about  1.2/e 
collisions.  The  rms  scattering  angle  at  this  point  is  of  the  order  ./e/2 
radians.  Thus,  when  the  exhaust  molecules  are  light,  the  atmospheric 
molecules  penetrate  the  exhaust  in  essentially  straight  lines  and  come  to 
rest  with  respect  to  the  exhaust  molecules  after  a  few  collisions .  Cie 
kinetic  energy  is  dissipated  in  heating  the  exhaust.  Until  such  time  that 
the  fractional  concentration  of  atmospheric  molecules  becomes  significant, 
the  increase  in  the  exhaust  gas  temperature  is  slight.  Because  of  the  strong 
rate  of  increase  of  exhaust  gas  density  with  distance,  the  bombarding 
molecules  will  al~  be  stopped  at  roughly  the  same  depth.  Thus,  in  this 
model,  we  expect  a  snowplow  action  in  which  ambient  mo1 .cules  are  swept 
up  into  a  thin  layer  located  some  distance  from  the  plume  edge.  Initially, 
the  density  of  ions  in  this  layer  increases  at  a  constant  rate.  This  process 
wild,  continue  until  the  deposited  energy  raises  the  temperature  of  the  mixture 
sufficiently  to  cause  an  expansion  of  the  gas  into  the  oncoming  flow.  This 
constitutes  the  development  of  the  shock  layer.  A  non-continuum  transient 
period  now  follows  as  the  shock  layer  becomes  fully  developed.  In  the  fully 
developed  state  we  expect  that  the  peak  electron  density  will  be  limited 
to  the  value  cor  sponding  to  a  strong  shock.  Thus,  according  to  this 
model  the  electron  density  may  go  through  a  maximum  as  the  flow  develops. 
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Thin  maximum  value  may  or  may  not  exceed  the  value  expected  behind  a  strong 
shock. 

The  opposite  situation  obtains  in  the  other  limiting  case  when  the 
molecular  weight  of  the  atmospheric  molecules  is  much  smaller  than  that  of 
the  exhaust  gases.  Here  again  there  is  little  energy  transfer  per  collision 
but  the  average  scattering  angle  is  large.  The  incoming  molecules  emerge 
from  the  exhaust  gas  structure  after  a  few  collisions  with  essentially  the 
same  energy  with  which  they  entered  but  are  traveling  in  random  directions . 

In  this  approximation  the  exhaust  gas  acts  like  a  rigid, diffusely  reflecting, 
adiabatic  wall.  Tlie  atmospheric  shock  layer  will  be  established  in  front 
of  it  in  a  few  coJlision  times.  In  this  limit ,  the  electron  density  is 
expected  to  rise  smoothly  to  its  maximum  value . 

In  cases  of  practical  interest  the  molecular  weights  of  the  exhaust 
gas  and  of  the  ambient  atmosphere  are  comparable.  Thus  we  expect  a  state 
intermediate  between  the  two  models  described  above.  Since  the  structure  of 
the  shcck  layer  is  considerably  different  in  these  two  limiting  models ,  it 
would  appear  that  a  fairly  detailed  description  of  this  non-continuum  regime 
is  required  to  obtain  a  reasonauly  accurate  description  of  the  flow.  Attempts 
in  this  direction  have  not,  in  general,  met  with  much  success.  Ihe  approach 
usually  invoked  is  to  assume  that  the  continuum  Navier-Stokes  equations 
apply  throughout  the  flow.  In  this  case,  even  though  the  distribution 
functions  are  not  Maxwellian,  a  one -fluid  model  is  used  in  which  a  unique 
temperature  and  stream  velocity  is  assigned  to  each  point  in  the  flow. 

The  collisionless  approximation  (defined  by  Eqs .  8,  9>  and-  10). 
is  equivalent,  near  the  plume  edge,  to  a  one -fluid  approximation  of  the 
heavy  atmospheric  molecule  case.  However,  this  simple  prescription  fails 
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~n  tlie  neighborhood  of  the  thin  layer  where  most  of  the  atmospheric 
molt  cyj.es  collect.  Rough  calculations  of  the  density  of  atmospheric 

molecules  in  this  layer  indicate  that  the  fractional  concentration  of  air 
here  is  of  the  or^er 


where  e  is  the  fractional  energy  loss  per  collision  (e  «  2  m/M),  y  is 
the  specific  heat  ratio,  R  the  exhaust  gas  radius,  X&  the  mean  free  path 
in  the  ambient  atmosphere,  n  the  number  density  in  the  ambient  and  n 

c 

the  number  density  o-  the  axis  of  symmetry.  This  expression  is  useful  when 
both  (y+1)/(y-1)  and  the  density  ratio  n  /n  have  large  values.  For  our 
conditions  we  expect  values  of  the  density  ratio  between  xcA  and  unity. 

For  y  =  1.4,  and  e  ^  0.3,  for  example,  the  fractional  concentration  in  the 
deposition  layer  is  then  of  the  order  0. I  to  0.5  times  (r/x  )5/6.  3hus 
we  expect  the  collisionless  approximation  to  be  useful  so  long  as  the  radius 
of  the  vacuum  boundary  is  smaller  than  the  mean  free  path  in  the  ambient 
atmosphere.  At  later  times  in  the  expansion  the  heating  of  the  exhaust  gases 
by  the  entrained  atmosphere  will  result  in  the  development  of  an  appreciable 
shock  layer  structure. 

This  procedure  has  been  used  to  provide  initial  profiles  for  a  small 
disturbance  calculation  of  the  flow  field  at  very  high  altitude  (290  km). 

In  a  small  disturbance  approximation,  the  flow  is  assumed  identical  to  an 


unsteady  cylindrically  symmetric  viscous  expansion.  A  solution  of  the  one- 
dimensional  unsteady  Navier-Stokes  equations  has  been  obtained  numerically 
using  a  finite  difference  procedure.  The  density,  velocity  and  pressure 
distribution  at  various  times  are  shown  in  Figs.  II-1.8,  19,  and  20.  The 
calculation  was  initiated  0.6  seconds  after  the  exhaust  gas  was  expelled 
from  the  nozzle  and  used  the  collisionless  representation  as  the  initial 
state.  At  this  time  the  radius  of  the  exhaust  plume  was  about  2  km.  Since 
the  mean  free  path  at  this  altitude  is  about  1.5  km,  the  collisionless 
approximation  is  expected  to  be  unreliable  after  this  time.  In  Fig.  11-21  we 
have  plotted,  versus  time,  the  radius  of  the  cylinder  containing  a  mass 
equal  to  that  of  the  exhaust  products.  Also  shown  in  Fig.  11-21  is  the  locus 
of  the  points  of  equal  exhaus t -a ii  mass  fraction.  The  dividing  streamline, 
which  bounds  a  mass  equal  to  that  of  the  exhaust  gases,  expands  rapidly  for 
20  seconds  at  which  time  it  attains  its  maximum  value  of  about  15  km.  After 
this,  it  contracts  lightly,  reaching  a  minimum  of  13  km  at.  50  seconds.  After 
hk  seconds,  the  exhaust  gases  are  so  diluted  with  air  that  nowhere  does  the 
fractional  air  concentration  drop  below  50$. 

Although  the  collisionless  representation  is  not  expected  to  be  valid 
when  the  plume  radius  is  greater  than  the  ambient  mean  free  path,  it  apparently 
gives  a  rough  approximation  of  some  features  of  the  flow  field  throughout  the 
early  part  of  the  expansion  even  when  the  radius  considerably  exceeds  the  mean 
free  path.  In  Fig.  11-22  we  have  plotted  the  velocity  and  density  profiles  at 
3.2  and  6.5  seconds.  References  to  this  figure  show  that  the  collisionless 
approximation  gives  quite  a  reasonable,  although  smoothed  out,  representation 
of  the  velocity  profile  even  at  these  relatively  late  times.  The  density 


distribution  is  less  well  represented,  especially  in  the  shock  layer  region. 
Nevertheless ,  since  the  gross  features  of  the  flow  are  reasonably  well 
reproduced,  the  very  simple  analytic  form  of  this  approximation  makes  it 
quite  useful,  for  eliciting  qualitative  and  semi  -quantitative  information 
about  the  interaction  of  the  ionized  species  with  the  neutral  wind  during 
the  early  part  of  the  expansion. 

So  far  we  have  only  considered  a  cylindrical  expansion.  It  is 
expected  that  similar  agreement  (and  lack  of  agreement)  may  be  expected 
in  the  more  general  axially  symmetric  case.  Even  the  general  three 
dimensional  problem  of  a  rocket  traveling  in  a  stratified  atmosphere  is 
trivially  formulated  in  this  collisionless  approximation. 


IV-53 


IV-ih.  Kinetic  Effects 


The  calculations  carried  out  at  180  km  indicate  that  the  high  tempera¬ 
ture  region  is  confined  to  e  Bmall  volume  in  the  bow  region  of  the  plume. 
Although  the  stagnation  temperature  of  the  air  stream  is  relatively  high 
(~  6500°K  at  180  km  altitude)  the  cooling  of  the  air  stream  by  the  outward 
diffusion  of  the  cold  exhaust  gases  maintains  the  gases  at  relatively  low 
temperatures.  The  high  temperature,  low  density  entropy  layer  that  would 
exist  adjacent  to  the  Jet  boundary  if  the  diffusive  effects  were  negligible, 
does  rot  occur  at  these  altitudes  (Bee  Figure  11-12). 

In  Figures  II -lb  and  17  we  show  the  tempera* >  -e  and  density  histories 
along  two  characteristic  streamlines.  The  streamline  locations  are  shown 
in  Figure  II -1.  Reference  to  these  figures  indicates  that,  even  for  fluid 
elements  wuich  enter  the  air  shock  layer  relatively  close  tc  the  stagnation 
point,  the  temperature  exceeds  2000 °K  for  less  than  a  second.  Since  the 
collision  frequency  in  the  nose  region  1b  less  than  100  sec-1  only  reactions 
having  fairly  high  rate  coefficients  and  low  activation  energies  can  be 
important. 

At  this  altitude  the  atmosphere  is  compared  primarily  of  and  0. 
Reference  HH  gives  the  mole  fractions  of  Ng,  0  and  0 ^  as  0.55,  0.39  and  0.06 
respectively.  Even  the  most  likely  ionization  reaction  between  the  air  con¬ 
stituents 


followed  by 


Ng  +  0  -  NO  +  N 


N  +  0  -  N0+  +  e 

will  produce  a  negligible  increase  in  the  electron  density.  The  mixing  of 
the  exhaust  constituents  with  the  ambient  atomic  oxygen  may  produce  signifi¬ 
cant  reaction. 


IV-51* 


\ 


Molecular  hydrogen  present  in  the  exhaust  gases  will  react  with  atomic 
oxygen  according  to 

Hg  +  0  -  OH  +  H 

i  p  3600  °K/T  3  ^ 

Kaskan  and  Browne  recommend  p  reaction  rate  of  3*3  x  10  e  J  **  enr/mol»-sec. 

Thus,  along  streamline  a  (see  Figure  II-l)the  fraction  of  the  hydrogen  that  reacts 
with  atomic  oxygen  is  estimated  to  be  one  or  two  per  cent. 

Observations  by  Rosenberg  indicate  the  CH  radical  is  present  in  the 
exhaust  products  of  Atlas  and  Titan  I  vehicles ,  It  may  form  part  of  one  of 
the  products  discharged  from  the  g68  generators.  Since  these  products  will  be 
concentrated  at  the  periphery  of  the  plume,  we  may  characterize  their  environ¬ 
ment  as  similar  to  that  on  streamline  a.  CH  can  give  rise  to  ionization  via  the 
reaction 

CH  +  0  -  CH0+  s  e . 

This  reaction  is  considered  to  be  a  source  of  chemionization  in  hydrocaii*on 

12  3 

flames.  Calcote  gives  a  value  of  about  2  x  10  cnr/mole-sec  for  the  rate 

M7 

coefficient . 

This  value  implies  that  10  to  20  per  cent  of  the  available  CH  radicals 
will  yield  free  electrons.  Hie  concentration  of  CH  in  the  exhaust  gases  in  not 
known.  In  the  At  "'as  engine  about  three  per  cent  of  the  total  mass  flow  is 
diverted  through  the  gas  generators.  Since  this  portion  of  fuel  flow  is  burnt 
at  very  rich  conditions  it  is  expected  that  a  variety  of  hydrocarbons  will  be 
present.  If  more  than  one  per  cent  of  this  mass  flow  results  in  the  CH  radical, 
the  resulting  contribution  to  the  electron  content  of  the  flow  field  may  be 
significant  at  this  altitude.  This  source  of  ionization  will  be  more  important 
at  lower  altitudes  and  less  at  higher  altitudes.  However,  it  is  difficult  to 
determine  its  absolute  importance  without  a  more  precise  definition  of  the 
exhaust  composition. 


IV-55 


Glass  quotes  a  value  for  the  i  ate  constant  of  acetelyne  with  atomic 
oxygen  at  low  temperatures  of  2  x  10  ^  cm^/sec.  Foi  this  value,  we  voujd 
expect  about  1^  of  the  available  acetylene  to  react  along  streamline  a. 
Since  only  a  fraction  of  the  flow  will  be  hydrocarbons  and  oifLy  a  fraction 
of  those  that  do  react  with  atomic  oxygen  will  yield  CH,  this  source  of 
chemiionization  may  rot  be  of  major  importance. 
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yL  POINT  BEING  CALCULATED 

O  POINTS  USED  IN  FIRST  ITBATlw 
•  POINTS  USED  IN  BOTH  ITERATIONS 


Fig.  IV-1.  FINITE  DIFFERENCE  PROCEDURE 
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V. 


THE  EFFECT  OF  THE  GEOMAGNETIC  FIELD 


( 


i 

* 


( 


v-l  Ion  Distribution  at  Early  Tine  -  Simple  Mouel 

The  discussion  in  Section  III-2  indicates  that  the  ion-neutral  mean 
path  may  be  appreciably  smaller  than  the  neutral-neutral  value..  In  this 
section  we  develop  a  very  simple  model  of  the  ion  motion  that  is  useful 
when  these  two  mean  free  paths  are  widely  different.  We  consider  here 
only  the  early  time  motion  of  a  cylindrically  symmetric  neutral  flow  in 
a  longitudinal  or  transverse  magnetic  field. 

In  Section  IV-3g  the  usefulness  of  the  collisionless  approximation  of 
the  neutral  motion  was  discussed.  There  it  was  indicated  that  this  approxi¬ 
mation  was  a  useful,  although  crude,  representation  of  a  one-fluid  model 
of  the  neutral  wind  during  the  initJal  stages  o''  the  expansion.  The  time 
histories  of  the  density  and  velocity  distributions  for  a  typical  high 
altitude  expansion  are  shown  in  Figs.  IV- 13  and  l4 . 

Significant  variations  in  the  neutral  flow  can  only  occur  over 
distances  greater  than  a  neutral-neutral  rean  free  path.  When  the  ion-neutral 
mean  free  path  is  much  smaller  than  the  neutral-neutral  value,  we  may  assume 
that  the  ion  motion  equilibrates  with  the  local  velocity  of  the  neutral 
wind.  To  a  first  approximation,  we  expect  that  the  ion  velocity  component 
parallel  to  the  field  lines  thus  will  remain  close  to  that  of  the  neutral 


wind 


(1) 


The  perpendicular  component  will  reach  a  quasi-steady  diffusion  value.  As 


discussed  in  Section  III-l  we  would  normally  expect  the  velocity  across  the 
field  lines  to  be  of  the  order 


V-l 


BLANK  PAGE 


for  a  finite  length  flow  field. 

Here  we  assume  that  the  ion  pressure  gradients  are  not  important  and  * 

that  polarization  electric  fields  do  not  give  rise  to  dominating  effects 
(such  as  anomalous  diffusion).  Using  the  collisionless  approximation  for 
the  flow  field  (Figs.  IV- 13  and  14)  to  specify  and  v  as  a  function  of 

position  and  time  we  have  evaluated  the  ion  density  distributions  corresponding 

* 

to  the  ion  velocity  field  given  by  Eqs .  1  and  2.  Here  both  a  longitudinal 
and  a  transverse  magnetic  field  have  been  considered.  The  density  distributions 
at  various  times  are  shown  in  Figs.  III-4,  5,  and  6.  Ion  densities  along  a 
line  parallel  to  B  and  along  a  line  perpendicular  to  B  are  shown.  In  the 
calculation  v  was  set  equal  to  6  x  10  sec-1  and  <juc  =  240  sec-1  where 
n  is  the  number  density  of  the  neutral  wind. 

The  calculation  was  initiated  at  a  time  t  =  0  assuming  an  initial 
uniform  ion  density  distribution.  A  schematic  representation  of  the  resulting 

I 

three  dimensional  ion  density  distribution  is  shown  in  Fig.  III-3. 


* 

In  a  previous  progress  report  the  calculation  was  carried  out  for  a  cylindrically 
symmetric  neutral  flow,  in  which  case  the  radial  electric  field  due  to  charge 
separation  resulted  in  a  much  lower  perpendicular  diffusion  velocity  (corres¬ 
ponding  to  the  value  given  in  Eq.  Ill -2  instead  of  l).  For  the  flow  fields  of 
practical  interest,  Eq.  V-2  is  considered  more  suitable. 


V-  2 


V- 2.  Linearized  Approximate  ->n  for  General  C  emetry 

In  the  linearized  approximation  the  equations  o?  motion  may  be  written 


in  the  form^jsee  Section  V-L)  ^ 

ri2S  <  ( e  ♦  *  O,-  -  V;  V{  (!) 


+  5aJL? )  + 

\  c  J 


(2) 


Here  a^  and  ag  are  the  accelerations  of  the  charged  particles  induced  by 
the  neutral  "wind."  In  terms  of  the  momentum  velocity  u,  defined  by 


£  »  (/VH;  4.  ^  **'*) 


(3) 


and  the  current  velocity  v,  defined  by 

v  »  OJ;.  -  ^ 


(4) 


these  equations  may  be  rec  st  in  the  form 

<5  *  V  **  «>c 


+  a,  -  v,4  + 


(5) 


and 

*  |»(4  -  (t-*)**)*-*,.  +  v«.  *-*.*).  (6) 


The  linearized  form  of  Poisson’s  equation  may  be  written  in  the  form 

v)  =  O  .  (7) 

Here  we  have  used  the  notation 

a  V  y  <8> 

a.  4T  *»€  «  A**;  ♦  '*'«) 

p1,  s 

<X  *  2*** 

V-3 


.  JM; 


^  *  (*\*Ml  1’**«>4)/(‘*C*'*0  (8) 

Vu  e  ^<^4*  *4.  v/v*a*O/0*Hi  ■+'**4.) 

Hv*  Ve-v>*  . 

In  the  general  case  the  driving  forces v  represented  by  the  accelerations 
and  a^,  include  the  ion  pressure  gradients,  and  the  ion- ion  and  ion-electron 
interactions  in  addition  to  the  neutral-ion  interactions . 

At  high  altitudes  where  the  ion-ion  collision  frequencies  are  comparable 
to  the  ion-neutral  values  and/or  at  late  times  when  the  neutral  drag  forces 
are  small,  these  contributions  to  the  acceleration  terms  can  be 
significant  and  should  be  included.  At  the  lower  altitudes  and/or  at  early 
times,  the  neutral  interactions  are  expected  to  dominate  (see  Section  III). 

For  the  present  we  will  neglect  these  ion-ion  contributions  aud  assume  that 
the  terms  and  a.^  are  known  functions  of  position  and  time . 

Even  with  these  rather  considerable  simplifications,  it  is  difficult  to 
obtain  detailed  information  concerning  the  reaction  of  the  ion  fluid  to  the 
neutral  wind  except  in  rather  simple  geometries.  Thus  we  expect  that  a 
mathematical  model  suitable  for  computation  of  an  actual  flow  field  may  be 
a  drastic  simplification  of  the  general  case.  In  order  to  gain  an  understanding 
of  the  importance  of  which  features  of  the  neutral  flow  and  neutral-ion  interaction 
which  must  be  included  in  such  a  model  we  consider  various  limiting  cases. 

In  Section  the  dispersion  relations  for  the  plasma  are  given  and 
1  he  various  characteristic  frequencies  and  propagation  modes  discussed.  In 
Section  the  nature  of  the  high  and  moderate  frequency  response  of  the 
plisma  to  the  neutral  wind  is  discussed  in  the  linearized  approximation  for 
some  simple  geometries.  In  Section  V-4  we  consider-  the  solution  of  the  non¬ 
linear  equations  of  motion,  again  for  some  simple  geometries  only. 


In  addition  to  the  equations  of  motion,  Eqs .  5,  6  and  7,  we  need 
the  two  Maxwell  equations 

<Uj*JL  i  (9) 

and 

cm*4.  C?t  *  (  £  ^  vo*  J /dv  (10) 

The  frequencies  we  ate  concerned  with  are  well  below  the  plasma  frequency 
U)p»  In  this  case  it  Is  usually  valid  to  neglect  the  displacement  current 
term  £  in  Eq.10.  lhis  is  of  considerable  advantage  since  it  materially 
reduces  the  complexity  of  both  the  analysis  and  the  numerical  evaluations. 
However,  when  polarization  charges  are  important  it  is  then  necessary  to 
introduce  Poisson's  eauation  (the  divergence  of  Eq.  10  )  explicitly: 

dU  (  £  +  £)  *  O  . 

Thus  Eq.  10  becomes 

0*4  ^  v/c1*  .  (1n 


The  approximation  of  the  magnetic  field  resulting  from  Eq.  11  is  one  which 
is  determined  solely  by  a  divergence -free  current  and  thus  cannot  be  used 
to  deduce  the  charge  distribution.  For  this  we  need  the  additional  equation 
7*  Equations  5,  6,  7,  9  and  11  form  a  complete  set  from  which 
solutions  for  the  electromagnetic  field  (£,  the  electric  current 

velocity  (v)t  and  the  mass  velocity  (u)  may  be  obtained  in  terms  of  the 
driving  forces  and  .  In  the  general  case  a^  and  depend  on  the  ion 
pressure  gradients  in  addition  to  the  neutral-ion  interactions.  When  these 
are  important  we  nee  ...ntroduce  the  continuity  equation  for  the  ion  fluid. 
This  may  be  written  in  the  (linearized)  form 

p  +  p  •  O  (12) 


where  p  is  the  mass  density  of  the  ionized  species.  Also  required  are 

energy  equations  for  the  ion  and  electron  temperatures  T.  and  T  . 

i  e 

For  the  low  frequencie-  of  our  principal  interest  (mi  v,  ,v  )  it 

i  e  ‘ 

will  usually  be  sufficiently  accurate  to  assume  that  both  the  pressures 
and  the  temperatures  are  isotropic  (Pj.  =  Jj  ,  Tx  =  ^  ),  On.  the  other 
hand,  because  of  the  large  ion-electron  mass  ratio, it  is  generally 
necessary  to  allow  for  differing  values  for  the  ion  and  electron 
temperatures . 

For  the  present  we  will  consider  only  the  approximation  in  which  the 

driving  frequencies  are  sufficiently  low  and  the  neutral  drag  terms 

sufficiently  strong  that  the  pressure  gradient,  terms  r .*?  nr  ^x\  »i-.ant, 

i.e.,  that  the  drying  terms  and  iv.  be  considered  known  a  priori  - 

This  is  formally  equivalent  to  a  cold  plasma  approximation  (if  =  =  0). 

As  discussed  in  Section  III-2  two  limiting  mode'1  a  are  of  primcipal 

interest;  the  identical  particle  model  (ifW)  (v10  =  o,  v.  =  v„.  ™  -  m  , 

ii.  1  c.  i  e' 

a^  =  0,  a  =l)  and  the  zero  electron  mass  model  (ZEM)  (r^/ra^  =  0,  or  =  0). 
For  the  identical  particle  :.odel  the  governing  equations  become 
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OJLt-£L  «. 
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The  elementary  single  species  magnetohydrodynamic  equations  result 


from  this  model  when  it  is  assumed  that  the  current  and  the  electric  field, 
reach  quasi-steady  state  values  (v  "w  0,  £  ~  0).  This  state  is  approached 
vhen  the  driving  frequencies  are  small  compared  to  both  the  plasma  frequency 
<Dp  and  the  cyclotron  frequency  u)c  -  In  this  approximation  the  electric  field 
is  given  by 

£ .  a  v  -  (17) 


and  the  equations  of  motion  may  be  reduced  to  the  form 

i5c)x  Hr  o',  -  V 

f 

and 


(18) 


uJc  e  cu-Jl  V>  T71 1  trJ  .  (!9) 

'  *  4  ^  f  c 

Historically,  these  equations  (or  the  more  general  higher  frequency  ones 
using  Eq.14  )  have  proven  quite  useful  for  describing  many  aspects  of  the 
motion  of  a  finitely  conducting  plasma  in  a  magnetic  field.  Hie  forms  given 
in  Eq.  18  and  19  contain  the  description  of  damped  Alfven  wave  motion  and 
diffusion  of  the  magnetic  field  through  the  ionized  fluid.  The  symmetry  in 
the  motion  of  the  positive  and  negative  ions  in  the  identical  particle 
model  may  be  expected  to  allow  the  formulation  and  solution  of  problems  for 
which  it  is  difficult  to  reduce  the  non-identical  particle  equations  to 
a  tractable  form.  This  applies  also  to  the  equivalent  non-linear  relations. 

Of  course,  this  model  can  only  be  expected  to  describe  some  of  the 
aspects  of  the  motion  of  a  real  plasma  in  which  the  positive  and  negative 


ion  are  different. 


A  second  model  which  exhibits  many  of  the  features  of  the  actual 


motion  which  are  lost  in  the  identical  particle  approximation  is  that  obtained 
when  the  electron  mass  is  allowed  to  approach  zero.  In  this  case  (ZEM)  the 
current  equation  (Sq  .14  )  reduces  to  the  form 

£  s  —  (  U  “  V*)  X  oSc  •  (2°) 


The  remaining  equations  are  unchanged.  this  model  the  conductivity 
parallel  to  the  field  lines  is  infinite  (i.e.,  =  0).  The  plasma  motion 

is,  in  general,  more  complicated  in  this  model  than  in  the  identical  particle 
model  because  of  the  lack  of  symmetry.  The  electrons,  now  having  negligible 
inertia,  respond  instantly  to  the  ion  motion  but  their  motion  is  typically 
very  different  from  the  positive  ion  motion.  In  the  absence  of  electric 
fields  the  electrons  move  strictly  along  the  magnetic  field  lines . 

If  it  is  known  that  polarization  charges  are  not  important  therefor 
frequencies  small  compared  to  the  plasma  frequency,  the  displacement  current 
in  Eq.15  may  be  neglected  and  the  linearized  equations  recast  in  the  form 


ft  *  ~v  wt)  X  Wc  +  a,  -  V, 


(21) 


and 


5$  *  ( &  *Iot)  —  ^  (wc.  7)  cu^St  .  (22) 


Some  aspects  of  the  difference  between  these  two  models  are  readily  seen  in 
the  corresponding  dispersion  relations.  These  are  tabulated  in  Table  V-l 
for  propagation  parallel  to  the  field  lines  and  for  propagation  perpendicular 
to  them.  If  the  cyclotron  frequency  is  much  larger  than  both  the  collision 
frequency  and  the  characteristic  fiequency  of  the  driving  forces,  the 
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dispersion  relations  for  the  two  models  arc  identical.  In  this  case  the 
magnetic  forces  are  sufficiently  strong  that  both  the  positive  and  the 
negative  ions  are  constrained  to  motion  along  the  field  lines  and  the 
difference  between  the  two  acd_-  „s  vanishes . 

Different  is  between  the  two  models  may  also  be  seen  jn  the  steady 
state  approximation  by  expiring  the  dependence  of  the  j'nass  velocity  on  the 
electric  field.  The  mass  velocity  u  may  be  expressed  in  terms  of  the 
electric  field  £  ana  the  driving  acceleration  te*m  a  .  For  the  identical 
particle  model  Eqs.l8  and  19  may  be  combined  to  yield  (when  u  =  0): 


(*x) 


STH 


=  ( ^  ^  a.\  _ 

\  *  *  ~,J  vi  + 


Tnc  zero  electron  mass  model  yields  the  expression 


ou 


ZKH 


(til'K  E  + 


au.  *** 


Bie  factors  of  2  appearing  in  Eq.  23  are  due  to  the  normalization  of  the 
electric  and  magnetic  fields  given  in  Eq.  8.  2  £  and  2  u>c  correspond  to 

the  same  electric  and  magnetic  fields  in  the  identical  particle  model  as 
do  f  and  tuc  in  the  zero  electron  mass  model#  In  Eq.  2k  two  tore  appear 
which  do  not  occur  in  the  identical  particle  model.  These  differences  have 
obvious  physical  interpretations.  First,  when  a  radial  driving  force  a  is 
applied  to  an  axia-iy  symmetric  aligned  field  configuration,  the  two  ions 
drift  in  opposite  azimuthal  directions  with  the  same  speed.  Only  when  the 


1 


* ;  - 


particles  have  the  same  mass  does  this  motion  give  rise  to  a  vanishing 
angular  momentum  of  the  plasma,  corresponding  to  the  first  acceleration 
term  in  Eq.  2h .  Second,  the  application  of  a  radial  electric  field  gives 
rise  to  opposing  radial  diffusion  velocities  for  the  two  ions,  in  this 
case,  when  the  particles  are  identical,  the  .  et.  radial  mass  flux  vanishes , 
as  indicated  by  the  absence  of  a  term  proportional  to  Ci  in  Eq.  23.  In 
the  zero  electron  mass  model  only  the  positive  ions  heve  a  non-vanishing 
radial  velocity. 


V-3  Quasi-Steady  Linearized  Approximation 

In  this  section  we  examine  the  nature  of  the  solutions  of  the  linearized 
equations  of  motion  in  a  steady  state  approximation  in  which  we  assume  the 
driving  forces  are  time -independent.  We  are  particularly  interested  in 
evaluating  „he  effects  of  the  polarization  electric  fields  induced  by  various 
types  of  neutral  winds  and  those  due  to  the  charge  separation  resulting  from 
the  mass  difference  of  the  positive  and  negative  ions.  We  will  first  consider 
the  identical  particle  model  and  then  the  non-identical  one. 

V-3a  Identical  Particle  Model 

When  the  ions  are  identical  except  for  the  sign  of  the  charge  the 
equations  of  motion  become 


u  •  %  +  a ,  - 

(1) 

V/%  *  £  4 

(2) 

v  (  £  f  v)  = 

o  . 

(3) 

We  assume  that  the  neutral-ion  interaction  is  also  independent  of  the 
sign  of  the  charge  (a^  =  0).  In  the  steady  state  approximation  these  equations 
become 

&  *  &>»*.  *  (“> 

v- 1  (v,  2  ♦  o.«“c  +  +  4^(5) 

.  o  (6) 

and 

ouA  £  *  ““  =  c  (7 J 
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Actually  it  is  not  consistent  to  set  curl  £  =  0  for  the  general  problem 
we  want  to  consider,  since  the  neutral  flows  of  interest  are  those  that 
tend  to  expand  the  plasma.  The  induced  emf  due  to  a  finite  dB/dt  gives 
rise  to  a  drift  velocity  u  which  in  turn  allows  the  field  lines  to  move 
wi+h  the  plasma.  However,  here  we  are  primarily  interested  in  the  effects 
of  the  polarization  fields.  Since  the  equations  are  linear  we  may  solve  for 
these  two  motions  separately  -  first  those  for  which  curl  £  =  0  and  second 
those  for  which  div  j?  =  0  and  then  superpose  the  results.  In  this  section 
we  are  concerned  only  with  tie  former  case.  Thus  we  may  derive  the 
electric  field  from  a  potential  cp: 

*  I  .  (8) 

Throughout  the  present  analysis  we  assume  that  the  magnetic  field 
is  c  nstant  and  aligned  in  the  z  direction.  Tire  component  of  the  current 
velocity  v  parallel  to  the  magnetic  field  may  be  derived  from  Eq.  5 

^  /(iv«)  • 

Using  Eqs.  5  and  8  ,  Eq.  6  may,  afcer  some  manipulation,  be  reduced 
to  the  fonn 

V*<P  +  +  4 

where  K  is  defined  by 


Since  the  neutral  drag  term  has  only  r  and  z  components,  the  right- 
hand  side  of  Eq.  9  may  be  written  in  the  form 


V-13 


- U  A  (  *Q>  ^c/y  \ 

Kt*  ^0  \  ±  'y'J 

where  a.  is  the  radial  component  of  the  drag  term.  Thus,  when  the 

JU 

neutral  flow  and  the  plasma  properties  are  axiaJly  symmetric  the  polari¬ 
zation  electric  fields  vanish.  In  this  case  the  analysis  described  in 
Section  v-6  may  be  applied  to  the  motion. 

An  azimuthal  variation  in  either  the  ion  density  or  the  neutral 
drag  may  result  in  a  non-zero  charge  density.  The  resulting  electric 
fields  may  considerably  alter  the  ion  motion  from  that  expected  in  the 
absence  of  polarization  fields . 

Substitution  of  Eq.  ?  into  4  yields  the  following  expression  for 
the  momentum  velocity  u: 


it  *  /  4-  1 

*  •  u; +  ^  i 

•1  i 

*4  </ 

(10) 

As  expected,  the  motion  parallel  to  the  magnetic  field  is  unaffected.  The 
motion  perpendicular  to  the  field  lines  is  the  sum  of  the  diffusion  term  uD 
in  the  direction  of  the  driving  force  and  a  term  due  to  the  polarization 
fields  Ug •  Jn  the  general  case  there  will  also  be  a  term  due  to  the  induced 
field  which  we  are  not  considering  here.  Thus  the  velocity  in  Eq.  10  is 
really  the  velocity  relative  to  the  field  lines.  The  diffusion  term  is 
the  same  as  that  given  in  SectionIII-1.  (Note  that  the  Larmor  frequency  in 
Eq.io  is  defined  in  terms  of  (n^  +  me)  whereas  in  Section  III-l  it  is  defined 
in  terms  of  m^  only). 

In  this  section  we  want  to  examine  the  polarization  term  Ug.  First 
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we  note  that 


dUo  c  *  (y«  t)  •  4.  /  (*-'  *  W«V)  •  <ii> 


Thu*?  the  notion  induced  by  the  polarization  fields  (V  x  ^  =  o)  is 
incompressible  (div  u  =  0). 

We  now  want  to  examine  the  nature  of  solutions  of  Eq.  9  in  the 
special  case  when  the  azimuthal  variations  of  both  the  electron  density 
(i.e.,  cOp)  and  the  perpendicular  acceleration  are  small.  Before  pro¬ 
ceeding  to  this,  however,  it  will  prove  convenient  to  express  the  results 
of  this  section  in  terms  of  the  current  density  J  and  the  parallel  and 
perpendicular  conductivities,  and  0^  .  The  current  parallel  to  the 
magnetic  field  is  given  by 

j  -  •  (12) 

Thus,  since  the  conductivity  is  defined  by  the  relation  a  =  ej/nC>  the 
conductivity  parallel  to  the  field  is  given  by 


From  Eq.  5 


3U  .  A  e*  cv,  . 

the  current  perpendicular  to  the  field  lines  is 


(13) 


X.  •  (£4.  *  ^5*)/ ♦  u'/4)  • 


Thus ,  the  perpendicular  conductivity  is  given  by 


a  4  eV  (  I  4  . 


(15) 


(14) 
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The  condition  that  div  ^  =  0  in  the  steady  state  (Eq.  u  )  may  he 
now  expressed  in  the  form 

<u(  t  *  %  «*)  *  -  <**(  r*  s,,,a‘A)  •  (i6) 

According  to  this  relation  the  force  exerted  by  the  neutrals  on  the 
ions  in  the  direction  perpendicular  to  the  magnetic  field  may  he  thought 
of  as  inducing  a  current  across  the  field: 

j  /  rA  /V,  •  (17) 

Since  the  total  current  must  he  divergence-free  these  wind-induced 
currents  cause  polarization  electric  fields  to  he  set  up  in  order  to 
create  additional  currents  ?  (  =  +  CJ.  €|  )  such  that 

div  (' Jq  +  J  )  =  0  (18) 

The  solution  of  Eq.  9  in  the  general  case  is  difficult  because  of 
the  anisotropy  of  the  conductivity.  However,  when  both  the  magnetic  field 
and  the  collision  frequency  are  constant,  the  anisotropic  Poisson's 
equation  resulting  from  Eq.  9  may  he  reduced  to  the  more  familiar  iso¬ 
tropic  form  by  a  simple  transformation  of  coordinates.  To  accomplish  this 
we  first  introduce  the  potential  cp  CC  =  7?)  with  the  z  51x18  aligned 
parallel  to  B. 

By  transforming  from  the  real  space  coordinates  (r,  0,  z)  to  a  set 
(r,  0,  z')  in  which  the  coordinate  parallel  to  the  magnetic  field  is 
shrunk  by  the  factor  vbj,  /o^  >  Eq.  9  may  be  reduced  to  the  isotropic 

form 

Vf)  *  -V  •(**  .  (19) 
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Here  V  is  the  gradient  operator  in  the  r,  9,  z'  coordinates  and 

i'*  *  i/fii  1 


In  order  to  gain  a  tetter  understanding  of  the  nature  of  the  solu¬ 
tions  to  Eq.19  we  now  limit  ourselves  to  the  case  in  which  the  frac¬ 
tional  azimuthal  variations  in  the  ion  density  (i.e.,  in  )  and  in  the 
perpendicular  acceleration  (a^  )  are  both  small.  Then,  retaining  only 
the  first  order  terms  in  these  quantities,  we  may  reduce  Eq.  19  to  tte 


v/Vf  * 


«( v)  - 1-  •>■*> 


where  ne  is  the  mean  value  of  the  ion  density  (averaged  over  the  azimu¬ 
thal  angle): 

is-  I  *•«*>•)  «*e 

Since  solutions  to  the  inhomogeneous  Laplace  equation  are  familiar 
and  relatively  easy  to  visualize,  it  is  useful  to  interpret  Eq,  21  as 
Poisson's  equation  for  the  (psuedo-)  electric  field  e'  ^  e'  =  v'cp)  due 

to  a  (psuedo-)  charge  distribution  p',  defined  by 


p's  _  <“<•  0  f  'W4a.i,  \ 

~  <nrv,v>  se V  M%aT^J 


The  psuedo-field  E  is  related  to  the  real  electric:  field  E  by 


E  =  E 
r  r 


K  —  V 

0  0 


(23) 


and 


E 


—  E 
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This  psuedo- charge  distribution  ,  p',  may  be  thought  of  as  a  superposition 
of  a  number  of  charged  rings,  on  each  of  which  jr  and  z  are  constant  and  on 
each  of  which  the  total  charge  is  zero.  At  each  point  the  charge  density 
is  proportional  to  the  magnetic  field  strength  (toj,  to  the  drag  term 
(" a  jv  V'  due  to  the  neutral  wind  ( thi  s  drag  term  is  roughly  of  the  order 

v  i^  iy 

of  the  perpendicular  component  of  the  velocity  of  the  neutrals  across  the 


field  lines)  and  to  the  fractional  variation  of  the  product 


.1  * 

e  1* 


Azimuthal  variations  in  n  or  a,  may  result  xrom  random  inhomo- 

e  1^ 

geneities  in  the  neutral  wind  or  in  the  ambient  ion  density  and  from 
atmospheric  stratification  of  the  neutral  or  the  ambient  iou  density. 

Quite  large  variations  in  these  quantities  may  be  expected  as  the  flew 
develops  if  the  axis  of  symmetry  of  the  neutral  flow  is  not  aligned  with 
the  magnetic  field.  In  general,  we  expect  non-zero  values  of  p'  whenever 
the  neutral  wind  velocity  across  the  field  lines  is  significantly  different 
from  zero. 

In  order  to  obtain  a  feeling  for  the  resultant  motion  of  the  ion  fluid 
let  us  consider  first  tha  special  case  in  which  the  neutral  flow  field  is 
axially  symmetric  with  the  axis  of  symmetry  parallel  to  the  magnetic  field 
/d9  =  (T)  and  in  which  the  ratios  of  the  characteristic  longitudinal 


0  a 


V-18 


dimensions  (parallel  to  the  field)  to  the  lateral  dimensions  for  both 
the  neutral  flow  and  for  the  ion  density  inhomogeneities  are  much 
larger  than  •  At  hi6h  altitudes  in  the  ionosphere  this  require¬ 

ment  is  expected  to  be  difficult  to  meet  in  practical  situations  because 
of  the  typically  large  ratios  that  usually  obtain  for  o^  /o^  (see 
Table  V-2).  However,  the  simplicity  of  this  case  warrants  its  treat¬ 
ment  first. 

Since,  in  this  case,  the  characteristic  longitudinal  dimension  is 
large  compared  to  the  lateral  dimensions,  even  in  the  r  ,  3  ,  z  coordinate 
system,  the  resultant  potential,  distribution  in  this  coordinate  system, 
as  well  as  in  the  real  space  coordinate  system,  is  approximately  two- 
dimensional  (r,  0  only). 

Let  us  ask  what  happens  to  randomly  distributed  inhcanogeneities  in 
the  ion  density.  The  total  electric  field  at  any  point  may  be  thought 
of  as  a  superposition  of  contributions  from  individual  irregularities . 
Suppose  that  the  ion  density  distribution  may  be  represented  as  a  sum  of 
contributions  of  the  form,  for  example 

/) f\c- ~  (24) 


*The  present  discussion  has  been  based  on  the  assumption  of  equal  masses 
for  the  positive  and  negative  ions.  However,  it  is  shown  later  that  the 
appropriate  parameter  for  non-equal  masses  is  still  the  ratio  of  the  two 
conductivities . 


V-19 


V-la.  Radial  Neutral  wind  interacting  with  plasma  inhomogeneities. 


Fig.  V-lh  Resulting  diffusion  of  the  inhomogenieties . 


and  that  the  characteristic  dimensions  (r  )  of  the  irregularities  are 

°i 

small  compared  to  those  of  the  neutral  flov  field.  As  a  result  of  the 
variations  in  ng  in  the  direction  of  a^  x  polarization  charge 
distributions  will  be  set  up  (see  Fig.  V-l).  The  resultant  electric 
field  near  the  center  of  an  irregularity  (r  =-  r^)  will  be  of  the  order 


•m 


X  Wj 
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(25) 


where  £"is  the  field  due  to  all  the  other  irregularities.  Since  these 
are  randomly  distributed  we  expect  0.  Equation  23  will  approximately 
represent  the  electric  field  in  the  immediate  vicinity  of  the  irregularity 


only.  Hie  drift  velocity  induced  by  this  field  will  he  of  the  order 


u  s  £***  as  5k 


■Me.  ^  **€ 
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Thus,  regions  in  which  the  ion  density  is  greater  than  the  average  tend  to 
move  in  the  direction  of  the  neutral  wind  and  those  having  a  lower  density 
move  against  the  wind  (:  te  that  the  velocity  is  not  sensitive  to  the 
dimensions  of  the  inhomogeneity,  nor  to  the  value  of  the  magnetic  field). 
Since  this  motion  is  incompressible  (Eq.  11 )9  the  motion  is  simply  a 
rearrangement  of  the  plasma,  the  regions  of  higher  density  being  blown 
outward  in  the  direction  of  the  neutral  wind.  A  mean  radial  velocity. 
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ur,  may  be  evaluated  from 


u. 


(27) 


Thus  the  mean  outward  velocity  is  given  by 


The  fact  that  the  drift  velocity  is  independent  cf  the  magnetic  field 
cannot  be  trie  in  the  U  lit  of  B  -  0.  Hie  restrictions  on  Eq.  26  are  that 
the  characteristic  dimensions  of  the  inhomogeneities  be  small  compared 
to  that  of  the  neutral  flow  field  but  large  compared  to  c.  Larmor  radius. 
Thus  when  the  Larmor  radius  is  comparable  to  the  flow  field  dimensions 
Eq.  28fail8.  When  the  characteristic  dimension  of  the  inhomogeneity  is 
comparable  to  that  of  the  neutral  flow  field  the  mean  outward  velocity 
cannot  be  obtained  from  Eq.  28  .  Tne  general  nature  of  the  motion  in 
this  case,  however,  should  he  similar;  the  high  density  regions  moving 
with  the  virid,  the  low  dei.sity  regions  against  it. 

Equations  10  and  28  thus  show  that  there  are  two  modes  by  which  the 
plasma  may  diffuse  across  the  field  lines.  In  the  first,  the  diffusion 
velocity  is  given  by 


(£9) 
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This  is  norma]  diffusion  in  which  the  diffusion  velocity  decreases  with 

_p 

increasing  magnetic  field  as  B  for  strong  fields.  In  the  second,  the 
diffusion  velocity  is 


This  is  the  so-called  anomalous  diffusion  and  is  due  to  varying  polariza¬ 
tion  electric  fields  set  up  as  a  result  of  inhomogeneities  in  the  plasma 
density.  In  this  case  of  cylindrical  symmetry  it  is  independent  of  B  so 
l^-ng  as  the  Larmor  radius  is  sufficiently  small. 

We  will  now  examine  the  case  wMch  we  expect  to  he  more  realistic,  in 
which  the  longitudinal  dimensions  of  the  neutral  flow  iield  and  of  the 
inhomogensities  are  small  compared  to  v/cr^  /oj.  times  their  lateral  dimen¬ 
sions.  In  the  r,  <d',  z  space  the  regions  of  non-zero  charge  density  are 
thus  disc-shaped,  the  thickness  in  the  direction  of  the  field  being  much 
less  than  the  radial  dimension.  The  magnitude  of  the  charge  density  is 
the  same  as  .efore  hut  now  the  volume  is  much  reduced.  Thus  we  expect 
weaker  electric  fields  and  reduced  diffusion  velocities.  In  order  to  esti¬ 
mate  the  magnitude  of  the  reduction  vc  compare  the  value  of  the  radial 
electric  field  on  the  mid- plane  of  disc  of  thickness  i  and  radius  R  to  ■‘•■he 
value  expected  for  infinite  cylinder  of  the  same  radius.  Rough  calculations 


indicate  that  a  reduction  factor  of  the  order  of 


/l  +  inSS') 
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is  expected  uhen  R  »  l .  If  the  actual  longitudinal  dimension  is  L,  then 


we  expect  the  radial  electric  field  to  be  reduced  by  a  factor  of  the  orde** 


/?r 

Fig.  V-2a  Charge  distributions  in  the  real  and  transformed  spa 
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In  the  r,  9*,  z  coordinates  this  radial  field  will  be  relatively  uniform 

(within  a  factor  of  2  to  4)  for  lz'|  £  R.  In  real  space,  therefore,  the 

/ff? 

radial  field  is  more  or  less  uniform  for  |z|  £  J  R.  In  other  words, 
the  high  conuuctivity  along  the  f-'jld  lines  causes  the  charge  distribution 


to  spread  out  more  or  less  uniformly  along  the  field  lines  over  a  width 


We  now  consider  what  happens  to  a  field-aligned  inhomogeneity  in  the 
plasma  whose  length  is  L  and  v'.dth  is  R,  Here  we  assume  tliat  the  neutral 
flow  also  has  at  least  the  length  L.  The  drift  velocity  Induced  by  the 
electric  field  will  be  of  the  order 

%  ( 0*  (33) 

where  g  -  =•,/  —  .  This  drift  motion  will  be  present  not  only  in  the 

,  fa? 

inhomogeneity  (|2|  s  1/2)  but  also  beyond  it  (out  to  2  =**  J  R).  Within 
the  inhomogeneity,  where  the  drift  velocity  is  correlated  with  the  ion 
density,  this  motion  results  in  a  mean  diffusion  velocity  across  the  field 


lines  given  by 


C,  *  f  (  14 


Beyond  the  inhomogeneity  <  z  <  y?»>  the  electric  fields  cause  the 
same  type  of  Incompressible  rearrangement  of  the  plasma  but,  because  the 
drift  velocity  is  uncorrelated  with  the  ion  density,  no  mean  transport  of 
ions  occurs. 

In  this  case,  where  the  longitudinal  extent  of  the  inhomogeneities  is 
not  too  large,  the  anomalous  diffusion  velocity  varies  as  B-1  for  strong 
fields  and  is  dependent  on  both  the  degree  of  inhomogeneity  in  the  ion 
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density  and  the  ratio  of  length  to  width  of  the  individual  irregularities . 

In  general  we  would  expect  the  irregularities  in  the  ambient  ionosphere 
to  be  field-aligned  but  we  do  not  presently  have  estimates  of  their  strength 
or  shape. 

In  Table  V-2,  we  have  listed  estimates  for  the  magnetic  field  dependent 
factors  for  the  normal  diffusion  and  the  anomalous  diffusion  velocities  as 
a  function  of  altitude.  For  the  normal  diffusion  term  the  factor  listed  is 


the  ratio  of  the  diffusion  velocity  to  the  velocity  that  would  obtain  in  the 
absence  of  the  magnetic  field.  The  same  quantity  is  listed  for  the  anomalous 
diffusion  velocity  in  the  special  cases  when  the  parameter  y  has  the  values 


1  or  10.  Here  y  is  defined  as 

c?7  tV-i*) 


Values  for  this  quantity  between  1  and  10  might  be  expected  to  be  induced  by 
a  neutral  flow  field  whose  longitudinal  dimensions  are  of  the  order  of  10 
times  its  lateral  dimensions  (a  typical  configuration)  and  whose  axis  of 


symmetry  is  not 


bh  the  magnetic  field.  In  this  case  we  would 


A  n 


e  V 

expect  values  for  f - )  of  the  order  of  unity  to  arise  as  a  result  of 

De 

the  predominantly  one -dimensional  motion  of  the  plasma  along  the  field  lines. 
Although  the  preceeding  analysis  is  not  valid  for  such  strong  variations  in 
the  plasma  density  we  assume  th^  effects  will  be  qualitatively  similar. 

Also,  when  the  neutral  flow  axis  and  magnetic  field  do  not  coincide,  polari¬ 
zation  fields  will  develop  due  to  the  second  term  in  Eq.  21 :  this  has  been 
ignored  here. 

Subject  to  these  approximations  we  would  expect  values  of  y  of  the 
order  of  unity  when  the  neutral  flow  axis  is  perpendicular  to  the  magnetic 
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field  and  of  the  order  of  10  when  the  flow  axis  is  close  to  but  not 
precisely  aligned  with  the  field.  When  the  angne  between  the  flow  axis 
and  the  field  is  less  than  the  width  to  length  ratio  of  the  neutral  flow 
it  is  rot  clear  what  order  of  magnitude  might  be  expected  for  ^A  ng/ne  ^ 

Table  1-2. Steady  state  ratios  of  diffusion  velocity  to  the  values  expected 


for  zero  magnetic  field 


Anomalous  diffusion 

Normal  Diffusion 

r,  *  r  “'ci  WT1 

Longitudinal 

Spreading 

Ratio 

s/ofl/qk 

Vb^/c^ 

L  V  u± 

J  J 

Altitude 

y  =  i 

y  -  10 

Ambient 

Shock 

Ambient 

tan 

100 

0.70 

1.0 

1.0 

1.4 

120 

.06 

0.64 

0.98 

17 

150 

.008 

.08 

0.01 

0.5 

125 

200 

(.0015) 

(.015) 

0.0004 

0.04 

(670) 

250 

(.0006) 

(.006) 

0.00006 

0.006 

(1700) 

300 

(.0002) 

(.012) 

0.000006 

0.0006 

(5000) 

-  the  ratio  of 


In  Table  V-2  we  also  list  che  spreading  ratio,  -/o^  /o^ , 
the  distance  along  the  field  lines  that  the  polarization  charges  are  spread 
to  the  lateral  dimensions  of  the  flow  in  the  steady  state.  Above  200  km 
altitude  the  values  are  of  the  order  of  500  o.'  greater.  These  very  high 
values  are  not  expected  to  be  obtained  in  practice  for  two  reasons.  First, 
since  the  lateral  flow  dimensions  at  200  km  altitude  are  of  the  order  of  a 
few  kilometers,  this  ratio  corresponds  to  spreading  distances  of  the  order 
of  thousands  of  kilometers.  Over  these  distances  along  the  field  lines  the 
collision  frequencies  cannot  be  assumed  to  remain  even  approximately  constant. 
As  the  field  lines  penetrate  to  lower  altitudes  the  perpendicular  conductivity 
increases^ effectively  tending  to  neutralize  the  charge  distributions  more 
rapidly.  In  the  other  direction,  as  the  field  lines  climb  to  higher  alti¬ 
tudes,  the  cross  conductivity  decreases,  reducing  the  leakage  current.  This 
altitude  effect  predominantly  determines  the  manner  in  which  the  charges  are 
finally  neutralized  far  from  the  neutral  flow  perturbation  but  is  not  expected 
to  have  much  influence  on  the  electric  fields  in  the  Immediate  neighborhood 
of  the  neutral  disturbance. 

A  more  important  effect  is  the  result  of  the  fact  that  the  neutral  flow 
is  not  a  steady-state  phenomena  and,  under  these  high  altitude  conditions, 
the  quasi-steady  approximation  may  fail  as  far  as  the  anomalous  diffusion  is 
concerned.  The  argument  is  as  follows:  We  visualize  the  neutral  motion  as 
a  primarily  lateral  expansion  with  characteristic  expansion  velocity  V  and 
starting  at  time  t  =  0.  Due  to  various  plasma  or  neutral  flow  inhomo¬ 
geneities  or  to  positive -negative  ion  asymmetries  (see  Section  ),  charge 
separation  occurs  within  the  confines  of  the  reutral  flow  field.  Because  of 
the  high  conductivity  along  the  field  lines  there  is  a  tendency  for  these 


V -28 


charges  to  spread  out  rapidly  in  the  direction  of  the  field.  We  expect 
that  this  spreading  is  accomplished  as  a  disturbance  propagating  a'*-  the 
Alfven  speed  along  the  field  lines.  The  nature  of  this  Alfvl’n  motion 
is  discussed  in  Section  V-6.  Thus  the  maximum  extent  over  which  the  charges 
may  bn  srread  by  the  time  t  is  C^t.  Since,  at  this  time,  the  lateral  flow 
dimensions  are  of  the  order  of  Vt  we  expect  that  the  maximum  value  of  the 
spreading  ratio  is  CA/v.  During  the  initial  expansion  of  the  neutral  flow 
we  expect  values  of  >'  between  1  and  h  km/sec  (say  a  mean  value  of  2  km/sec). 
Thus,  maximum  spreading  ratios  during  the  firsL  several  seconds  of  the 
motion  are  expected  to  be  of  the  order  of  100.  (Here  we  are  referring  only 
to  the  spreading  of  disturbances  whose  lateral  dimensions  are  comparable 
to  those  of  the  neutral  flow  field.  Smaller  disturbances  will  attain  their 
quasi-steady  state  values  more  quickly.) 

According  to  the  above  argument  we  expect  that  the  appropriate  values 
for  the  anomalous  diffusion  velocities  above  an  altitude  of  150  km  are  not 
as  low  as  those  shown  in  Table  V-2  but  may  remain  more  comparable  to  the 
values  at  150  km.  Thus  it  appears  that  below  aoout  200  km  the  normal  diffu¬ 
sion  mechanism  can  account  for  most  of  the  motion  across  the  field  lines  in 
the  region  of  the  neutral  flow.  Above  about  200  km  the  electric  fields  due 
to  polarization  charges  may  contribute  much  of  the  cross  field  motion.  The 
relative  importance  of  each  mechanism  depends  on  several  parameters  (shock 
strengths  and  velocities,  ambient  ion  density,  angle  between  trajectory  and 
magnetic  field,  degree  of  inhomogeneity,  etc.)  so  that  these  values  can  only 
be  expected  to  be  indicative.  More  detailed  examination  of  specific  situa¬ 
tions  will  he  required  to  get  quantitative  results. 
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The  present  considerations  indicate  that  at  some  altitude  between  120 
and  150  km  the  motion  of  the  ions  across  the  magnetic  field  begins  to  be 
appreciably  reduced  below  the  neutral  velocity.  Above  20C  km  altitude  the 
lateral  velocity  of  the  ions  may  be  expected  to  be  between  1  and  10jt  of  the 
neutral  velocity  during  the  first  few  seconds  or  tens  of  seconds  of  the 
neutral  motion. 

At  lute  times,  when  the  neutral  density  is  comparable  to  the  ambient 
density  and  the  polarization  fields  have  reached  their  quasi- steady  values, 
the  cross-field  motion  will  be  even  further  suppressed. 


V-jb.  Linearized  Steady  State  Approximation  for  Non- Identical  Ions 


When  the  positive  arid  negative  ions  do  not  have  identical  properties, 
charge  separations  can  result  from  this  asymmetry.  In  this  section  we 
develop  tuo  steady  state  equations  in  the  approximation  that  the  electron 
mass  is  small  compared  to  the  ion  mass.  For  small  values  of  me/m  ,  Eqs.  V-2-5 
and  6  may  be  recast  in  the  form  (for  u  =  v  =0) 


4  r  (  ^ <Vr)/^  (l) 

and 

Here  we  have  defined  a  source  term  according  to 


This  term  represents  a  current  source  and  is  non-zero  in  general  because 
of  the  asymmetry  between  the  interactions  of  neutrals  with  the  ions  and 
neutrals  with  the  electrons . 

Equation  2  may  be  rewritten  to  express  the  current  velocity  v  in 
terms  of  the  driving  forces : 


where  f  is  an  effective  electric  field  defined  by 


F?  I  +  'ixtZS t/»,  -  ft*. 


(5) 


Here 


<£.  v,  r 

4r(L  ^  [ 


\  +  gN  ^  /«£ _ 

(  \4  +  (V»AJ«) 


;  (6) 
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g-  .  uJi 
*•  '  <nr 


l  it  r  */ui*  - — — 1 

rC  ^  ^  |+  4  (V,/uic)  J 


(7) 


and 

\  *  *  (8) 

The  steady  state  condition  that  div  7=0  and  curl  £=0  yield  an  equation 
for  the  electric  potential  due  to  the  polarization  charges : 


-  in  G  (9) 


where  the  source  function  tr  is  defined  by 


Q  ■  ^  c  -  \)  +  -  ?»  ^  b)  -  <V  « 


*8 


(10) 


As  was  done  in  Section  V-3&  we  assume  first  that  the  ratio  of  the  parallel 
conductivity  on  to  the  perpendicular  (Pedersen)  conductivity  op  is 
independent  of  position.  We  also,  as  before,  transform  to  the  space 
r' ,  9 , ,  z  *  in  which 


r  =  r 


6=9 


,  F*' 

z  =V5" 


(ii) 


If  V'  represents  the  gradient  operator  in  the  r' ,  Q',  z'  space, we  may 
rewrite  Eq.  in  the  form 


+(V'<P  s 


«  v*<s 


(12) 
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If  we  assume  that  the  spattal  variations  in  the  Pedersen  and  Hall 
conductivities  are  small,  then  the  first  order  approximation  to  the 
potential  equation  becomes : 


V'<P  -- 


dux 


*(13) 


'Vhen  both  the  plasma  and  the  driving  neutral-ion  interaction  are  axially 
symmetric  (we  always  assume  that  the  driving  forces  have  r  and  z  components 
only),  Eq.13  reduces  to  the  form: 


-  -U  'JLuo 


( 


7 


(14) 


2  2 

Throughout  most  of  the  ionosphere  the  term  0  ^^2^0  is  KmB  1 1  co“P&red 

to  unity  (see  Table  III-l)*  In  this  case  ^  VlaP  .  and 

o-a  — — 


q:  jg  <*>?_  v, 

MTTC 


(15) 


The  ratio  of  perpendicular  to  parallel  conductivities  is  given  by 


f  >>■ 

<■*£♦  V" 


(16) 


Since  we  expect  values  for  to  be  smaller  than  those  for  a^  by  a  facte 
of  the  order  of  Jae/mi,  the  source  function  G  in  Eq.  10  is  approximately 


* 

This  assumption  is  made  throughout  the  remainder  of  this  section. 


represented  by 


For  the  axially  symmetric  case,  the  pol  ential  equation  then  tak  ,-s  the 
form: 

v'*<p  *  -  rr  ••j.  "  %)  .  n«> 

ITie  first  term  on  the  right-hand  side  of  Eq.  ,\8  results  from  the  fact 
that  the  ions  move  more  easily  across  the  field  lines  than  do  the  electrons . 
The  second  term  is  independent  of  the  magnetic  field  and  is  a  result  of  the 
difference  bet.veen  the  acceleration-  imparted  by  the  neutrals  to  the  ions 
and  to  the  e '  ;rons . 

The  drift  motion  due  to  t:»-  radial  variations  of  the  driving  term 
may  be  evaluated  as  follows.  Due  to  this  '  m,  the  electric  field  has 
value,  vi+uin  the  confines  of  the  neutral  flow  field,  of  the  order  of 

f  W *  -  -C  Q, 

cr  w*>  *  t 

where  f  is  the  reduction  factor(defined  in  Section yjT-ja  due  to  the  spreading 
of  the  charge  distribution  along  the  field  linet  * 

$  &  4  0* 


Using  Eqs .  1,2  and  we  may  i'OUghly  evaluate  the  magnitude  of  the  outward 


*  (19) 


We  may  estimate  the  outward  drift  to  be 


Thus,  when  the  reduction  of  the  polarization  field  by  anisotropic 
conduction  along  the  field  lines  is  considerable,  (f  «  l)  the  outward  drift 
motion  is  determined  by  the  ions  alone  and  has  the  same  value  as  in  the 
identical  particle  model  (Section  V-3a).  For  cylindrical  symmetry, 
however,  there  is  no  conduction  along  the  field  lines  (f=l).  In  this 
case  the  outward  drift  velocity,  although  only  evaluated  by  Eg. 
to  order  of  magnitude,  is  much  smaller  (c.f.  Eq.  III-1-2)  due  +o  the  fact 
that  the  electrons  now  have  to  be  dragged  across  the  field  lines . 

When  the  variations  of  the  plasma  density  in  space  are  appreciable, 
additional  polarization  fields  will  arise.  As  was  shown  in  Section  V-3a 
these  can  result  in  an  appreciable  increase  in  the  rate  of  diffusion  of 
the  plasmc  across  the  field  lines.  We  will  not  attempt  a  detailed 
evaluauion  of  these  motions  here  but  simply  identify  the  terms  giving 
rise  to  the  drift  motion.  This  case  is  somewhat  different  than  when  the 
ions  are  identical  since  here  a  polarization  field  exists  even  when  the 
plasma  is  uniform. 

In  the  linear  approximation  we  may  separate  the  electric  field  into 
two  terms  £0(=^90)  and  £ .  j£o  is  the  field  that  would  exist  in  a 
uniform  plasma  and  is  the  contribution  due  tc  variations  in  the  plasma 
density  and  (azimuthal  only)  in  the  neutral  flow.  Here  cp0  may  be  evaluated 
from 

V/Z(fi  *  -  +  s.  V.  X.  (21) 

•  ^  * 
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The  potential  due  to  the  variations  in  plasma  density  is  determined  by 


2 _«!**<; 


V  ‘m  v,  J  ~ 


<19 


^■55 

a;  'br 


P»,  -^5 

■5?  *y‘ 


The  first  term  in  Eq.  21  is  the  only  one  appearing  in  the  identical  ion 
model  (Eq.  9  ).  The  other  terms  result  from  the  fact  that  the  positive 
and  negative  ions  are  different.  The  first  two  terms  in  Eq.  21  give  rise 
to  electric  fields  which  cause  predominantly  radial  drift  motions.  The 
third  and  fourth  term  result  ir.  predominantly  azimuthal  drifts  and  the 
last  term  does  not  give  rise  to  ~  drift  motion  (to  a  first  approximation). 

More  detailed  considerations  similar  to  those  in  Section  v_3&  would  be 
required  to  evaluate  the  effect  of  these  additional  terms.  In  any  case 
it  is  clear  that  the  electrons  affect  the  plasma  motion  primarily  through 
the  polarization  fields,  and  that  the  importance  of  the  electron-neutral 
and  electron-ion  interactions  depends  on  the  extent  to  which  the  electrons 
can  neutralize  the  polarization  charges  by  conduction  along  the  field  lines . 
In  general,  this  question  can  only  be  answered  by  considering  the  geometry 
and  time  dependence  of  the  specific  problem  of  interes-  in  some  detail. 

In  the  following  paragraphs  we  restate  this  conclusion  in  more  general 
te  nns . 

The  first  term  on  the  right-hand  side  of  Eq.  21  is  the  (pseudo)  charge 
density  that  results  from  the  fact  that  the  electrons  are  more  strongly 
bound  to  the  field  lines  them  are  the  ions .  In  the  cylindrically  symmetric 
case,  the  charge  neutral  condition  requi.es  that  the  radial  velocities  of 
ions  and  electrons  be  approximately  equal.  Rather  strong  radial  fields  must 
then  be  set  up  to  pull  the  electrons  across  the  field  lines  with  the  ions. 


When  the  flow  dimension  in  the  dx  ’ection  parallel  to  the  field  is 
not  too  large,  the  electrons  do  not  necessarily  have  to  move  with  the 
ions.  In  this  case,  when  no  azimuthal  polarization  electric  fields  are 
present,  the  motion  across  the  field  lines  is  primarily  of  the  ions  only. 
The  electrons  move  relatively  freely  along  the  field  lines  in  an  attempt 
to  maintain  their  density  approximately  equal  to  the  ion  density. 

When  the  distance  the  electrons  have  to  move  to  compensate  for  a 
given  lateral  ion  motion  is  not  too  large  (i.e.,  small  compared  to  Ja^ /ap 
times  the  ion  displacement),  the  forces  required  to  move  the  electrons 
are  small  compared  to  those  required  for  the  ions.  Since,  in  the  absence 
of  an  externriiy  applied  electric  field,  the  driving  forces  exerted  on  the 
electrons  are  small  compared  to  those  exerted  on  the  ions,  the  motion  of 
the  plasma  is  largely  controlled  by  the  properties  of  the  heavy  ions  and 
their  interactions  with  the  neutrals  and  is  insensitive  to  the  electronic 
properties.  However,  in  certain  situations  the  electrons  can  have  a 
considerably  stronger  influence  on  the  plasma  motion  than  simply  the  mass 
ratio  would  indicate.  For  example,  when  the  symmetry  of  the  neutral  flow 
field  is  such  that  the  electric  fields  induced  by  the  ion  motion  are 
nearly  perpendicular  to  the  magnetic  field  lines,  relatively  large 
electric  fields  have  to  be  set  up  in  order  to  get  significant  motion  of 
the  electrons.  Also,  when  the  distances  the  electrons  have  to  travel  in 
order  to  compensate  for  a  given  ion  motion  are  much  xarger  than  the  ion 
displacements,  the  electrons  must  be  accelerated  to  considerably  higher 
velocities  than  do  the  ions.  These  effects  will  tend  to  increase  the 
influence  of  the  electronic  properties  on  the  motion  of  the  plasma. 


Formulation  of  the  Nonlinear  Equations  for  Simplified  Q 


retries 


In  the  region  of  the  iorosphere  under  consideration,  neutral  particles 

U 

are  about  10  times  more  numerous  than  charged  particles.  On  the  other 
tend,  the  termor  cyclotron  frequency  of  each  charged  species  is  considerably 
higher  than  its  collision  frequency.  In  this  situation,  since  the  earth's 
magnetic  field  is  rather  weah,  a  one-fluid  hydrodynamic  calculation  of  the 
motion  of  the  neutral  particles  may  be  carried  out  without  regard  to  the 
currents  or  ponderomotive  forces  which  this  motion  induces.  The  ions  and 
electrons,  however,  describe  a  number  of  Larmor  cycles  in  the  geomagnetic 
field  between  collisions,  and  their  motion  is  therefore  strongly  affected 
by  the  magnetic  field.  The  following  simple  example  provides  some  insight 
into  the  situation. 

Consider  a  charged  particle  with  charge  q  =  ±e,  immersed  in  a  uniform 
magnetic  field  and  a  bath  of  neutral  particles.  If  collisions  are  infre¬ 
quent  compared  with  the  rate  of  Larmor  gyration  in  the  magnetic  field,  we 
may  view  the  collision  as  a  displacement  of  the  guiding  center  of  the 
particle.  A  collision  will  move  the  guiding  center  by  at  most  2  rL 

(r  -  mVC/e^  is  the  termor  radius).  We  expect  the  average  "drift"  velocity 
L 


Fig.  1  Iod  in  a  Neutral-Particle  Bath 


of  the  guiding  center  to  be  of  the  order  of  (t//a>  )  V^,  where  u  is  the 
collision  frequency  with  the  neutrals.  The  collisions  will  predominantly 
o-cur  in  the  portion  of  the  cycle  when  the  ion  is  moving  against  the 
neutral  stream,  and  the  collisions  will  tend  to  make  the  particle  drift  to 
the  right  (in  Fig.  l)  if  it  is  positively  charged  and  to  the  left  if  it  is 
negative.  A  calculation  of  the  drift  motion,  taking  averages  over  all 
collisions  (Appendix  II ),  will  reveal  that  the  lateral  average  velocity 
(in  the  direction  q  X  B)  is  and  the  drift  parallel 

to  V  is  at  the  speed  v  {vc  +  co^)"  V^.  The  drift  velocity  is  always  per¬ 

pendicular  to  the  direction  from  which  the  neutrals  appear  to  come,  but  to 
the  ion  in  its  drifting  motion  the  apparent  direction  of  the  neutral  stream 
changes;  hence  the  component  parallel  to  arises.  The  above  simple  pic¬ 
ture  neglects  electric  fields. 

Since  the  sign  of  the  component  of  drift  motion  perpendicular  to  is 
dependent  on  the  sign  of  the  charge,  an  electric  current  will  arise.  This 
current  will  cause  a  change  in  B  (though  slight),  end  the  time-rate  of  this 
change  will  induce  an  electric  field,  according  to  Ampere's  law.  This 
electric  field  will  bring  about  an  additional  drift  of  the  charged  particles, 
according  to  the  well-known  expression 

Vd  =  c  E  X  B/B^  , 

which  is  independent  of  charge.  The  simple  picture  of  figure  1  will  hold  if 
we  perform  a  Lorentz  transformation  to  the  frame  of  reference  in  which  E  =  0 
if  E  does  not  vary  significantly  over  a  larmor  radius. 

In  the  following  sections,  we  will  undertake  a  detailed  solution  of  the 
ion  and  electron  motion,  including  the  influence  of  the  magnetic  and  induced 
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electric  fields  and  the  neutral  flow  field,  in  the  vicinity  of  a  missile 
traversing  the  ionosphere.  The  philosophy  will  be  to  understand  one 
simple  geometry,  namely  that  of  cylindrically  symmetric  neutral  flow  with 
the  geomagnetic  field  along  the  axis  of  symmetry.  We  will  also  discuss 
preliminary  attempts  to  apply  these  results  to  more  general  neutral  flow 
fields,  for  example  the  spherically  symmetric  motion  of  the  neutral  gas 
in  a  uniform  magnetic  field. 

In  the  analysis,  it  will  be  assumed  that  ions  and  electrons  collide 
only  with  neutrals.  This  approximation  will  be  justified  on  the  basis  of 
the  relative  number  densities  and  cross  sections.  It  will  also  be  shown 
that  the  relative  change  of  the  geomagnetic  field  induced  by  the  flow  is 
very  small.  Any  ionization  or  recombination  occurring  during  characteristic 
flow  times  will  be  neglected,  as  they  have  little  bearing  on  the  plasma 
motion  and  are  very  difficult  to  include. 

A  number  of  other  simplifications  arise  from  the  large  time  scale 
peculiar  to  the  flows  under  consideration.  For  example,  since  the  mean  free 
time  and  mean  free  path  of  the  neutrals  are  of  the  order  of  0.1  to  1  sec  and 
1  km,  the  passage  of  the  neutral  flow  past  a  point  in  s^ce  excites  plasma 
oscillations  (whose  frequency  is  of  the  order  of  10  mc/sec)  only  very  weakly . 

In  addition  to  these  effects,  it  is  expected  that  the  electron  motion 
is  coupled  to  the  ion  motion  through  the  magnetic  field,  effectively 
limiting  the  speed  at  which  divergence-free  current  disturbances  in  the 
electron  fluid  can  propagate  through  the  plasma.  This  limiting  velocity 
is  expected  to  be  of  the  order  of  the  ion  Alfven  speed. 
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Basic  Equations 


The  motion  is  governed  by  the  following  equations  (using  c.  g.  s. 
Gaussian  units)  }  with  subscripts  i  and  e  referring  to  ions  (of  various 
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The  collision  integral  for  momentum  exchange  of  the  Jth  species  with  all 
other  species  is;1 

%  s2~  zrh  k-  V  f  <$>  too  s  /  H/  ^ 

k+4  •  V 


(8) 


where  C,  is  the  velocity  of  a  single  particle  of  the  kth  species  and  f  is 
the  velocity  distribution  function  for  the  kth  species.  C  represents  the 
velocity  after  a  collision.  The  momentum  transfer  cross  section  used  here 
is  defined  as 
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(9) 


(T  —  2tT  j  rr( 9  )( 1-cos  9  )  sin  9  d9 
o 

is  the  mean  velocity  of  the  Jth  species  and|p,  is  its  pressure  tensor. 

We  will  consider  only  cases  where  there  is  a  single  ion  species,  fcr  sim¬ 
plicity.  The  energy-conservation  equation  is  omitted  since,  in  all  cases 
considered,  the  pressure  gradient  term  in  Eqs.  (6)  and  (7)  will  be  neglected 
in  comparison  with  x,he  momentum-exchange  integrals  I_^,  Ie*  In  this  approxi¬ 
mation  Eqs.  (l)  to  (7)  form  a  closed  system.  No  chemical  reactions  have 
been  included.  However,  charge  exchange  between  an  ion  and  a  parent  atom 
(same  species),  which  is  the  dominant  charge -transfer  process,  will  be 
included  as  part  of  the  momentum-exchange  integral,  since  I'or  3uch  a  pair 
the  results  of  charge-exchange  and  momentum-exchange  are  equivalent. 

Formidable  difficulties  arise  in  the  evaluation  of  the  collision  inte¬ 
gral  in  (8).  The  distribution  functions  f  ,  f  can  hardly  be  called 

J  K 

Maxwellian  with  any  confidence  since  we  must  allow  for  the  possibility  of 
qui^e  large  ion  "slip"  relative  to  the  neutrals;  a  three-fluid  model  using 
Maxwellian  distributions  could  be  used  as  a  good  approximation  in  the  event 
that  the  fluid  were  in  quasi-equilibrium.  However,  the  criterion  for  the 
validity  of  this  approach  is  not  clear  when  the  mean  relative  velocities 
are  not  small  compared  with  thermal  velocities.  It  would  be  very  inaccurate 
to  approximate  the  f  's  with  6-func+ions,  since  the  thermal  velocity  of  the 

V 

molecules  url  ions  behind  the  shock  is  of  the  same  order  as  the  velocity  of 

propagation  of  the  shock  into  the  undisturbed  surrounding. 

Even  if  the  f  ,'s  were  known,  ct  is  'isually  not  a  convenient  function  of 
0  m 

the  relative  speed  of  the  interacting  particle*. .  The  integration  of  Eq.  (8) 


See  Section  III-2. 


\ 


is  quite  complicated  even  if  we  assume  some  r  potential  law  between 

particles.  Therefore  we  will  simply  lump  all  of  the  uncertain  factors  in 
the  integral  in  (8)  into  a  collision  frequency,  i>  ,  which  we  may  assume  lo 
be  constant  to  a  reasonable  approximation: 


'-V  tv  ^  v.v  ( VXi ) 

(p  is  proportional  to  the  density  n  . ) 
jk  k 


(10) 


Equation  (10)  is  exact  for  an  interaction  potential  which  varies  as 


r  ,  (the  so-called  "Maxwellian-molecule"  force  law).  In  this  case  a 


m 


-1 


varies  as  ^ V  ) 

rel  ,  and  the  integrations  over  arbitrary  f  ,  f  become  trivial. 

C  * 

2  -4 

Dalgarno  has  shown  that  the  r  interaction  potential  is  a  rather  good 

approximation  for  ion-neutral  encounterr  at  low  energies  (~  1  eV).  For 

hard-sphere  potentials  (n  =  <*■ ,  independent  of  energy),  ZhdanovJ  has 

shown  that  for  the  quasi-equilibrium  case  [  )  «  VThe™a  J 

collision  integral  I  again  takes  the  fvro  shown  in  (10). 

Since  a  large  part  of  the  subsequent  analysis  will  involve  the  simple 

case  of  flow  which  is  cylindrically  symmetric  about  the  geomagnetic  field 

direction,  we  now  obtain  the  form  of  Eqs.  (l)  to  (7)  in  this  special  case. 

d  _  d 


In  cylindrical  coordinates,  with  §T  =  5^  ~  and  dropping 


1  \  ‘  ‘  e  V  «• 

0  s  WKu;  -«V*)  4  i  ^ 
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(12 ) 


(13 ) 
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Here  u,  v  are  radial  and  azimuthal  components  of  velocity.  The  neutrals 

have  only  a  radial  velocity^.  All  species  of  neutrals  are  assumed  to 

have  the  same  mean  velocity,  when  averaged  over  a  few  collision  times  or 

mean  free  paths.  Since  Ez  =  0,  there  being  no  motion  along  the  z-direction, 

the  r  &  9  components  of  Eq.  (l)  yield  SB ^/dz  -  0,  dB^/dt  =  0;  hence 

B  =  EL  =  0  for  all  time, 
r  y 

We  have  elected  to  work  with  the  equations  of  motion  of  the  ions  and 
electrons,  since  it  is  more  direct  than  to  think  of  such  motion  in  terms  of 
an  ''Ohm’s  Law"  which  Includes  Hall  current  and  ion  slip.  The  Ohm's  Law 
approach  results  in  a  simplification  only  if  a  steady  state  is  assumed.  It 
is  shown  subsequently  that  E-^s.  (6)  and  (j)  lead  to  the  usual  form  of  the 
generalized  Ohm's  law  in  a  steady,  uniform  state. 
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Before  proceeding  to  numerical  calculations,  we  will  learn  what  we 
can  from  a  linearized  approximation  to  the  governing  equations.  We  will 


consider  the  general  axially  symmetric  case,  since  the  assumption  of 
cylindrical  symmetry  is  needlessly  restrictive  in  the  linearized  approxi¬ 
mation  and  since  we  will  want  to  discuss  more  general  axial  symmetry  later. 

The  object  of  the  present  section  is  to  determine  the  accuracy  of 
several  levels  of  simplicity  in  the  governing  equations.  This  is  done  in 
the  hope  that  it  will  be  possible  to  neglect  te’nos  which  are  related  to 
high-frequency  oscillations  of  the  plasma,  and  thus  make  possible  numerical 

#5  -.3 

calculations  v:ith  reasonably  long  time  steps.  For  example,  if  n  -  10  cm  , 

T  -1 

the  plasma  angular  frequency  is  5*5  x  10  sec  ,  whereas  we  will  be  con¬ 
cerned  with  flow  times  of  the  order  of  one  second. 

Assuming  a  constant  density  of  neutrals,  n^,  constant  ion  and  electron 
collision  frequencies  with  neuorals,  v  and  v and  treating  vg, 

(n  -  neQ),  (n^  -  niQ)  and  (B  -  Bq)  as  first-order  quantities  and  neglecting 
V*jpwe  obtain  from  Eqs.  (l)  to  (7)  in  the  axially  symmetric  case  the 
following  equations  (with  cylindrical  coordinates  and  d/S9  =  0): 


- 

lls  * 

1 

—  Vt 

«.  -  s 

Tr 

M. 

1 

r 


(20) 


(21) 


(22) 


V- 45 


(  U  -  -U€\  4> 

4.  ^ 

(24) 

«- 

\  *■  */ 

«■  'd-r 

/  V  1  Ert 

*y» 

c  • 

C 

(^_-ut)+  i  _. 

(25) 

£(rd»)  * 

“‘)  *  k  ^ 

( 2<  ) 

®£  +-s(i 

•  ?rM- 

^.v\  «„ 

’S' 

(27) 

&Jis  +  <v\«/.i 
\ r 

4  _  n 

(28) 

1».  £.(««• 

^4*)  -  V;(^-V„) 

(29) 

_ 4( 

«3*  *1'' 

1  «- 

)-*  <?*-*-) 

(30) 

Taking  tbe  curl  of  Eqs.  (23)  to  (25)  and  d/St  of  (26)  to  U8),  we 
obtain 
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(33) 


Thus  n  is  uncoupled,  as  are  is.  and  a  ,  and  there  remain  9  scalar  equations 

1  6 

in  the  unknowns  v,  ,v  ,  and  E. 

i '  e  ‘ 

The  plasma  motion  may  be  viewed  as  a  forced  oscillation  with  a  forcing 
function  V^t)  independently  prescribed.  1  problems  of  concern  here, 
significant  variations  of  the  flow  field  require  significant  fractions  of 
one  second  to  take  place;  thus  we  expect  the  Fourier  components  of  V  (t) 
at  tne  plasma  frequency  to  be  small.  An  important  simplification  will 
therefore  be  the  assumption  that  n^  =  n&,  which,  as  we  shall  see,  elimi¬ 
nates  the  unimportant  high  frequency  oscillations.  In  this  analysis 
we  restrict  ourselves  to  the  condition  that  the  radial  and  longitudinal 
motions  of  the  positive  and  negative  ions  are  identical.  This  result 
follows  from  the  assumption  n  =  ng  when  the  particles  are  identical  in 
general  or  when  the  driving  interactions  have  cylindrical  symmetry  if 
the  particles  have  different  masses.  In  the  general  axially  symmetric 
configuration  for  differing  particles,  longitudinal  electron  currents  may 
radically  alter  the  coupling  of  the  ions  with  the  magnetic  field.  However, 
fc  -  the  present  analysis,  we  are  concerned  with  examining  the  non-linear 

effects  and  we  shall  make  the  assumption  that  u  =  u.  =  u  and  w  =  v,  =  w 

e  i  e  i 

^Eqs.  29,  30)  in  order  to  get  a  tractable  first  approach,  'f’his  permits  only 
azimuthal  currents,  which  do  not  change  the  Coulomb  E-field.  In  order  that 
this  model  be  self-consistent,  it  becomes  nezessary  to  drop  Eqs.  (33)  and 
(35),  and  to  determine  E  ,  E  by  the  requirements  that  u  =  u  and  w  =  v 
in  uhe  r,  z  components  of  Eqs.  (31)  and  (32).  We  then  have  four  coupled 
equations  (assuming  v^  =  0,  a  non-rotating  neutral  flow  field): 
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Xu  addition,  there  are  the  uncoupled  equations: 


(38) 


(3°) 


(l+  ^)\*^  *  -  (y;  +  ***v«/«>;)(w-  (*o) 

fc.\  =  T  +*(w— O)  Cl) 

E<-  '  ”  ^  +  T  ( it  +v‘(“'u^)  (',2) 


In  the  identical  particle  model,  the  electrons  and  ions  have  equal 

masses,  and  equal  but  opposite  charges.  In  this  case  (from  (38) 

and  (39) ),  =  ve  =  v,  and  one  more  dependent  variable  is  eliminated. 

Furthermore,  substituting  Eqs .  (37)  and  (40)  into  (4l)  and  (42)  this  limit 

yields  E  =  E  =  0.  We  obtain  a  similar  form  of  the  equations  when 
r  z 

me  /  m^,  by  introducing  the  variables: 

AT  = 

V  S  ^  +  '**«**/*< i  (43) 

Then,  from  (38)  and  (39)> 
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where  m  -  i/e'  =  (m^  +  sy^/n  and  v ±'  =  (my  + 

From  (45),  we  estimate  v  to  be  of  the  order  of 


**■/*"  v/ 


m  y  }/m. 
e  e  ' 


Since  we  expect  «s  (a^m^)2,  the  term  involving  v  in  (44)  is  smaller 

_x_ 

than  the  term  involving  v  by  (tn^/m.  )2.  Thus  for  nig  «  m  ,  eq  uation  (44) 
takes  the  form: 


>r  ^  (p  _  uB.\  .  o'/tr 

•R  ~  /**.>**«  —  /  *» 


(46) 


Equation  (46)  has  the  same  form  as  does  the  combination  cf  (38)  and  (39)  for 
ions  of  equal  mass  and  opposite  charge.  The  resulting  simplification  v.ll 
be  used  extensively. 

We  now  have  three  models  of  the  linearized  problem,  which  we  wish  to 
compare  in  detail..  These  are  represented  by: 

a)  Equations  (31)  to  (35)*  the  general  problem,  which  we  shall  be 
able  to  handle  only  for  cylindrical  symmetry. 

b)  The  electrical  neutrality"  model,  n^  =  n^,  which  employs 
Eqs.  (36)  to  (39),  and  welch  holds  for  arbitrary  axial  symmetry 
for  identical  particles  in  general  and  for  non-identical  particles 
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when  it  can  be  shown  that  the  and  z  motions  are  the  same. 


c)  The  'very  unequal  mass"  model  (in  «  nr  and  n  *  n  ) ,  in  which 

e  i  i  e 

Eqs.  (38)  and  (39)  are  replaced  by  (46),  and  which  has  the  same 
mathematical  form  as  the  "equal  mass"  model;  it  also  permits 
a  formal  axial  symmetry. 

We  now  wish  to  express  the  equation  for  the  electric  field  in  terms 

of  a  source  function.  To  do  this,  we  first  perform  a  Fourier-transformation 
in  time  oi  the  equations  governing  each  model.  In  each  of  the  three  models, 
we  substitute  into  Eq.  (34)  from  the  other  pertinent  equations  to  obtain  a 
differential  equation  in  space  for  the  transform  which  has  the  form 

V  €  +  t  (  k  (  *o  *  j  •  ft  C  ***) 

where  ^ 

The  functions  A(w)  and  k*"(to)  are  shown  below  for  each  model. 

(a)  General  Case  (cylindrical  symmetry). 
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In  the  numerical  analysis,  we  will  be  restricted  to  the  cylindrically 
symmetric  case  except  for  several  simple  examples  discussed  in  Section 
For  simplicity  we  will  again  take  n^,  and  i>Q  as  constants,  but  here  we 
will  include  the  nonlinear  terms  in  the  equations. 

The  equations  employed  are  Eqs.  ( lU)  to  (22).  We  know  that  the 
equations  are  hyperbolic  in  nature  and  therefore  that  solutions  may  be 
obtained  by  marching  forward  in  time.  This  is  shown  by  combining  Eqs.  (l) 


to  (?)  to  obtain: 


Vlv-t)  -V't  *  . 


In  the  cylindrically  symmetric  case  at  hand,  for  one  ion  species  present, 
this  equation  takes  the  form  (for  the  azimuthal  component): 


HTt€ 


'  (oj) 


The  treatment  of  the  equations  numerically  is  considerably  different, 
depending  on  whether  or  not  the  displacement  current  term  in  Eq.  (67)  is 
neglected.  To  gain  insight  into  this  question,  we  refer  to  the  linearized 
approximations  studied  in  the  preceding  section.  Requiring  all  the  dependent 
variables  to  -vary  as  exp(iKr  -  and  taking  the  limit  of  large  Kr,  we 

find  that  the  dispersion  relations  for  the  three  models  considered  are 
given,  respectively,  by  Eqs.  {l\Q) ,  (50)  and  (52). 
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From  Eq.  ('  w)»  we  see  that  if  we  use  Model  'V* ,  oscillatory  solution? 

are  possible  for  u)  a)  ^  uip  +  u)euJp,  and  in  a  narrow  band 

uJp  -  tdfUJp  ~  ic  uip  *  iif  .  (We  have  used  the  fart  that  u  i  /«^.  «  1.) 

On  the  other  liand,  if  we  use  the  dispersion  relations  (50)  and  (5i  )> 


again  neglecting  collisions,  we  obtain 


1* 


(69) 


We  find  oscillatory  solutions  for  u)  ~  and  u)  ~  uip*  If  the  displacement 
current  is  neglected  ("l"  in  Eqs .  (68)  and  (69),  all  solutions  are  expo¬ 
nentially  decaying  for  oj  >  in  Models  (b)  and  (c)  (electrically'  neutral 
plasma),  although  there  remains  a  narrow  band  of  oscillatory  solutions  near 
uj  -  alp  in  Model  (a) . 

Thus  if  ve  employ  tl«e  assumption  of  electrical  neutrality  and  neglect 
the  displacement  current,  only  the  frequencies  OJp  are  significantly 
misrepresented.  Since,  as  we  have  shown  in  Section  ,  phenomena  of  interest 
occur  at  characteristic  frequencies  far  below  we  may  neglect  the  displace¬ 
ment  current  with  very  little  error.  For  example,  for  co  «  (68)  or  (69) 

yields  n2  *  6  x  id  for  ng  **  10'  ,/cm^  and  B  =  0.3  gauss  (typical  values  in 
the  ionosphere). 

It  proves  very  convenient  to  omit  the  displacement  current  in  performing 
numerical  computations.  If  we  do  sc,  we  cannot  arbitrarily  specify  the 
electric  field  at  some  initial  time,  but  rather  we  must  calculate  it  from 


(70) 


where,  as  in  Section  V-  ,  v  -  v  -  v^.  Set 
equations  are  obtained  by  the  arguments  use 


-  v  .  Jetting  u.  ~  u  •  u,  the  remaining 
e  1  e 


d  to  obtain  Eq.  ( lli  )  •  wp  -  irld 


<5/  "  AV\' 


Combining  Eqs.  ( lr»)  and  (£0),  and  approximating  the  rentrifugai  term  by 


u.'&jl  _  ***>  v — C  **• "  w  •  ( 72 ) 

'Tt  *  ~  ^  >  ”  ^c 

To  close  the  system,  we  use  Eqs.  ( la )  and  (ltt). 
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Sample  Calculation 


The  procedure  outlined  in  the  preceding  section  has  been  incorporated 
into  a  n'imerieal  program  for  the  evaluation  of  the  current  distribution  and 
plasma  motion  induced  by  a  cylindrical  neutral  wind.  The  principal  objec¬ 
tive  is  to  be  able  to  study  the  non-linear  motion  of  the  plasma  during  the 
early  stages  of  the  expansion  of  a  dense  cloud  of  gas  when  the  collision 
frequencies  are  moderately  large.  We  expect  the  maximum  field  disturbances 
to  occur  when  the  neutral  ion  collision  frequencies  are  comparable  to  the 
n  cyclotron  frequencies.  The  calculation  is  set  up  to  take  a  neutral 
wind  calculated  in  the  small  disturbance  approximation  and  will  evaluate 
the  plasma  motion  in  a  parallel  magnetic  field  under  the  influence  of  this 
driving  force.  In  attempting  to  treat  the  neutral  flow  field  corresponding 
to  a  vacuum  expansion,  difficulties  were  experienced  because  of  the  severe 
ion  density  gradients  that  result  at  the  edge  of  the  sweeping  region  in  this 
situation  (see  Fig.  III-1*). 

This  problem  has  not  been  adequately  resolved  to  date.  A  simpler  but 
rather  hypothetical  situation  has  been  treated.  Here  the  neutral  flow 
velocity  was  chosen  to  have  a  form  similar  to  that  obtained  in  the  colli¬ 
sionless  representation  but  the  ion-neutral  collision  frequency  was  taken 
to  be  constant.  Specifically  the  neutral  radial  velocity  was  assumed  to 
have  the  form 


A(r/R 

V  =  (r/t)/^l  +  e 


1) 


j 


where  R  =  V  t  and  A  is  typically  a  number  between  5  and  10.  Typical  results 
o  o 

are  shown  in  Figures  V-3  to  V-6.  In  this  example  the  calculation  was  started 
after  the  neutral  flow  had  reached  a  certain  size  (i.e.,  the  interaction 
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between  the  ions  and  the  neutrals  was  switched  on  suddenly  alter  the 
neutral  flow  had  developed).  The  current  distribution  and  elec.,.!  and 
magnetic  fields  undergo  a  transient  stage  as  the  quasi-steady  solution  is 
approached.  Th4 s  sudden  transient  motion  generates  an  Alfv^n  wave  which 
propagates  away  from  the  source  of  the  disturbance  as  a  compression  pulse 
in  the  magnetic  field.  (See  Figure  V-3. )  Of  course  this  pulse  is  solely 
due  to  the  switching  on  of  the  interaction  and  has  no  real  physical 


counterpart . 


r,°I  itionsalp  of  F resent  Model  to  General  Ohm's  Law 

Neglecting  pressure  gradients,  we  may  write  the  electron  anu  ion 
momentum  equations  (6)  arid  (7)>  employing  (1C),  as  follows; 


In  the  steady  uniform  state,  the  left-hand  sides  of  the  above  equations 
vanish.  It  is  simple  to  :  ow  that  the  remaining  terms  lead  to  the 
customary  stead;  ta^e  Ohm's  Lav  for  a  weakly  ionized  gas,  including 
Hall  effect  and  ion  slip.  To  do  this,  ve  introduce  the  notation; 

'TH  U  *  4H€'  -  +  4*4  i 

tr  =  Vi  -  V€ 

-  **i  *  **€ 

<mv,  *  **iK  v« 

=  4*i  %  +  **•  V4 

Equations  (75)  and  (76)  may  be  written  (with  d/dt  =  0): 


c  t r«_a  f  mi>)  ( %  -u)  +  fa-V  ■>***;  v-  *Q 

Cm 


*  «!  '**€**1  C. 

Noting  that  J  =  en  v/c.  and  substituting  for  tT  in  (79)  from  (78) 


as  a  function  of  V  and  V^,  we  obtain 


1  - 1  mza&m — (m+  v*  **  ]-/•  Aj*t+  &$*$)*  b_ 

V  *<rrr3  '4^7* 


? s 


***Cm 


Z(Vt.-K)  _  *** 
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(75) 


(76) 


(77) 


(78) 


(79) 

(80) 


(90) 


where 


This  is  completely  equivalent  tr  the  well-known  result  (e.g..  Ref.  11, 
appendix),  where  the  result  is  given  as 


$  *  <r[a  +  K/*8  -  (9J.) 

and,  to  dominant  terms  (and  using  our  own  notation), 


CT  a  *X/*< 


i 


The  foxlowing  useful  expressions  for  the  particle  drift  velocities 

in  the  steady  uniform  state  are  obtained  from  (75)  and  (76)  in  the 

same  manner,  with  subscripts  j|  and  J,  inferring  to  components  parallel 

— * 

and  perpendicula'  tc  Bq: 


) 


(92) 


(93) 


It  is  usually  convenient  to  rewrite  Squati.on  (9l)  in  terms  of  the  Hall 
and  Pedersen  conductivities,  and  o^ .  Here 
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where  th'1  conductivities  are  given  by 


<r  (  >+  r\T3x)  _ 

(<  3-XS*!\,'+(j<.<r8)v 

‘  X  A> 

(u  r\&)  +  (kgtb) 


C“  *■  • 

\\ 


V - 5 .  Nume r ic al  Values  oi'  Cross  Sections  and  Collision  Frequencies 

We  consider  first  the  encounters  between  ions  and  neutrals  *  Following 
Dalgarno,^^  the  interaction  potential  at  large  separations  is 


4>(r)  *  - 


and  the  corresponding  momentum -exchange  cross  section,  defined  by  Eq.  (9)> 


c r  * 

m 


I 

2,2/ 1 r  V*) 


where  a  is  the  polarizability  of  the  neutral  atom,  p  =  mima(mi+  ma)  1  is 
the  reduced  mass,  and  V  is  the  relative  speed  of  the  ion  and  neutral 
particle.  The  collision  frequercy  is  thus  independent  of  temperature  or 


relative  velocity: 

*  2.2/ir  (<**//*•)  • 

At  the  lower  altitudes  of  interest,  say  150  km,  the  ions  and  neutrals 
are  mostly  diatomic. ^  Taking  them  to  be  N2,  Og,  N*,  02  and  NO  ,  the 
reduced  mass  and  polarizability  vary  only  slightly  between  the  various 
possible  ion-neutral  pairs,  and  we  find 


[diatomic  moleci  les 
and  ions] 


where  n  is  now  the  total  mrrber  density  of  neutrals,  and  we  have  used 
a 

-24  3 

a  =  1.5  x  10  cm  . 

At  the  higher  altitudes  of  interest,  say  300  km,  the  ions  and  neutrals 
are  mostly  monatomic N,  0,  N  ,  0  (mostly  0  and  0  ).  Again  p  and 
a  are  very  insensitive  to  the  choice  of  pairs,  and  we  find 


JA  ~  S,(o 


[monatomic  molecules 
and  ions] 


( 


Since  V  is  in  i  he  range  1  to  >  km/ sec,  we  are  dealing  with  cross  sect  ions 

of  the  order  of  2  to  o  x  10  ^ '  cm’  . 

There  is  in  addition  an  effective  momentum  exchange  between  ions  and 
neutrals  brought  about  by  charge  exchange.  This  is  much  more  complicated, 
inasmuch  as  there  is  often  an  exchange  between  unlike  particles.  However, 
exchange  between  unlike  particles  appears  to  be  unimportant  for  the  following 
reasons. 

(8' 

Following  the  discussion  by  Hasted'-  is  a  maximum  when  the  relative 

velocity  corresponds  roughly  to  resonance; 

^  (  ) 
where  the  length  parameter  a  *  7  x  10  cm  and  AE  is  the  energy  difference 
between  the  two  ionization  states.  For  V  «  Vves>  ac  obeys  approximately 
the  relation 

Cg(4*Ak)  (-A  )*5l/4&  V)  (7) 

-14  2 

where  r  (max)  is  usually  of  the  order  of  10  cm  . 
c 

For  the  species  N,  0,  0o,  NO,  there  are  ten  possible  different 

C.  C- 

charge -exchange  reactions  involving  unlike  particles  (each  reaction  is  two- 
way).  The  largest  | AE |  is  ueV,  for  (NO^,  N^)  or  (NO,  N*),  and  the  smallest 
| AE |  is  0.6?  eV,  for  (o'*,  N^)  or  (0,  II  ),  Thus  the  smallest  resonance- 

n 

velocity  V  is  that  for  O.oy  eV,  whicn  is  1.2  x  10  cm/sec.  Wr‘  consider 
only  velocities  of  at  most  several  km/sec,  which  by  (j)  lead  to  <r  several 
orders  of  magnitude  lower  than  c^,  for  unlike  particles. 

For  charge -exchange  interaction  between  an  ion  and  a  parent  molecule, 

V  =0.  The  only  such  interaction  we  need  consider  is  (0,  G),  since  they 
are  the  only  like-pair  in  which  the  ion  and  molecule  ure  both  dominant  in 
the  total  nvmiber  density  of  ions  and  molecules  respectively,  this  occurring 
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above  roughly  l8o  km.^  ^  The  experimental  measurements  of  St  ebb ings, 

(9) 

Smith  and  Erhardt'’  '  indicate  that  the  following  empirical  formula  holds 
from  about  0.1  eV  to  10'  eV: 

<£  *  ^ *  to  * -  0  63  *to  *  At,  e]  G**t 

where  £  is  the  energy  of  the  primary  particle  in  eV.  For  V  =  10^  cm/sec 

jc  p 

(0.08  eV),  o  ^  x  10  ‘  cm  (to  within  probably  a  factor  of  3)« 

For  the  sample  calculations,  which  are  performed  for  the  vicinity  of 
200  km,  we  therefore  take 

)/•  &  mO 

-/ 


(8) 


(9) 


Elec iron -neutral  momentum  exchange  cross  sections  may  be  obtained  from  the 

tables  of  Ref.  5  which  indicate  that  the  cross  sections  at  relevant  energies 

-15  2  l/° 

are  all  of  the  order  10  J  cm  .  Thus  v  is  of  the  order  (m  /m  )  '  v.  for 

e  '  i'  e'  i 

equilibrium  and  ~  for  cold  electrons  and  ions  moving  relative  to  the 

neutrals.  The  same  result  is  obtained  if  one  employs  the  model  of  Moiseivitsch10, 
which  is  equivalent  to  Eqs.  (l)  to  (3)  for  the  low  energies  under  considera¬ 
tion.  Again  we  obtain,  independently  of  considerations  of  temperature  or 
re  lat  ive  ve  loc  it  ie  s : 


A 


**  *  ('mi  /T.  )  Vj- 


(10) 


ter 


or 
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